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frequency  conversion  in  the  infrared.  Detailed  studies  have  shown  excellent 
theoretical  and  experimental  agreement  in  the  Raman  interactions. 
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ABSTRACT 


A program  of  basic  research  in  laser  physics  and  laser  techniques 
has  been  carried  out  in  order  to  help  meet  Air  Force  mid-term  and  long- 
term requirements  as  outlined  in  the  AFSC  Research  Planning  Guide.  A 
Fast  Hankel  transform  algorithm  has  been  developed  and  applied  to  non- 
linear interactions  and  to  diffraction  biased  resonators  studies.  Pico- 
second spectroscopy  via  the  induced  grating  approach  has  been  successfully 
used  to  study  fast  relaxation  processes.  Optical  parametric  amplifiers 
and  oscillators  have  been  shown  to  be  useful  sources  of  coherent  radiation 
over  the  0.4  ym  to  4 Um  spectral  range.  Finally,  stimulated  Raman  scat- 
tering and  four-wave  mixing  has  been  used  for  efficient  frequency  conver- 
sion in  the  infrared.  Detailed  studies  have  shown  excellent  theoretical 
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REVIEW  OF  ACCOMPLISHMENTS 


I.  INTRODUCTION 

This  annual  report  covers  research  work  carried  out  by  Professor 
A.  E.  Siegman  in  the  areas  of  optical  beam  and  resonator  calculations 
and  by  Professor  R.  L.  Byer  in  the  areas  of  optical  parametric  amplifier 
and  oscillator  studies  and  stimulated  Raman  scattering  studies. 

During  the  past  year  several  important  results  have  been  obtained 
under  this  research  program.  In  particular,  advances  have  been  made  in 
the  application  of  the  Fast  Hankel  transform  algorithm  to  the  problems  of 
nonlinear  interactions  and  diffraction  biasing  to  prevent  locking  of  ring 
gyroscopes.  Picosecond  spectroscopy  techniques  have  been  extended  to  the 
laser  induced  grating  technique  and  to  photo-acoustic  spectroscopy. 

Detailed  theoretical  and  experimental  results  were  obtained  for  both  opti- 
cal parametric  amplifiers  and  optical  parametric  oscillators.  Finally, 
stimulated  Raman  scattering  was  studied  and  applied  to  problems  of  effi- 
cient frequency  conversion  in  the  infrared  spectral  range. 

These  studies  have  led  to  publications  which  are  listed  below  and  to 
the  completion  of  four  Ph.D.  theses. 

II.  OPTICAL  BEAM  AND  RESONATOR  CALCULATIONS 

During  the  past  year  very  substantial  steps  were  accomplished  in  com- 
pleting and  testing  the  novel  Fast  Hankel  transform  algorithm  that  we  first 
Invented,*  and  in  applying  it  to  optical  calculations,  particularly  to 
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optical  resonator  and  optical  beam  propagation  calculations  of  importance 
to  high  power  laser  applications.  Extensions  to  the  algorithm  itself 
included  a simple  "lower  end  correction"  technique  to  compensate  for  the 
finite  truncation  of  the  transform  integral  at  the  lower  end  of  the  inte- 
gration range;  more  accurate  ways  of  handling  the  discrete  Fast  Fourier 
transform  steps  involved  in  the  calculation;  and  better  understanding  of 
the  choice  of  sampling  intervals  in  the  calculation.  These  improvements 
brought  the  accuracy  of  the  transform  close  to  the  limits  of  single  pre- 
cision arithmetic  on  our  machine  (fractional  error  ~ 10  ^).  In  addition, 
a two-stage  FHT  algorithm  requiring  roughly  twice  the  storage  and  compu- 
tation time  was  developed  which  reduces  the  error  of  the  algorithm  itself, 

-12 

within  its  range  of  validity,  to  below  ~ 1’’  which  is  probably  better 
than  will  ever  be  needed. 

In  addition  substantial  understanding  was  developed  and  testing 

carried  out  on  the  ranges  of  applicability  of  the  transform,  especially 

in  optical  beam  calculations,  in  order  to  overcome  the  brute-force,  cut- 

and-try  methods  too  often  used  in  these  calculations.  Some  of  these  re- 

2 

suits  were  reported  at  the  1978  Annual  Meeting  of  the  OSA  and  a major 

publication  summarizing  all  of  this  work  to  date  will  be  completed  and 

3 

submitted  for  publication  shortly. 

Work  on  the  application  of  fast  transform  methods  to  nonlinear  beam 
4 

calculations  was  also  carried  to  completion.  These  results  will  be  pub- 
lished in  the  near  future."* 

There  is  still  useful  work  to  be  done  on  this  topic.  In  particular, 
with  the  FHT  algorithm  itself  now  essentially  completed  it  would  be  very 
useful  to  test  and  verify  the  hueristic  "hand  waving"  criteria  generally 
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used  to  select  number  of  points  and  sampling  size  in  optical  diffraction 
calculations,  so  as  to  define  the  best  compromise  between  accuracy  and 
minimum  number  of  points.  The  trade-off  between  "two-step"  and  "one-step" 
transform  propagation  methods  needs  to  be  tested,  especially  the  alleged 
sensitivity  to  under-sampling  of  the  latter.  Many  useful  special  resonator 
types,  such  as  hole  coupled  and  soft  edged  unstable  resonators,  need  to  be 
studied.  We  are  continuing  to  work  on  all  of  these  topics. 

Finally,  a rather  different  and  unusual  technique  of  "diffraction 
biasing"  to  obtain  non-reciprocal  propagation  and  freedom  from  locking  in 
ring  laser  gyroscopes  was  invented  and  is  being  submitted  for  publication. 
Whether  this  technique  will  be  of  practical  importance  remains  to  be  seen, 
but  it  does  represent  a quite  different  and  novel  attack  on  a long  standing 
problem  of  substantial  importance  to  Air  Force  requirements. 

III.  PICOSECOND  SPECTROSCOPY 

Some  very  interesting  measurements  on  the  picosecond  spectroscopy  of 
ultrafast  relaxation  times  in  solids  were  completed  during  this  year,  using 
the  novel  transient  laser-induced  grating  technique  developed  in  our  group. 
The  measurements  demonstrated  in  particular  a novel  photo-acoustic  pressure- 
induced  amplitude  grating  effect  in  molecular  crystals,  of  a type  which 

does  not  appear  to  have  ever  been  observed  previously.  This  work  was  re- 

7 8 

ported  at  the  1978  OSA  Meeting  and  a complete  publication  is  now  in  press. 

A preprint  of  this  article  is  attached  as  Appendix  I of  this  report. 

An  excellent  Ph.D.  dissertation  on  this  topic  was  completed  by  Dr.  Jose 

9 

Salcedo  at  the  end  of  this  year.  Dr.  Salcedo  will  be  staying  on  as  a post- 
doctoral fellow  to  provide  very  valuable  assistance  in  continuing  this 
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research,  at  no  cost  to  this  contract. 


IV.  OPTICAL  PARAMETRIC  OSCILLATOR  AND  AMPLIFIER  STUDIES 

During  this  past  year  two  extensive  theoretical  and  experimental 
studies  were  completed  on  optical  parametric  amplification^  and  optical 
parametric  oscillators.^ 

The  gain,  conversion  efficiency  and  tuning  behavior  of  KD*P  and  LiNbO^ 

OPA’s  were  studied  in  detail.  The  experimental  results  were  shown  to  be  in 

excellent  agreement  with  theory.  These  results  underscored  the  importance 

of  OPA's  as  efficient,  widely  tunable,  coherent  radiation  sources  across 

the  visible  and  near  infrared  spectral  regions.  The  work  was  part  of  the 

12 

Ph.D.  thesis  of  Dr.  Richard  Baumgartner  and  is  to  be  published  in  the 

June  issue  of  the  Journal  of  Quantum  Electronics.  A preprint  of  the  paper 

is  attached  as  Appendix  II.  This  work  presented  the  first  quantitative 

study  of  parametric  amplifiers  since  their  discovery  in  1965. 

A second  detailed  study  of  the  threshold  and  linewidth  of  a LiNbO^ 

optical  parametric  oscillator  was  completed  recently.  This  study,  which 

13 

fulfilled  part  of  the  Ph.D.  requirement  for  Dr.  S.  J.  Brosnan,  for  the 
first  time  compared  detailed  theoretical  models  for  0P0  operation  with 
careful  measurements.  The  results  were  submitted  for  publication  and  have 
been  accepted  as  an  invited  paper  for  the  June  issue  of  the  Journal  of 
Quantum  Electronics.^  The  threshold  and  conversion  efficiency  behavior 
of  an  0P0  is  now  well  characterized.  In  addition,  linewidth  studies  led 
to  successful  operation  of  the  0P0  at  8 mJ  per  pulse  at  10  pps  in  a stable 
single  axial  mode.  A preprint  of  this  paper  is  attached  as  Appendix  III. 
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STIMULATED  RAMAN  STUDIES 


During  the  past  year  we  completed  deta_  ed  studies  of  stimulated 
Raman  processes  for  conversion  of  tunable  radiation  sources  to  new  fre- 
quencies. The  work  considered  both  stimulated  emission  processes  and 
four-wave  mixing  interactions  and  the  close  coupling  of  the  two  phenomena. 

We  have  theoretically  and  experimentally  verified  the  pump  linewidth 

14  15 

independence  of  stimulated  Raman  scattering.  ’ This  work  was  prepared 
for  publication  and  is  to  appear  in  the  July  issue  of  the  Journal  of  Quantum 
Electronics.  The  theoretical  study  showed  in  an  elegant  manner  via  coupled 
mode  analysis  the  important  role  four-wave  mixing  plays  in  the  stimulated 
Raman  scattering  process. 

The  threshold/bandwidth  studies  were  undertaken  in  support  of  infra- 
red rotational  stimulated  Raman  studies  in  hydrogen  gas  at  10.6  yin.  In  a 
very  important  experiment,  we  were  able  to  achieve  40%  photon  conversion 
efficiency  from  10.6  ym  to  16.95  ym  by  stimulated  Raman  scattering. ^ ^ ^ 
This  experiment  demonstrated,  for  the  first  time,  that  efficient  conversion 
of  high  energy  and  peak  power  infrared  sources  was  possible  by  the  stimu- 
lated Raman  scattering  approach.  The  results  of  the  rotational  Raman  ex- 
periment have  been  published  and  are  attached  as  Appendix  IV. 

It  should  be  noted  that  for  the  first  time  it  is  possible  to  conceive 
of  arbitrarily  large  energies  available  over  the  8 ym  to  20  ym  spectral 
region  by  rotational  Raman  scattering  of  a pressure  tunable  CO^  laser  source 
in  H^,  D^  and  HD.  Thus,  Raman  scattering  should  play  an  increasingly 
important  role  in  the  generation  of  tunable  high  average  power  infrared 
radiation . 
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In  a separate  study  we  quantitatively  compared  in  detail  four-wave 

mixing  theory  and  experiment  for  application  to  a widely  tunable  infrared 

13 

source.  This  work  was  also  part  of  Dr.  S.  J.  Brosnan's  Ph.D.  thesis. 

The  thesis  is  now  being  prepared  as  a Ginzton  Labor;  ory  report  and  as 

a technical  report  for  this  program.  The  four-wave  mixing  studies  showed 

that  continuously  tunable  coherent  radiation  is  now  available  over  the 

extended  4 - 18  ym  spectral  range  from  a LiNbO^  0P0  pumped  mixing  cell. 

19 

The  results  are  being  prepared  for  publication. 

A key  technical  innovation  of  the  Raman  studies  was  the  invention  and 

20 

demonstration  of  a multiple  pass  Raman  cell.  The  cell  allowed  signifi- 
cant reduction  in  Raman  pump  power  by  providing  25  passes  in  a re-focusing 
geometry.  The  technical  paper  describing  the  multiple  pass  cell  is  being 
prepared  for  publication.  The  cell  is  also  the  subject  of  a patent  appli- 
cation under  Air  Force  Office  of  Scientific  Research  support. 
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I.  R.  SALCEDO  and  A.  E.  SICCMAN.  *lxu>w.  if>:k 


Abttrmcl-k  liMlient  piling  of  tingle!  electronic  excited  itilci  il 
produced  k»  • penticenc-dopaJ  p-terphcny|.  moWculir  crystal  by 
optical  absorption  from  two  crossed  time-coincident  picosecond  ex* 
dtitioa  pokes  it  532  nm.  The  diffraction  properties  of  this  volume 
piling  ire  probed  by  I weak,  variably  delayed.  Bragg-male  bed  pico- 
second probe  pulse.  At  high  excitation  intensities,  a strong  oscillatory 
behavior  in  the  time-dependent  scattering  efficiency  is  observed  super- 
imposed on  the  exponential  decay  pattern  of  the  excited-state  paling. 
We  attribute  the  oscillatory  behavior  to  a thermal  paling  which  induces 
Coherent  microwave  acoustic  phonons.  These  in  turn  modulate,  at  the 
sound  frequency,  the  optical  absorption  properties  of  the  pcntaccne 
molecules  in  the  excited  state.  An  acoustically  induerd  •mphtude 
grating  effect  is  thus  obtained,  in  contrast  to  conventional  acousto- 
optic phase  paling  effects. 


f.  Introduction 

WE  report  the  observation  of  strong  photoacoustic  ef- 
fects in  molecular  crystals,  and  particularly  the  observa- 
tion of  an  acoustooptic  amplitude  grating  effect,  obtained 
while  applying  the  picosecond  transient  grating  method  to  the 
study  of  energy  migration  processes  in  molecular  crystals. 

The  details  of  our  transient-grating  method  and  of  the  en- 
ergy migration  studies  have  been  given  elsewhere  [1],  (2), but 
it  may  be  appropriate  to  briefly  review  the  transient  grating 
concept.  Two  coherent,  timc-coincidcnt  picosecond  excita- 
tion pulses  are  transmitted  through  an  cxpciitnental  sample 
which  is  absorbing  at  the  excitation  wavelength,  as  in  Fig.  I. 
In  the  experiment  described  here  the  absorbers  consist  of  or- 
ganic molecules  (pcntaccne)  imbedded  in  a transparent  crystal- 
line host  (p-terplicnyl).  Optical  absorption  in  the  interference 
pattern  between  the  two  excitation  beams  produces  a volume 
grating  pattern  of  electronic  singlet  excited  states  of  the 
pcntaccne  molecules.  This  cxcitcd-statc  grating  then  acts  as 
an  optical  absorption  grating  or  “transient  hologram.”  The 
time-dependent  decay  of  the  excited-state  grating  is  probed 
by  Bragg  diffraction  of  a weak,  variably  delayed  picosecond 
probe  pulse.  The  intensity  of  the  diffracted  probe  pulse  is 
monitored  as  a function  of  probe  pulse  delay.  In  the  ease  at 
hand,  the  grating  pattern  decays  through  a combination  of 
excited-state  decay  and  cxcitcd-statc  transport  (which  washes 
out  the  grating  fringes).  For  purely  diffusive  cxcitcd-statc 
energy  transport,  the  decay  is  exponential,  and  the 'decay  rate 
gives  the  transport  parameters  directly.  A diffusion  coef- 
ficient for  the  singlet  excited  electronic  states  of  pcntaccne 
in  p-terpheny!  has  thus  been  obtained  [2] . 
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When  we  performed  llic  same  experiments  at  optical  excita- 
tion intensities  above  ~ 500  MW/cm1,  a strong  well  defined 
pscillatoiy  behavior  in  the  scattering  efficiency  versus  decay 
was  observed.  In  this  paper  we  interpret  this  oscillatory  be* 
havior  as  due  to  strong  photoacoustic  effects.  At  high  excita- 
tion intensities,  a significant  thermal  grating  is  sii|«cr imposed 
on  the  excited-state  grating.  This  thermal  grating  produces  a 
transient  disturbance  which  can  be  described  as  a "frozen" 
stress  pattern  plus  two  countcrpropagaling  coherent  micro- 
wave  acoustic  waves  (3).  The  density  modulation  associated 
with  this  disturbance  modulates,  at  the  microwave  sound  fre- 
quency, predominantly  the  optical  absorption  cross  section  of 
the  pcntaccne  excited  electronic  singlet  slates,  most  likely  by 
pressure  tuning  the  absorption  line.  This  in  turn  modulates 
the  effective  depth  of  the  excited-state  grating  and  thus  its 
scattering  efficiency,  producing  an  amplitude  grating  effect. 
To  our  knowledge,  such  photoacoustic  amplitude  grating  ef- 
fects have  not  been  previously  observed.  They  provide  inter- 
esting information  on  phonon  processes  and  on  pressure  tun- 
ing in  the  crystals,  in  addition  to  preserving  the  excited-state 
transport  information. 

II.  Experimental  Ouservations 
A.  Apparatus 

The  experimental  apparatus  for  these  studies  is  outlined  in 
Fig.  2.  A continuously  pumped,  high  repetition  rate,  Q- 
switched  and  actively  mode  locked  NdrVAG  laser  (4)  pro- 
duces bursts  of  mode-locked  pulses,  from  which  single  pulses 
are  selected  by  a UNbOj  Pockcls  cell  in  a Blumlcin  configura- 
tion. These  pulses  are  efficiently  doubled  to  532  nm  by  a 
temperature-tuned  90°-phasc -matched  CD#A  crystal,  and  split 
into  two  strong  excitation  pulses  (^300  kW  peak  power  each) 
and  a weak  variably  delayed  probe  pulse  (fc20  KW).  All 
pulses  are  linearly  polarized  in  the  plane  of  the  figure,  and 
are  Gaussian  in  time  and  in  space  with  a duration  [full  width 
at  half  maximum  (FWllM)]  of  50-75  p$.  Dvic  to  the  par- 
ticular optical  absorption  properties  of  pcntaccnc  at  532  nm, 
excitation  and  probe  pulses  arc  chosen  to  be  at  the  same  wave- 
length. To  meet  the  Bragg  condition,  the  probe  pulse  is 
brought  from  the  back  of  the  crystal  along  the  reverse  path  of 
one  of  the  excitation  pulses,  but  with  its  polarization  rotated 
by  90°.  The  first-order  Bragg-diffracted  pulse  preserves  this 
polarization  and  is  scattered  along  the  reverse  path  of  the 
other  excitation  pulse.  A polarization  beam  splitter  then  de- 
flects the  signal  into  a photodctcctor  connected  to  a lock-in 
amplifier.  The  lock-in  output  drives  the  y axis  of  an  X-Y  re- 
corder, while  the*  axis  is  driven  by  Cite  motorized  probe  pulse 
delay  line.  For  slow  delay  scan  rate  and  high  laser  repetition 
rate  (100  Hz),  data  gathering  is  essentially  continuous  in  a stro- 
boscopic sense."  Each  experimental  plot  contains  typically 
50  000-100  000  laser  shots  and  takes  10  min  to  complete.  Ob- 
served scattering  efficiencies  arc  typically  several  percent,  y ield- 
ing high  signal  to  noise  ratio  in  the  measurements  (.S/A'  = 1 00  is 
easily  achieved).  To  Outlier  enhance  the  signal  discrimination, 
every  other  pair  of  excitation  pulses  is  chopped  and  phase 
sensitive  detection  is  used.  A microscope  objective  can  also 
be  used  to  image  the  interaction  volume  onto  a vidicon  tube 
with  taigc  magnification,  allowing  the  display  of  the  grating 
fringes  on  a television  monitor,  "lliis  provides  a very  accurate 
tool  to  align  all  three  beams  in  the  crystal. 
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B.  Results 

Results  showing  diffusive  energy  transport  at  low  excitation 
Intensity  have  been  previously  reported  and  studied  |2|.  Here 
we  concentrate  on  the  high  excitation  intensity  results,  which 
show  the  oscillatory  character  of  the  timc-dc|icndciit  scatter- 
ing efficiency.  Fig.  3 shows  this  oscillatory  dependence  as  ob- 
served in  four  experiments  from  a series  performed  under  iden- 
tical conditions  except  for  different  fringe  spaeings  A,  which 
are  adjusted  by  changing  the  angle  between  live  two  excitation 
beams.  In  all  these  results  the  grating  k vector  is  parallel  to 
the  crystal  b axis.  The  time  period  T of  the  oscillatory  modu- 
lation changes  with  fringe  spacing  A with  the  linear  depen- 
dence shown  in  Fig.  4.  A straight-line  fit  yields  a slope  of 
2.63  X 10’  ctit/s,  which  we  take  to  be  the  velocity  of  sound. 
The  velocity  of  sound  in  the  same  crystal  along  the  same  di- 
rection was  independently  measured  to  be  v,  = 2.65  t 0.05 
em/s.  confirming  the  acoustic  origin  of  the  oscillations. 

Further  examination  of  the  shape  of  the  oscillatory  decays 
In  Fig.  3 leads  us  to  conclude  that  the  observed  effects  must 
result  from  thermal  excitation  of  an  acoustooptic  amplitude 
grating.  This  grating  must  be  produced,  we  believe, by  acoustic 
pressure  tuning  of  the  excited  molecular  state  absorption  in 
the  sample.  The  theoretical  analysis  supporting  these  con- 
clusions is  outlined  in  the  following  section. 


III.  Tiikoklticai.  Formulation 
A.  rune  Dependence  of  the  Scattering  Efficiency 
The  geometrical  and  parametric  dependence  of  the  scattering 
efficiency  of  volume  liologiaius  is  well  known  |S],  (6|.  Here 
we  arc  mainly  interested  in  the  exponentially  decaying  and 
oscillating  time  dependence  of  the  observed  scattering  from 
the  volume  grating  patterns  produced  by  thermal  plus  excited- 
state  gratings  in  our  experiments.  In  these  experiments 
scattered  probe  signals  were  observed  primarily  at  the  first- 
order  Bragg -matching  angle.  However,  the  volume  grating 
pattern  also  contains  higher  spatial  harmonics,  particularly 
at  the  second-order  Bragg  angle,  which  should  also  be  observ- 
able. In  general  the  diffracted  probe  signal  intensity  /„(t)  at 
the  nth  order  Bragg  matching  angle  is  given  by 


'i„(r)=J  lp{t‘-  t)\l'(i')*X„(t')\'  dt'  (I)  ^ 

where  the  • denotes  convolution,  and  /p(t'  - t)  and  /e[t  ) arc 
the  probe  and  excitation  pulse  intensities,  with  the  probe 
pulse  being  delayed  in  time  by  an  amount  t.  Out  of  the  total 
complex  susceptibility  x(*.0  »»»  the  sample  flic  relevant  part  - 
is  the  spatial  Fourier  component  of  order  w,  i.c.,  the  compo- 
nent X«(0c*P  (M»).  which  corresponds  to  spatial  vector 
rk  = #i2»r/A,  where  A is  the  grating  spacing.  As  written  here. 
X(x.t)  is  the  Green's  function  or  the  unit  impulse  response 
produced  by  an  excitation  pulse  that  is  a unit  impulse  in  time. 

In  this  paper  we  will,  in  fact,  treat  the  excitation  pulse  /,(/) 
as  a unit  impulse,  since  the  *'•70  ps  excitation  (and  probe) 
pulses  arc  short  compared  to  the  **  I ns  oscillation  and  decay 
times  in  these  experiment*,  the  scattering  susceptibility  X»«(0 
is  in  grucial  complex,  i.c.,  x»» ~ x)(  4/X«*  *ilh  *hc  ic.il  pail 
being  responsible  for  phase  grating  or  index  of  refraction 
grating  effects,  and  the  imaginaiy  part  being  responsible  for 
amplitude  or  ahsoiplion  grating  effects. 

Hie  physical  process  that  produces  the  coherent  acoustic 
*avcs  of  interest  hcic  is  a thermally  induced  acoustic  grating 
which  is  superimposed  on  the  excited-stale  grating  |?|  The 
thermal  giating  and  the  associated  acoustic  waves  can  he 
analyzed  stalling  from  the  linearized  hydrodynamic  equations 
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Ikit  AT and  A p aic  llw  kxal  tcni|Wialiiie  and  density  change,. 
ft  Is  Ihc  c<|uilihfiuitt  density,  i)  is  the  vcl petty,  p is  llic  cipan- 
tion  coefficient,  X is  the  llterinal  conductivity,  a is  the  optical 
absorption  coefficient,  and  is  the  local  optical  electric  field 
in  the  interaction  volume.  Elcclrostriction  can  be  neglected, 
and  the  transient  solution  for  the  density  variation  A p can 
then  be  readily  obtained  (8)  for  an  cacitalion  pulse  very  short 
compared  to  the  acoustic  limes  of  interest: 


v/fce>s  t'-f-y 


Ap(jr.  f)  « A pm„  [e',/T"  cos  cut  - 


cos  (2»x/A) 


0) 

Here  rse  b Ihc  acoustic  wave  attenuation  time  or  sound  absorp- 
tion lime,  rth  is  the  thermal  conduction  time,  and  <o  * 2 nJrK 
b the  acoustic  frequency  where  ta  * A/e,,  with  u,  being  the 
acoustic  wave  velocity.  For  Ihc  acoustic  frequencies  generated 
in  these  experiments  (<o/2n  300-800  MHz),  the  acoustic 

decay  time  b r#c  ^-50-100  ns,  and  the  thermal  conduction 
decay  time  can  be  ignored. 

A more  complete  analysis  would  take  into  account  the 
Gaussian  transverse  spatial  dependence  of  the  excitation 
beams,  as  well  as  the  finite  duration  of  the  excitation  pulse 
/,(/).  In  particular  there  will  be  an  additional  effective  decay 
of  the  acoustic  signals  as  the  finite  trains  of  left-  and  right- 
going acoustic  cycles  move  across  each  other  and  eventually 
cease  to  overlap  spatially  with  each  other.  Since  our  gratings 
typically  had  at  least  20  fringes  across  the  half  width  of  (!>c 
beam  spot  size,  however,  this  had  little  effect  during  the  5 or 
6 acoustic  cycles  typically  followed  in  our  experiments. 

Although  we  will  speak  loosely  of  pressure  tuning  effects, 
we  attribute  the  observed  effects  in  our  experiments  to  in- 
duced changes  in  the  local  optical  properties  of  the  medium 
caused  primarily  by  local  density  changes  rather  than  either 
local  temperature,  local  stress,  or  pressure  effects.  Fig.  5 is 
a plot  of  the  excited  state  grating  pattern  Nx(x%t-  0),  the 
ground  state  molecular  density  Af0 (x,te  0),  and  the  instanta- 
neous local  density  change  Ap(.v,  t),  with  the  latter  quantity 
being  plotted  at  successive  instants  of  time  through  one  com- 
plete acoustic  cycle,  neglecting  damping  or  decay  effects.  The 
acoustic  density  grating  pattern  is  in  phase  spatially  with  the 
excited-state  grating  pattern,  as  well  as  with  Ihc  sinusoidal  part 
of  the  ground  state  molecular  density.  Because  of  the  instan- 
taneous nature  of  Ihc  thermal  excitation  in  bulk,  both  the  den- 
sity and  the  pressure  must  be  described  as  having  a static  or 
•‘frozen**  grating  component  as  well  as  two  countcrpropagating 
traveling-wave  components.  (The  “frozen”  components  will  of 
course  decay  away  w ith  a much  longer  thermal  relaxation  time.) 
Note  also  that  the  density  variation  with  time  at  any  one  point 
In  spact  arc  “one-sided,”  i.c.,  ihc  local  density  swings  up  and 
back,  or  down  and  back,  but  not  both  up  and  down  at  any  one 
location.  The  molecular  densities  A'0(v)  and  A'i(x)  decay 
slowly  due  to  cxcitcd-statc  relaxation  and  diffusion. 

B.  Phase  Versus  Amplitude  Crating  effects 

Conventional  acoustooptic  interactions  in  solids  nearly  al- 
ways occur  through  a phase  grating  mechanism,  caused  by  a 
change  in  local  index  of  refraction  associated  with  the  acoustic 
disturbance.  Examination  of  the  theoretical  expressions  and 
of  Fig.  5 shows,  however,  that  this  cannot  he  the  primary  ex- 
planation for  (Mir  observations,  and  that  we  arc  dearly  observ- 
ing an  acouslooptic  amplitude  grating  effect. 

In  our  results  the  instantaneous  scattering  efficiency  is  seen 
to  be  reduced  nearly  to  zero  at  each  odd  half  acoustic  cycle 
after  Ihc  excitation  pulse.  This  is  also  the  instant  when  the 
density  grating  teaches  its  maximum  amplitude  in  spatial 
phase  with  the  cxcitcd-statc  grating.  The  implication  is  that 
the  density-induced  scattering  at  those  instants  must  add  in 
opposite  phase  (not  in  quadrature)  to  Ihc  cxcitcd-statc  in-, 
duccd  grating.  We  believe  the  excited  slate  grating  must 
lutcly  be  an  amplitude  grating,  and  hence  the  density -induced 
giatingiiuist  be  an  amplitude  rather  than  a phase  grating  aho. 
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Equation  (3)  allows  that  the  density  grating  consists  of  two 
traveling-wave  components  of  relative  amplitude  ♦ J each,  and 
a fixed  or  constant  component  of  relative  amplitude  - 1.  The 
two  traveling-wave  components  decay  with  the  acoustic  decay 
time  t0C.  which  lasts  25-50  acoustic  cycles,  while  the  fixed 
component  decays  with  the  thcimal  decay  time  rth,  which  is 
much  longer.  Each  of  these  components  will  induced  a 
scattered  signal  field  component  at  the  Bragg  angle,  with  cor* 
responding  amplitudes  and  phases.  Note  that  the  three  den* 
sity  grating  components,  and  also  the  three  vectorial  scattered 
wave  components,  sum  to  exactly  0 at  / = 0 since  the  density 
modulation  starts  from  an  initial  value  that  is  identically  zero. 

FiE.  6 shows  vcctorially  in  the  complex  plane  the  scattered 
vector  field  amplitudes  contributed  by  each  of  these  three  den- 
sity wave  components,  plus  the  larger  component  contributed 
by  (he  cxcitcd-statc  grating,  at  three  instants  of  time  corre- 
sponding to  0,  } , and  I acoustic  cycle.  The  sketches  arc 
drawn  with  vertical  arrows  representing  the  excited  state  grat- 
ing. assumed  to  be  a pure  amplitude  grating.  The  heavy  arrow 
in  each  case  b the  vector  sum,  or  the  instantaneous  vector 
amplitude  of  the  total  Bragg-scattcrcd  signal  field.  The  length 
of  this  vector  squared  would  give  the  total  Bragg-scattered  in- 
tensity at  that  instant. 

The  upper  sketch  assumes  that  the  density  waves  induce  a 
phase  grating  effect,  i.c.,  the  three  acoustically  induced  com- 
ponents are  90*  out  of  phase  with  the  exCited-state  grating 
component;  while  the  lower  sketch  assumes  the  density  wave 
components  are  in  phase  with  the  excited-state  component. 
The  two  smaller  density  wave  components  rotate  in  opposite 
directions  in  the  complex  plane  at  the  acoustic  frequency, 
since  they  represent  oscillating  or  traveling-wave  terms.  Note 
that  the  observed  amplitude  variation  of  the  total  scattered 
signal  with  time  can  only  be  explained  by  the  lower  sketch. 
An  acoustic  phase  grating  would,  in  fact,  cause  the  total 
scattered  intensity  to  have  an  amplitude  variation  of  opposite 
sign. 

C.  Detailed  Analysis  of  the  Scattering  /efficiency 

Leaving  aside  inessential  constants  and  geometrical  factors, 
the  scattering  susceptibility  Xn(0  can  be  written  as 

X«(0*/ »/.#(«.  oj  (4) 

where  Nlfi  and  o,  p are  the  density-modulated  molecular  state 
populations  and  cross  sections,  respectively,  and  where  the 
subscript  n indicates  that  only  the  cxp(/>i2irjc/A)  spatial  com- 
ponent is  to  be  considered.  The  sum  in  the  present  ease  ex- 
tends over  two  levels,  the  ground  state  i ~ 0 and  the  first  ex- 
cited singlet  level  # =*  1 . The  population  and  cross  sections  may 
be  written  as 

♦ Ap(x./)/p0) 

o i',(x,l)  * o,-  ♦ (da, /Op)  Ap(r,  /)  (5) 

where  A',(v./)  and  </,  without  arguments  give  the  unperturbed 
values.  The  values  both  of  cr,  and  of  their  derivatives  may  in 
general  be  complex,  with  the  ical  part  corresponding  to  ab- 
sorption and  the  imaginary  part  corresponding  to  phase  shift 
or  index  of  refraction  effects. 

3hc  gcncial  expression  for  Ap(v,f)  has  been  given  in  (3). 
The  populations  of  (he  ground  and  excited  levels  produced 
by  the  excitation  beam  delta  functions  arc  |2) 

A'lO.O  - } AT"  <•-'"(!  - c'1’"'  co$(2<i.«/A)| 
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where  Nq  without  argument's  is  the  total  density  of  absorbing 
molecules  in  the  crystal,  A, is  the  initial  optically  induced 
excited-stale . density  at  the  grating  peaks,  r is  the  excited 
State  fluorescent  lifetime,  and  (l/r  ♦ k*D)  - K/2  is  the  decay 
rate  for  the  cxcitcd-statc  fringe  pattern  due  to  excited-state 
relaxation  plus  excited-state  diffusion  with  diffusion  coef- 
ficient D. 

We  now  define  the  normalized  maximum  density  variation 

U-APm.,/Po  «|  (7) 

and  Ihe  oscillatory  modulation  component 

*1(1)  * «P ( l/r.c)  cos  <j/  - I . (8) 

We  also  define  the  cross  section  sensitivity  parameters,  or 
normalized  density  derivatives 

5,  “ (3o0/3p)p0/o0 

S,  “(do, /dp)po/o,  (9) 

which  could  be  complex  or  even  imaginary  if  the  density 
modulation  produces  an  index  of  refraction  grating.  We 
then  obtain  the  n - I Bragg  angle  component  of  the  scatter- 
ing susceptibility  from  (4),  (8),  and  (9)  as 

-/X.(')  = (/V.m,72)(o,  - o0)e-<|/'“,«' 

+ 00,(5,  +l)W"572)e-"’d/(r). 

♦ 900(5,  + I)  (/V,  - (/V,m,,/2) e"/')y lf(r).  (10) 

This  expression  contains  three  terms.  The  first  is  tire  basic 
excited-state  transient  grating  signal  (2J  unperturbed  by 
acoustic  effects.  It  gives  the  pure  exponential  decay  observed 
at  low  intensities,  as  determined  by  the  excited-state  lifetime 
r and  the  diffusion  coefficient  0.  The  second  and  third  terms 
describe  how  the  first  excited  state  and  ground-state  diffrac- 
tion effects  are  modulated  by  the  thermally  induced  density 
variations  or  coherent  acoustic  waves.  In  the  (5|4  I)  factors 
In  the  second  and  third  terms,  the  factor  of  I represents  simple 
density  modulation  of  the  local  molecular  density,  while  the 
Si  factors  represent  density  modulation  of  the  molecular  cross 
sections  oi , as  in  (5).  It  will  turn  out  that  the  magnitudes  of 
the  Si  seem  to  be  substantially  larger  than  unity  to  fit  the  ob- 
served results. 

The  resulting  diffracted  signal  intensity  is  then  given  by 


/,(/)>/  /(,(»'-  OlfW 

+/(,«*'>  *f(<') 

♦ ^on/Vo/A'!""  - e-''l’)*l(l')  |‘  Jr'  (II) 


where  we  have  normalized  the  initial  scattering  efficiency  to 
unity.  The  oscillatory  terms  have  magnitudes 


t|(5p  ♦ I)Qq 

Aq  m 

0 1 - 00 

I)(5,  + Do, 

A | c 1 

0,  “ 0O 


(12) 


Note  that  squaring  this  cxpiession  not  only  squares  each  of  the 
•hove  terms,  hut  most  impoi lairtly  leads  to  heating  effects 
between  the  basic  excited-state  grating  and  the  two  acoustic 
terms,  because  the  probe  pol.cwidth  of  ^70  ps  is  short  com- 
pared to  the  lime  constants  involved  in  Xi(0. wc  *3n  simplify 
(II)  XX 

/,(0*=|r<'r,,)'  + a, Af(/) 

♦ /I.PAWA',"’"  - lf(»)|’ 

*Ao(1NoIK"  - C,I')*1(i)c-'kI,)'.  (13) 
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llcre  we  have  also  ignored  products  of  order  rj*  in  lt»c  second 
line,  since  r/1  « I. 

Fig.  7 »Im>wx  l.rj.  (13)  I'iullcil  fi*r  ly|Mc.i1  values  of  the  rele- 
vant lime  constants  as  we  slowly  Turn  on”  Ihe  excited-slate 
density-modulation  term  A£  and  (lie  ground-slate  density 
modulation  term  A p independently.  The  experimental  param- 
eters employed  arc 

jVP**//V o ■ 03 
r*  9.5  ns 
r#f  * 100  ns 
A K 5 |/m 

f,«2A5X  10*  cm/s 
/>»0.5  cm V*. 

It  clearly  seems  that  the  shape  of  our  experimental  results  can 
be  matched  only  by  assuming  (hat  density  modulation  of  the 
excited  state  is  the  dominant  mechanism  causing  the  acoustic 
effects.  Tliis  is  supported  by  the  data'shown  in  Fig.  8,  ob- 
tained from  the  third  plot  from  the  top  in  Fig.  3.  The  plot 
shows  both  the  exponential  decay  of  the  peak  amplitudes, 
which  gives  the  diffusion  coefficient  l),  and  also  the  decay  of 
the  dcplh-of-nmdulation  amplitudes  i.e.,  the  difference  be- 
tween the  fitting  exponential  and  the  observed  amplitude  at 
lire  minima.  The  decay  rates  arc  essentially  the  same,  which 
says  that  tire  oscillatory  term  decays  at  the  same  rate  as  the 
grating  component  of  /V,(.r),  not  A'0(x).  Only  if  density 
modulation  of  the  excited  state  is  the  dominant  mechanism 
can  this  depth  of  modulation  decay  be  explained.  For  ground 
state  modulation  the  decay  would  have  to  be  significantly 
slower. 

D.  Comparison  to  experimental  Results 

Fig.  9 repeats  the  experimental  result  of  Fig.  3,  second 
curve  from  the  top  (A  = 5 33  pm),  along  with  a semilog  plot 
of  the  peak  amplitudes  in.  this  curve,  showing  their  almost 
exact  exponential  decay.  At  the  peak  times,  the  instantaneous 
density  change  Ap(x,t)  is  zero,  and  lienee  the  observed  expo- 
nential decay  rate  for  Ihe  peaks  should  be  (2 It  + 2k*  O).  The 
observed  results  arc  in  fact  in  agreement  witlr  low-intensity, 
nonoscillating  measurements  which  give  r = 9.5  ns  and  D = 0.5 
cm1/5  for  this  particular  sample.  Also  shown  in  Fig.  9 is  a 
theoretical  fit  using  (13)  for  these  experimental  conditions, 
with  A0  ~ ‘0.02  and  A \ - + 0J7.  Evidently  it  is  density  mod- 
ulation of  the  excited  rather  than  the  ground  level  that  is 
dominant. 

From  all  our  other  measurements,  we  can  estimate  that 
o,/o0  *-«5,  and  from  the  heat  deposition  and  thcrmoacoustic 
properties  of  the  material,  w-c  can  estimate  that  in  the  high  in- 
tensity experiments,  the  normalized  maximum  density  change 
corresponds  to  i|~*  I0"3.  From  (12)  we  then  obtain  the 
rough  values 

(S04|)~-lOO 

(St  ♦ I) -4500 


at  the  peak  of  the  acoustic  wave  (where  A p/p0  <0).  There  is 
Independent  support  for  the  assumption  tlul  do0/«fy>  < 0,  if 
we  consider  the  ground  stale  absorption  spectrum  of  pcnla- 
cnie  in  the  vicinity  of  5.12  mil,  and  a Malic  pressure  tuning 
study  of  pent  acme  hues  [|0|.  No  data  on  the  pressure  tuning 
of  the  fust  singlet  excited  slate  of  pent  ace  nc  seems  to  he  avail- 
able, Iml  our  results  imply  that  the  net  cl  fee  ( at  5.12  urn  must 
be  of  opposite  sign,  and  several  times  larger,  llic  agreement 
shown  in  I ig  9 is  good,  especially  if  we  look  at  Ihe  simplicity 
of  the  theoretical  model,  In  which  we  did  not  include  highly 
nonlinear  effects  potentially  present  in  the  interaction  volume, 
such  as  the  following. 
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1)  In  tlic  lnf.1i  compressibility  regime  produced  by  I lie  luxe# 
induced  acomlic  waves  (excess  pleasure  * 50  aim).  Ilie  sound 
vctstcily  can  no  longer  be  coii\itlcicd  eoi»\la.»l  |ll|,aml  ibe 
acoustic  waves  can  no  longer  be  rigorously  treated  by  a linear* 
I zed  model. 

2)  A linear  dependence  of  cross  section  on  density  has  been 
assumed,  when  in  more  rigor  we  should  use  the  exact  nonlin* 
ear  absorption  lincsltapc  around  S32  nm. 

However,  we  believe  including  these  nonlinear  cITects  would 
only  change  Ihc.fit  in  Fig.  9 by  a negligible  amount,  and  in  the 
right  direction  for  an  even  more  exact  Tit. 

. IV.  Conclusion 

The  present  paper  reports  the  observation  of  amplitude  grat- 
ing effects  produced  by  coherent  acoustic  waves  induced  by 
optical  absorption  of  picosecond  pulses.  These  pulses  fr^Ti 
overlapping  transient  excited-state  and  thermal  gratings  in  the 
bulk  of  a molecular  crystal.  A simple  model  introduced  to 
treat  the  coherent  acoustooptic  interaction  fits  the  experi- 
mental results  with  quite  good  agreement.  Tire  results  yield 
information  about  a variety  of  processes,  including: 

1)  The  nature  and  rate  of  energy  migration  for  the  molec- 
ular excited  states,  obtained  by  studying  the  decay  of  the 
oscillatory  signal  peaks. 

2)  Information  regarding  high-frequency  acoustic  phonon 
processes  in  the  crystal,  including  acoustic  velocity  and  atten- 
uation. By  changing  the  grating  orientation  in  the  crystal, 
anisotropic  and  nonlinear  acoustic  effects  might  also  be 
directly  studied. 

3)  Information  regarding  the  pressure  or  density  modulation 
of  the  ground-state  and  cxcitcd-statc  absorption  in  the  crystal. 

The  amount  of  information  simultaneously  obtained  from 
measurements  such  as  these,  concerning  both  molecular  and 
crystalline  properties,  is  sizable.  We  hope  this  process  can  be 
naturally  extended  to  other  materials  as  well,  providing  a 
unified  understanding  of  liicir  properties. 
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Hf  I.  Transient  grating  concept  and  geometry.  Probe  pul>c  it  orth**c* 
onatiy  polarl/cd  relative  to  excitation  pubes,  and  signal  it  extracted 
with  rube  polarizer. 


Fig.  2.  Transient  grating  experimental  setup. 


Fig.  3.  High  excitation  Intensity  transient  grating  results  showing  the 
dependence  of  the  oscillation  period  T on  the  grating  frir  t spacing 
A,  for  four  experimental  results. 


n*.«-  Linear  dependence  of  A versus  T. 

i 


Fig.  5.  Density  wave  dynamical  behavior.  Frihge  peaks  are  at  x/A  * 0, 
1,2,*  • • . Excited  and  ground  state  population  spatial  distribution  at 

1-0. 


Fig.  6.  Graphical  (phasor)  description  of  acoustically  induced  phase  or 
amplitude  grating  added  to  excited  slate  amplitude  grating.  Resultant 
lime  dependences  are  radically  different.  Only  acoustically  induced 
amplitude  grating  effects  can  explain  observed  time  dependence. 


Fig.  7.  Excited  state  (upper)  modulation  and  ground  state  (lower) 
modulation,  (or  several  modulation  “depths”  A\  and  /t0.  Although 
behavior  may  look  similar,  the  oscillation  decay  is  radically  different, 
being  much  slower  for  ground  state  modulation. 


Fig.  8.  •-Peak  amplitudes  for  third  plot  in  Fig.  3 (A  * 4.7  pm),  show- 
ing excellent  exponential  decay.  Decay  time  constant  gives  D = 0.5 
cm*/*  for  this  particular  sample,  for  diffusion  atonr.Jfrc  b axis 
X- Depth  of  modulation  points  obtained  from  ininjpuiiii  oF  same 
plot.  Observed  decay  can  only  be  explained  if  an  acoustic  modula- 
tion of  excited  state  optical  properties  is  assumed.  Ground  state 
modulation  would  produce  a much  slower  decay. 


'Wl  / oni'TY) 


Fig.  9.  Theoretical  fit  to  second  experimental  plot  in  f ig.  3 (A  *»  5.3 
m pm),  for^o  * -0.02  and  /l|  * *0.3 7. 
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OPTICAL  PARAMETRIC  AMPLIFICATION 
Richard  A.  Baumgartner  and  Robert  L.  Byer 


I.  INTRODUCTION 

Experiments  in  nonlinear  optics  were  performed  in  1961 

soon  after  the  demonstration  of  the  laser. ^ These  experiments 

were  possible  due  to  the  increase  in  power  spectral  brightness 

made  possible  by  the  laser.  The  first  mixing  experiment  involving 

2 

three  optical  frequencies  was  performed  by  Wang  and  Racette  in 
1965.  Prior  to  that  experiment  the  possibility  of  parametric 

gain  in  a three  frequency  process  was  considered  theoretically 

3 4 5 

by  Kingston,  Kroll,  Akhmanov  and  Khoklov,  and  Armstrong, 

0 

Rloembergen,  Ducuing  and  Pershan  (ABDP). 

Parametric  oscillation,  an  important  extension  of  para- 
metric amplification,  was  first  achieved  in  1965  by  Giordimaine 
7 

and  Miller.  Parametric  oscillation  is  quite  useful  for  generation 

of  widely  tunable  coherent  radiation.  The  progress  in  parametric 

amplification  and  oscillation  has  been  the  subject  of  review  papers 

by  Harris, R Smith, ^ and  Byer.*^ 

Since  the  original  experiment  of  Wang  and  Racette,  laser 

sources  for  pumping  parametric  processes  have  improved  considerably. 

Pulsed  lasers  are  now  available  with  energies  greater  than  0.5  J/pulse 

at  repetition  rates  greater  than  ten  pulses  per  second  and  operating 

with  either  a TEM  transverse  mode  profile  or  various  plane  wave 

oo 

unstable  resonator  beam  outputs.  Recent  LiNbO^  parametric 
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I amplification  experiments  demonstrated  that  small  signal  pains 

of  50  and  saturated  conversion  efficiencies  of  20f  are  available. 
Consequently,  parametric  amplification  is  an  important  technique 
J for  amplification  and  power  generation  over  wide  frequency  ranges. 

L 

The  detailed  understanding  of  the  parametric  amplification  process 
is  also  useful  for  the  design  of  optical  parametric  oscillator 
tunable  sources. 

This  paper  examines  the  theoretical  solutions  of  parametric 

amplification  valid  in  the  depleted  pump  regime.  The  solutions 

are  compared  with  experimental  measurements  for  KD*P  and  LiNbO,^ 

optical  parametric  amplifiers  (OPA).  The  theoretical  derivations, 

0 

based  on  previous  work  by  ABDP,  Boyd  and  Kleinman  (BK)  and 
1 2 

Harris  and  computational  methods  are  considered  in  Sections  II 
and  III.  The  theoretical  treatment  includes  the  transition  from 
parametric  amplification  solutions  with  pump  depletion  to  solutions 
valid  where  pump  depletion  is  not  important.  A parametric  amplifier 
general  solution  program  with  and  without  time  dependence  is 
discussed  in  Section  III. 

The  experimental  measurements  are  presented  in  Section  IV. 

The  experiments  include  small  signal  gain  measurements  of  a 

* 

355  nm  pumped  KD  P OPA  and  its  angle  tuning  curve.  The  gain  and 
energy  conversion  efficiency  of  a 1064  nm  pumped  LiNbO^  OPA  is 
also  presented.  The  experimental  results  are  compared  with  the 
computer  solutions  presented  in  Section  III. 
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The  theoretical  and  experimental  results  and  comparisons 
presented  in  this  paper  are  the  first  detailed  study  of  the 
optical  parametric  amplifier.  Together  they  show  that  OPA ' s 
are  useful  devices  for  amplification  of  tunable  coherent  radiation 
in  the  optical  spectral  repion . 

I I . THFORY 

A.  Second  Order  Nonlinear  Interactions 

The  second  order  nonlinear  interaction  of  lipht  beams  is 
characterized  by  the  generation  of  a nonlinear  polarization 

P3^w3^  = eo  wl  ' w2  ^ : F1(“1)  e2^U)2^  C1) 

•a. 

where  P^(to^)  is  the  nonlinear  polarization  at  frequency  u)g  , 
e is  the  dielectric  constant,  >(-<!)_,  w,  , w0)  is  the  nonlinear 
susceptibility  tensor  and  and  are  the  interacting 

laser  fields. 

The  second  order  nonlinear  susceptibility  is  responsible 
for  sum  and  difference  frequency  generation  and  for  parametric 
amplification  in  the  interaction  of  two  laser  fields  in  a nonlinear 
medium.  Parametric  amplification  involves  three  frequencies 
related  by 

u>3  = + w2  • (2) 

The  field  at  is  assumed  to  bo  the  most  intense.  Energy  at 

either  o>,  or  w,,  incident  on  the  nonlinear  crystal  is  amplified 
with  the  remaining  non-incident  field  being  generated.  For  an 
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efficient  transfer  of  energy  to  occur  between  the  waves,  momentum 
or  phase  velocity  phase  matching,  defined  by  the  k vector  relation 


+ k0 


(3) 


must  be  accomplished,  where  | k^  | = 2ir  and  n^  is  the  index 

v of  refraction  in  the  crystal.  For  uniaxial  crystalline  media 

13 

with  birefringence,  phase  velocity  matching  can  be  achieved 
by  angle  tuning,  crystal  heating,  or  electro-optic  effect.  In 
the  present  experimental  work  Type  I angle  phase  matching  was 
used  with  the  signal  and  idler  (ujg)  fields  polarized  as 

ordinary  and  the  pump  field  polarized  as  extraordinary  waves 

* 

in  the  negative  birefringent  LiNbO^  and  KD  P crystals.  In  negative 
uniaxial  crystals  the  Type  I phase  matching  condition  for  collinear 
propagating  waves  is  given  by 


A3/ne^?3,e^  + 


where 


ne( >3,e) 


2 2 
sin  0 cos  0 


vv  vv 


(4) 


Here  9 is  the  propagation  angle  relative  to  the  optic  axis  and 
n (>)  and  n (1)  the  principle  extraordinary  and  ordinary 

G O 

crystalline  indices  of  refraction. 
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A phase  velocity  mismatch  factor 


Ak  = k3  ‘ kl  ~ k2 


is  a modification  of  Eq.(3)  when  phase  matching  is  not  quite 
achieved.  Possible  causes  of  non-zero  Ak  as  well  as  finite 
crystal  acceptance  angle  and  bandwidth,  are  discussed  in  Appendix  I 


The  pump  beam,  propagating  as  an  extraordinary  ray  in  the 
nonlinear  crystal  has  a power  flow  direction  at  an  angle  B to 
the  wavefront  normal.  The  Poynting  vector  walk-off  angle  for 
uniaxial  crystals  is  given  by 


tan  6 = 


2 2 

sin  0 cos  e ( n - n ) 

m m e o 

2 2 ^ 2 . 2 . 

n cos  f + n sin  f 
e m o m 


with  a positive  angle  B for  positive  birefringent  crystals. 
The  presence  of  B causes  a slight  adjustment  in  the  observed 
value  of  the  phasematching  angle  f given  by 


= e 


+ 6 


calc 


The  specific  nonlinear  polarization  relations  derived  from 


Eq.(l)  are 


P1  2 co  def f E2  E3 


P2  2 eo  def f E1  E3 


P3  2 eo  def f E1  E2 


(8a) 


(8b) 


(8c) 
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whore  the  notation  has  been  simplified  by  dropping  the  explicit 
frequency  dependence  and  the  susceptibility  has  been  written  in 


terms  of  an  effective  nonlinear  coefficient.  For  LiNbO^  with 
3 m point  group  symmetry 

deff  = d31  sin  9m  "d22  cos  9m  sin  3 ^ (9) 

. o * 

which  is  maximized  for  <}>  = -90  . For  KD  P with  42  m point  group 
symmetry  the  effective  nonlinear  coefficient  for  Type  I phase- 
matching is 

deff  = "d14  sin  9m  Sin  2 ^ (10) 


which  is  maximized  for  c|>  = -45°.  A derivation  of  the  effective 

nonlinear  coefficient  used  here  is  presented  in  Appendix  2 of  BK.  * 

The  coordinate  rotations  necessary  to  relate  d^j.  to  the  d^ 

coefficient  along  the  crystalographic  axes  are  discussed  in  Zornike 
14 

and  Midwinter. 


B.  Coupled.  Wave  Equations 

The  derivation  of  the  coupled  set  of  wave  equations  proceeds 
from  Maxwell's  equations 

'VxE  = - ||  FxH  = «T  + (11a) 

f-D  = 0 = 0 (lib) 
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_L 


and  tho  constitutive  relations 


D = eq  eF  + P 


C 12a) 


(12b) 


The  linear  polarization  is  contained  in  c so  that  P contains 
only  the  nonlinear  polarization  terms.  The  medium  is  assumed  to 
be  magnetically  inactive.  The  fields  obey  the  wave  equation 


2 9 

3 E . 3E  3 E 

3zX  ° 3t  ° ° 3t2 


= U 


32P 


°3t2 


(13) 


and  are  assumed  to  be  monochromatic  plane  waves  propagating  in 
the  near  field  in  the  z direction.  The  nonlinear  interaction 
occurs  over  distances  and  times  that  are  large  with  respect  to 
the  individual  sinusoidal  variations  of  optical  fields  so  that  an 
envelope  representation 


Vz.t)  "EmRc(Elnle«V  - Vb 


P.(z,t)  H„(Pim e3(6mt  “ V^) 

x ' ’ ' ^m  e lm 


(14a) 


(14b) 


is  used.  The  subscript  i represents  the  direction  of  polarization 
x or  y and  m the  frequency  of  interest  1,2  or  3.  The  Fqs.(14a 
and  14b)  are  substituted  into  Eq.(13)  and  the  slowly  varying 
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envelope  approximations  are  applied  leading  to 


9F 


m 


9z 


+ (n  /c) 
m 


3E  V] lco) 

— — + a E = -J_°— « p 
mm  m 


(15) 


9t 


2n 


m 


where  n = ( e ) " and  a = i V oc/n  is  the  electric  field 
loss  factor.  The  slowly  varying  envelope  approximation  is  quite 
valid  at  the  wavelengths,  intensities  and  pulse  lengths  relevant 
to  the  present  experimental  investigations.  The  envelope  represent- 
ation applied  to  the  nonlinear  polarization  components  given  by 
Eq.(8)  gives 


P1  - ’VWsV""1® 


(16a) 


P2  “ 2f0deffEIE3e'MkZ 


(16b) 


P3  = 2e0de(fElE2e+WkZ 


(16c) 


These  polarization  components  substituted  into  Eq.(15)  produce  the 
coupled  set  of  equations 


5E1  n1  dE.. 
cTz-  + c~  cHT  +0flEl 


(17a) 


+ 


n2  aE2 
c bt 


+«2E2 


-j 


^%£5e*e,o-3',1iz 

n2c  1 3 


(17b) 


dE„  n„  dE~ 

5TT  + c"  3t~  +p3E3  “ 


. E E e+-1dk7' 
n^ccos  /? 


(17c) 
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The  experiments  operate  in  frequency  regions  where  crystal 

2 

losses  are  negligible  i.e.  ~ <*3  = 0.  The  cos  B 

factor  in  Eq.(17c)  accounts  for  the  Poynting  vector  walkoff 

at  angle  B of  the  pump  field. 

An  improvement  in  the  coupled  equations  range  of  validity 

and  a mathematical  convenience  occurs  if  one  accounts  for  the 

slight  crystal  dispersion  contained  in  e . The  frequency 

dependence  of  e contained  in  the  definition  of  the  linear  part 

1 5 

of  the  polarization  is  written  as 

PLIN(t,z)  = e0  f*  x (w)E(to,z)e^a>tda)  (18) 

oo 

where  the  field  is  represented  as 

E(t,z)  = A(t,z)e+'j(wt  “ kz)  . (19) 

The  Fourier  transform  of  Eq.(19)  with  a Taylor  series  expansion 
of  k and  > , about  the  center  frequency  u>o  , when  substituted 
into  Eq.(18)  gives  to  the  first  order 

Plin(t,z)  = cQy  (w0)  A(t',z) 

where 

t'  = t - (n/c) [ 1 + ( o)/  [ 2e  ] ) z 

Previous  Ak  derivations  arising  from  e+l^a!t  factors 

remain  the  same  while  envelope  quantities  now  are  dependent  on 
t'  rather  than  t.  Thus  the  phase  velocity  terms  1/v  = n/c 

r 
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are  replaced  by  the  group  velocity  terms 


l/vp  = 3k/3u)  = ( n/c)  [ 1 + (w/[2e]) 


(20) 


in  Fq.(17).  Group  velocity  considerations  are  primarily  of 

interest  for  operation  of  parametric  amplifiers  with  short  pulses 

1 6 

in  the  region  of  the  tens  of  picoseconds . The  group  velocities  for 
the  coupled  waves  in  the  present  experiments  are  equal  to  within 
2 % so  that 


(21) 


is  a good  approximation  for  nanosecond  pulses. 

The  change  from  phase  velocity  to  group  velocity  in  the 

left  side  of  Fq.(17)  makes  a reduction  to  ordinary  differential 

1 7 

equations  possible  using  a transformation  indicated  by  Scott. 
The  transform  equations  are 

z -*■  r = v. 

t t = t - z/v 

K 

_3_  J 1_  3_ 

3 z 3 r v 3 x 


(22a) 

(22b) 

(22c) 


-> 


- 31 


3 

3t 


3 

3t 


( 22d) 


with  a resultant  simplification  of  the  coupled  equations  to 


. U|j-def  f * -jAkr 
J ^ c K2E3e 


J <E,e“^kr 

n2  C 


J -^d^f  f2  E.E2e+'1Akr 
n_c  cos  B 1 ^ 


The  use  of  e '^i  with  i = 1,2  and  3 accomplishes 

a rectangular  to  polar  dependent  variable  conversion.  The  polar 
form  equations  that  result  are 


dPl  - Wldeff  „ . „ 

5? P2p3sln0 


(23a) 


W2deff 

" ^p3Sin0 


(23b) 


W3  deff 


2“  Pi  PoSind 

n^ccos  /3  x A 


(23c) 


M = + “eft;  m3  P1P2  "2  U1  *2  A. 

5r  c Vcos2/?  A,  n2  p2  ~ SI  ^T)C' 


(23d) 


where 

9 = Ak  + <j)g  - <}>2  - 

An  invariant  for  the  parametric  conversion  process  representing  the 
power  flow  per  unit  area  parallel  to  the  direction  of  propagation  is 

W = (eQc/2)  (n.^2  + n2p 2 + ng  p2  cos2p)  (24a) 

or 

W = I1(0)  + I2(0)  + I3(0)  cos2  3 . (24b) 


Introducing  a set  of  normalized  dependent  and  independent  variables 
given  by 


U3 


y3n3 cos2 

4 k W / 


P3  = [^T3  C°s2  W)]' 


and 


f = 


4 deff  1T(:,r  r 


(eo  j ^ 2 ^ 3 n]n2n3  cos  6) 


(25a) 


(25b) 


(25c) 


( 25d) 


where 

AS  = Akr/ F (25e) 
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results  in  the  normalized  coupled  equations. 


‘1  _ 


- u2ug  sin  G 


(26a) 


2l  - - u,u„  sin  0 


(26b) 


. du„ 

— - = + u,u0  sin  0 

ac  1 2 


(26c) 


. — = AS  + 


/uiu; 

\ ua 


U2U3 


U1U3 


cos  0 


(26d) 


The  phase  equation  may  also  be  written  as 


d0  cos  0 d 

— = AS  + — [£n(u,u9u„)] 

d£  sin  0 d K 


The  phase  equation  can  be  integrated  with  a substitution  from 
Ea.  (26c.)  to  pive 


cos  0 = (T-  £ AS  u.^ ) /( u-j^Urj ) 


where 


= u^O)  u2(0)  u3(0)  cos  G + a AS  u^(0) 


is  a constant  of  integration 


r 


Energy  conservation,  power  conservation,  and  momentum 
conservation  for  the  three  wave  interactions  are  given  by 
the  three  auxiliary  equations 


and 


0)3  = u»l  + u>2 


1 = ^U^  + W2U  l + IDgUg 


K 

Ak  = - 
r 


dQ 

d? 


_ /ulu2  _ _ ulu3\ 

\ U3  U1  U2  / 


cos  9 


I ( 2) 

i c 28) 

1(26 d) 


Equation  (28)  is  a normalized  version  of  W.  The 
relative  phase  between  the  waves  is  related  to  momentum 
conservation  through  Ak  by  Eq.(26d).  Additional  invariants 
over  E are  derived  by  substitution  of  Eq.(2)  into  Fq. (28)  giving 


m., 


2 2 

u , + u„ 
1 3 


mr 


2 x 2 

U2  + U3 


m. 


2 2 
U1  U2 


(29a) 

(29b) 

(29c) 


which  arc 


the 


Manley-Rowe 


18 


relat ions . 
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I 


I 


The  solution  of  the  normalized  coupled  equations  is  based 
on  integrating  Eq.(26c)  after  substituting  the  integration  constant 
f and  the  Manley-Rowe  relations.  The  resultant  integral  is 


(30) 


The  cubic  polynomial  in  the  denominator  radical  suggests  the  use 
of  elliptic  integrals  as  solutions.  However,  several  algebric 
manipulations  are  required  to  produce  a standard  elliptic  integral 


form.  The  denominator  polynomial  has  roots  that  may  be  ordered  as 
2 2 2 

u^c  > u^b  > u3a  0.  Two  algebraic  substitutions 


(31a) 

(31b) 


convert  the  integral  expression  to  standard  elliptic.'  integrals. 
Bv  partitioning  £ into  ( f + £ ) - KQ  the  elliptic  integral 
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can  be  written  with  zero  lower  limits  as 


dy 

f(i  - y2)(i  - yV))" 


y(O) 


Jo 


dy 

fd  - y2T(i  - y2y2)  r 


The  inverse  operation  of  the  elliptic  integral  is  a Jacobian 
elliptic  function.  The  above  equation  re-written  with  the 
Jacobian  elliptic  functions  becomes 

y(S)  - sn  ^ (U3C  - U3a)J  ( c +co),Y) 

VC°)  - snr(u=o  - (?d)iY] 

The  peneral  solution  for  the  normalized  u values  is 


(32a) 

(32b) 


-Sc  r-> 

2 

= u3a 

+ <4,  - 

2 , 
“3a’ 

2r  , 2 

Sn  ^U3c 

o 

= [m2  - 

Ug(  K))  = 

-S(0) 

+ Ug(0  ) 

u*(  o 

=-[m1  - 

u2(  C)]  = 

U^fO) 

+ u2(0) 

u3a)''(£  + 50>'y) 

2,  t-X 

Ug(  O 


where  F 


r(sin_1[y(0)  ] ,y)/(u2c 


and  ?((/), y) 


standard  elliptic  integral . 


(33a) 

(33b) 


(33c) 


is  a 


C.  Special  Case  Solutions 

In  this  section  the  solutions  to  the  coupled  equations 
given  by  Fqs.  (33a,  33b  and  33c)  are  investigated  for  boundary 
conditions  appropriate  to  up-conversion,  sum  and  second  harmonic 
generation,  mixing  and  parametric  amplification.  Appropriate 
approximations  allow  the  elliptic  functions  to  be  reduced  to 
circular  and  hyperbolic  functions  in  the  non-pump  depleted  regime. 

The  frequency  up-conversion  process  is  described  by  the 
input  conditions 


causing 


u2(0)  = 0 (34a) 
u2(0)  < < u2(0)  (34b) 
u?„  = 0 and  r = 0 


where  is  the  sum  frequency,  is  the  frequency  of  the  highest 

intensity  wave,  and  is  the  lowest  frequency  which  is  up-converted. 
The  up-conversion  solutions  are 

u3(0  = u2(0)  sn2[u1(0)(f  + fo),y] 

u2(  0 = u2(0)  (l-sn2[u3(0)(  E + f^.y]) 

u2(f)  = u2(0)  -u2(0)  sn2fu1(0)(f  + K0),  Y] 

2 2 2 

Equations  (31)  and  (34)  yield  y = U2(0)/Uj(0)  <<  1 making 

2 

possible  the  approximation  sn(u,y)  ~ sin  u for  y <<  1 . 


Since  y(0)  = 0 from  Eqs.(32)  and  (34)  and  = 0,  inverting 

the  normalized  variables  given  by  Eq.(25)  yields  the  up- 
conversion  solutions 


I9(0)  ? 

I3(r)  = (u3/o)2)  sin  (r/A  ) 

cos  1 u 


(35a) 


and 


I9(r)  * I„(0)  cos^  r/ 1 ) 


u ' 


(35b) 


Ij(r)  - I1(0) 


(35c) 


where 


l/*u  = u -i  ( 0 ) f/r  = 


TT  deff[8I1(0)]1 

2 r 

Go  3 n1n2n3  cos  '&  c)2 


For  sum  generation  we  assume  the  input  conditions 


Uj(O)  = u2(0) 


(36a) 


u3(0)  = 0 


(36b) 


resulting  in  a denominator  polynomial  inside  the  radical  for  Eq.(30) 
of 

2 2 ^ 
t>3(m1  - u3) 

The  integration  simplifies  to 


r = L 


u3(r ) d u3 


-i. 


u3(0)  ml  " U3 


2 = []/u1(0)l  tanh  U;J(  0/^(0)  ] 
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t . • 


with  the  result  that 


M°)  2 

I3(r)  = (lOg/u,)  tanJT(r/£  ) 

cos  B 

I2(r)  = I2(0)  sech2  (r/£u) 

I]L(r ) = I1(0)  sech2  (r/tu) 

with 

I2(0)/I1(0)  = co2/oi1 


(37a) 


(37b) 


(37c) 


Sum  generation  reduces  to  second  harmonic  generation  if  the 
further  restrictions  w2  = 05 j » “3  = 2uj  and  1^(0)  = I2(0)  = 

\ I (0)  are  applied  yielding 


I 


I 


2w 


(r) 


0) 


(r) 


I (0)  2 

2~  tanh  (r/*SH) 

cos  B ••• 


Iu(0)  sech  (r/agH) 


(38a) 

(38b) 


l/£ 


sn 


’ deff 

nl^l  ^nl  C eo)2  cos  ® 


This  is  the  familiar  SIIG  solution  which  is  more  directly  derived 
through  an  energy  conservation  argument. 
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D.  Parametric  Amplification 

Parametric  amplification  of  energy  at  w by  a higher 
intensity  at  o>3  is  represented  by 


L 

U2(0)  = 0 

(39a) 

u2(0)  < u2(0) 

(39b) 

with 

U3a  = 0 * u3b  = m2  = u3(0)  > U3C  = m1  = u2(0)  + ui(°) 

and 


r = 0 

Since  y(0)  = 1 from  Eq.(31a),  the  elliptic  integral  becomes 
F[sin_1  y(0)  , y]  = F(ir/2,y)  = K(y) 

and  consequently 

= K(y)/[u2(0)  + u^(O)]1 

where  K(y)  is  a complete  elliptic  integral. 

The  identity  sn[u  + 2K(y),y]  = sn(u,y)  transforms  Eq.(33a) 
to 

1^(0  = u2(0)  sn2[u3(0)(?  - f0)/y,y]  (40a) 

y2  = u2(0)/[u2(0)  + u2(0) ] 
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(40b) 


For  convenience  since  6 is  a small  angle  we  set  cos  B = 1. 
The  parametric  amplification  solutions  which  follow  from  Eqs. 
(25),  (33)  and  (40)  are 


I2(r)  = I j (0)  + (aJ1/(o3)  I3(0) 

(l  - sn2  [(r  - ro)/«.,Y]j 

(41a) 

I2(r)  = (w2/w3)  I3(0)  (l  - 

sn2  [(r  - ro)/£,Y]) 

(41b) 

I3(r)  = !3( °)  an2  [(r 

- ro)/£,Y] 

(41c) 

where 


•n  cl  81  (0)  a . A 

1/!  ' ,,3(o)  r/(Yr)  = + Il<0)  -»/tI»<0> 

and 

rjl  = P3(o)  fQ/Y  = K(y)  = h *n  fl6[]  + I3(0)  w1/I1(0)  w3]J 

where  Fc].(25d)  and  the  above  expression  for  £o  with  an  expansion 
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of  K(  ) using  logarithm  have  been  used. 

For  small  values  of  f where  pump  depletion  is  not  important 

and 

o o o 

U1  • U2  < < U3  (42) 

further  equation  development  is  useful.  A combination  of 
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2 2 

Eqs . ( 25 ) and  (40)  results  in  (w/u)/  I/O)  = I/O)  Y /(I  - y ) 
which  modifies  Eq.(41)  to 

I/r)  = 1 1 ( 0 ) + [y2/(l  - Y2)]l1(0)  (l  - sn2[(r  - r(>)/l,y'])  (43a) 

I2(r)  = (w2/Wi)[y2/(1  - Y2)!  I/O)  (l  - sn2[(r  - r/^.y])  (43b) 

I3(r)  = I3(0)  sn2  [(r  - r//*^]  (43c) 


Identities  of  Jacobian  elliptic  functions  followed  by  the 
2 

approximation  y =1  permit  the  transformation 

[y2/(1  - Y2)]  (l  - sn2[(r  - ro)/i.,Y)j)  = y2  sn2(r/Z)/dn2(r/Z) 

- s.inh2  (r/£)  . 

The  result  is  the  familiar  relations  describing  amplification  in 
the  parametric  approximation  given  by 


- T1 

(0) 

cosh2( r / Z ) 

(44a) 

I/O) 

( ^ 

/ojg)  sinh2(r/fc) 

(44b) 

3(r) 

= 

I3(0) 

( 44c ) 

This  result  agrees  with  the  non-depletcd  pump  parametric 
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amplification  relations  previously  derived  by  Harris  and  Bycr . 

The  parametric  amplification  case  with  Ak  f 0 for  the  range 

Ala  < 2 (45) 

2 2 2 

and  Uj  » u2  < < u3  aIso  interest. 
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The  denominator  polynomial  of  Eq.(30)  for  Ak  f O becomes 

u3^m2  “ u3 ) " u2)  - (AS/2)2  (m2  - u2) 

2 2 2 

Except  for  the  change  from  u^a  =0  to  = (AS/2)  the 

development  proceeds  in  the  same  manner  as  the  Ak  = 0 solution. 
The  t / H in  Eq.(44)  becomes 


r/£  = (u2(0)  t f/r]2  - [Ak/2]2)2 


Introducing  the  parametric  gain  coefficient  rQ  defined  by 


Fo  = u3(0)  [ r/r] 


’Mff  8 V°> 

eo  >1  Vln2n3C 


(46) 


allows  the  solutions  to  be  written  in  a previously  derived  form 
as 

Ij(r)  = Ij (0)  cosh2  [T2  - (Ak/2)2]*  r (47a) 
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I2(r)  = (u2/«oi)I1(0)  sinh2  [ r2  - (Ak/2)2]*  r 


(47b) 


I3(r)  = I3(0)  (47c) 

III.  COMPUTER  METHODS  AND  EXAMPLE  SOLUTIONS 

Evaluation  of  Eq.(26)  was  accomplished  by  digital  calculation. 
A program  was  developed  to  calculate  1^  , 1^  , 1^  and  Q assuming 
plane  waves  with  user  supplied  initial  values.  To  initialize  the 


A 
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program,  values  for  crystal  length,  wavelengths,  refractive 

indicies,  nonlinear  coefficient,  and  walkoff  angle  were  entered. 

The  finite  difference  technique  used  in  the  computer  solution 

provided  a smooth  transition  between  regions  of  no  significant 

pump  depletion  and  the  high  conversion  efficiency  regime.  A 

computer  solution  based  on  the  evaluation  of  elliptic  integrals 

for  the  general  analytic  solution  had  difficulties  in  the  region 

20 

of  small  pump  depletion  as  noted  by  Bey  and  Tang. 

A computer  calculation  for  a parametric  amplifier  with  only 

a signal  wave  input,  in  addition  to  the  pump,  is  plotted  in  Figs. 

la  and  lb.  Figure  la  represents  the  relative  photon  fluxes  versus 

the  normalized  conversion  length  parameter  £ = r / H with  the 

associated  relative  phase  plotted  in  Fig.  lb.  These  plots  are 

2 

solutions  of  Eqs.(41).  The  P^  values  are  the  u^(f)  values 
2 

normalized  by  Ug(0). 

The  relative  signal  to  pump  wave  intensity  ratio  is  arbitrarly 
Chosen  in  a range  relevant  to  practical  parametric  amplifier 
...  taations.  At  the  optimal  conversion  point  , two  practical 

examples  of  the  pump  intensity  help  calibrate  £ . The  value 
of  dej. used  for  LiNbO^  is  5.58  • 10  ^ m/V  = 1550  nm  , 

>2  = 2400  nm)  and  for  KD*P  is  5.01  • 10  ^ m/V(^  = 580  nm  , 

?>2  = 913  nm)  Appendix  II  discusses  a Miller's  A derivation  of 
the  d..?  values  needed  to  calculate  d^^  in  Eqs.(9)  and  (10). 


- 45  - 


A computer  simulation  of  parametric  amplification  for  the 
interesting  case  where  both  P^O)  and  P2(°)  inputs  in  addition 
to  P_(0)  are  applied  to  the  parametric  medium  is  plotted  in 

O 

Fig.  2.  The  input  ratios  were  arbitrarly  chosen  to  be 
P1(0)/P3(0)  = 1.89-10-3  and  P2(0)/p3(0)  = 1.59-10-3  which  are 
near  typical  experimental  values.  Figure  2a  is  the  magnitude 
of  the  relative  photon  fluxes  while  Fig.  2b  plots  the  relative 
phase  from  an  arbitrary  initial  value  of  tt/4.  A study  of 
Eq.(26d)  for  various  input  conidtions  suggests  a first  order 
approximation  of  the  phase  as  a linear  function  of  £ from  the 
initial  value  9(0)  to  -u/2.  The  value  of  ^ , the  approximate 
point  where  9 first  becomes  -tt/2  , may  be  estimated  as 

r = -[it/ 2 + 9(0)] 

1 d£ 

After  the  exponential  conversion  region  the  pump  is  depleted  and 
the  phase  changes  to  + v/2  to  begin  the  sum  conversion  process 
which  re-generates  the  pump.  In  Fig.  lb  the  initial  phase  correction 
process  is  immediate  at  f = 0 since  Pg^l)  ~ 0 in  contrast  to 
Fig.  2b. 

A second  program  was  developed  to  describe  the  parametric 
interaction  for  Gaussian  time  envelope  inputs  for  the  three  waves. 

The  inputs  include  user  supplied  pulse  widths,  relative  time  offsets 
and  energy  fluence  values.  The  assumption  of  equal  group  velocities 
for  each  field  permits  the  inputs  to  be  divided  into  time  intervals 
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with  appropriate  average  intensity  values.  These  intensities 

become  inputs  for  the  first  computer  program.  The  sequence  of 

outputs  are  summed  to  provide  an  output  energy  per  unit  area,  E(  Q . 

A result  for  this  type  of  calculation  as  a function  of  £ is 

shown  in  Fig.  3.  The  peak  input  intensities  for  the  Gaussian 

__  o 

envelope  fields  are  in  a ratio  Wgl  (pk^ajjl^pk)  = 1.41  x 10  . The 

result  gives  conversion  efficiencies  that  are  more  realistic  for 

actual  pulsed  parametric  amplifier  operation.  The  Q i values 

2 2 

are  time  integrated  values  of  u.(f)/u„(0).  In  Section  IV  another 

JL  O 

computer  calculation  is  considered  which  also  takes  into  account 
the  transverse  spatial  variation  of  the  input  waves. 

IV.  PARAMETRIC  AMPLIFIER  EXPERIMENTS. 

A.  Introduction 

Parametric  amplifier  experiments  were  carried  out  using 

KD*P  and  LiNbO^  crystals  which  provide  parametric  amplification 

over  a tuning  range  that  extends  from  460  -1400  nm  and  1400-4000  nm. 

The  KD*P  crystal  was  pumped  with  the  third  harmonic  of  a Q-switched 

Nd:YAG  laser  while  the  I.iNbO,,  crystal  was  pumped  directly  at  1064  nm . 

o 

*2 

order  for  both  crystals.  The  crystals  were  available  in  large  sizes 
of  approximately  1 . 5 cm  diameter  by  5 cm  long  which  greatly 
facilitated  the  parametric  amplifier  measurements. 

The  parametric  gain  measurements  utilized  KD*P  crystals  due 
to  the  excellent  optical  quality  of  available  crystals.  The  LiNbO^ 
OPA  measurements  were  carried  out  in  conjunction  with  a remote 
air  pollution  monitoring  program  which  required  high  energy  tunable 


)g)  was  of  the  same 


The  parametric  gain,  proportional  to  ^ 
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The  LiNbOy  OPA  served 


21  22 

radiation  in  the  near  infrared  region.  ’ 
as  a final  conversion  stage  for  a 1064  nm  pumped  LiNbO^  parametric 
oscillator  source.  The  experiment  emphasized  the  overall  conversion 
efficiency  to  tunable  output. 

B.  KD*P  OPA  Experiments 
* 

The  KD  P OPA  measurements  were  conducted  to  verify  calculated 

OPA  gain  values  and  the  predicted  angle  tuning  curve.  Figure  4a 

shows  a schematic  of  the  experimental  configuration  used  to  measure 

parametric  gain.  Figure  4b  shows  the  experimental  set  up  which 

used  a dye  laser  input  to  verify  the  KD  P tuning  characteristics. 

In  both  cases,  the  laser  source  consisted  of  an  unstable  resonator 

Nd:YAG  osc il 1 ator /ampl if ier  system  which  generated  up  to  700  m J , 

23 

7 nsec  1064  pm  pulses  at  ten  pulses  per  second. 

jJ( 

For  the  KD  P gain  measurements  5f  of  the  1064  nm  beam 

was  selected  by  a beam  splitter  and  transmitted  through  two  polarizers 

for  variable  intensity  control.  The  beam  was  then  transmitted  through 

a two  to  one  beam  reducing  telescope  which  was  carefully  adjusted 

* 

to  provide  a well  collimated  idler  (m^)  beam  incident  on  the  KD  P 
OPA  crystal.  The  majority  of  the  1064  nm  energy  was  doubled  in 
a 2.5  cm  KD  P Type  II  angle  phasematched  crystal  producing  250  mJ 
pulses  at  532  nm . The  remaining  1064  nm  energy  was  summed  with 
the  532  nm  output  in  a 5.0  cm  KD  P Type  II  angle  phasematched 
crystal  to  produce  60  mJ , 6.5  r>s  pulses  of  355  nm  energy.  A prism 
disperser  separated  the  wavelengths.  The  horizontally  polarized 
355  nm  pumping  beam  was  reduced  in  diameter  and  collimated  by 
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1.7  to  1 telescope.  The  pump  and  input  idler  beams  were 

♦ 

combined,  passed  through  the  KD  P Type  I angle  phasematched 

parametric  amplifier  crystal,  and  separated  with  a prism  prior 

to  detection.  Small  uniform  intensity  regions  with  an  area  of 
—3  2 

7.85  x 10  cm  were  selected  using  apertures  to  allow  accurate 

24 

measurements  of  the  signal  and  idler  fluences.  The  corresponding 
input  fluences  for  the  pump  and  idler  waves  were  also  measured. 

j 

The  gains  obtained  from  the  measurement  of  plane  wave  small 
area  cases  are  compared  with  theory  in  Fig.  5.  The  dotted  line 
is  the  intensity  gain  calculated  using  the  non-depleted  pump 
approximation  given  by  Eq.(47).  The  solid  line  corresponds  to 
intensity  gains  calculated  with  the  input  idler  and  pump  fluences 
as  input  information  for  the  time  dependent  computer  calculations 
described  in  Section  III.  The  values  of  the  calculation,  including 
time  dependence  and  pump  depletion,  are  in  agreement  within  the 
experimental  error  of  the  measurements.  Pump  depletion  becomes 
more  significant  at  higher  intensities  consistent  with  the  increasing 
separation  betv.en  observed  results  and  the  non-depleted  pump  wave 
approximation . 

A slightly  modified  experimental  arrangement  shown  in  ^ig.db 

* 

was  used  to  perform  the  KD  P tuning  curve  verification.  The  532  nm 

energy  remaining  after  the  summing  crystal  was  used  to  side  pump 

25 

a dye  laser  oscillator.  ' The  pulsed  dye  laser  provided  an  input 
signal  wave  source  for  the  parametric  amplifier.  The  dye  laser 
generated  5 ns  pulses  of  0.2  to  0.5  mJ  energy  were  tunable  from 
550  nm  to  680  nm.  The  idler  wave  generated  in  the  parametric 
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amplifier  tuned  from  1000  nm  to  742  nm.  The  355  nm  pump  beam 
was  reduced  by  a 1.7  to  1 telescope  before  combining  with  the 
dye  laser.  The  spatial  overlap  of  the  dye  laser  Gaussian  mode 
and  the  unstable  resonator  pumping  beam  produced  modest  signal 
wave  gains  of  1.5, quite  adequate  for  tuning  curve  verification. 

Figure  6 shows  a comparison  of  the  calculated  and  measured 

phase  matching  points  for  KD*P.  The  phase  matching  curve  was 

2 6 

determined  from  a Sellmeier  expression  for  the  n and  ne 
indicics  of  refraction,  (see  Appendix  II),  combined  with  the  phase 
matching  condition  given  by  Eq.(4).  The  small  walk-off  angle 
correction  given  by  Eq.(7)  was  included. 

Computer  calculations  indicate  that  a substantial  amount  of 

* 

output  energy  should  be  available  from  a 5 cm  KD  P parametric 
amplifier  pumped  with  355  nm.  A 110  mJ  pumping  beam  with  an 
unstable  resonator  transverse  profile  and  a Gaussian  pulse  width 
of  6.5  ns  was  modeled  by  a set  of  three  concentric  rings  of 
graduated  intensities.  The  2 mJ  dye  input  pulses  at  580  nm  with 
a Gaussian  pulse  width  of  5 ns  used  the  same  ring  model.  The 
initial  calculations  ignored  pump  beam  walkoff.  A second  calculation 
re-examined  the  energy  contribution  to  the  output  from  the  most 
intense  dye  ring  taking  into  account  the  walkoff  displacement  of 
the  pumping  from  the  signal  and  idler.  Without  walkoff  calculations 
yielded  26.7  mJ  signal  wave  and  15.6  mJ  output  at  an  idler 
wavelength  of  913  nm.  This  is  an  overall  conversion  efficiency 
of  38 . 5£ . The  central  dye  ring  contributed  17.4  mJ  of  the  signal 
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and  10.3  mJ  of  the  idler  output  energy.  With  Poynting  vector 
vvalkoff  included,  the  central  dye  input  ring  contributed  15.1  mJ 
of  signal  and  8.9  mJ  of  idler  output  energy  for  an  efficiency  of 
20%. 

These  calculated  results  show  that  a carefully  designed 

* 

KD  P parametric  amplifier  can  convert  a significant  fraction  of 
the  input  pump  energy  to  tunable  signal  and  idler  output  energies. 

C.  LiNbOg  OPA  Experiments 

The  LiNbOg  OPA  experiments  were  designed  to  measure  the  OPA 

gain  and  energy  conversion  efficiency.  Figure  7a  and  7b  shows  a 

schematic  of  the  LiNbOg  OPA  experiments  with  a LiNbOg  0P0  and 

stimulated  Raman  scattering  input  sources.  The  angle  tuning  curve 

for  1.064  pm  pumped  LiNbOg  extends  from  1400-4000  nm  as  shown  in 

Fig.  8.  The  design  and  operational  characteristics  of  the  LiNbOg 

27 

0P0  are  considered  in  detail  in  a recent  paper.  For  the  present 

OPA  measurements  the  LiNbO,,  0P0  was  tuned  to  an  operating  wavelength 

o 

of  1900  nm  for  comparison  with  measurements  which  used  the  stimulated 
Raman  converter  fts  the  input  source. 

The  LiNbOg  parametric  oscillator  operated  at  56  mJ  input 
1064  nm  energy  at  two  times  above  threshold.  The  singly  resonate 
parametric  oscillator  utilized  a double  passed  pump  beam  to  reduce 
threshold  and  improve  the  output  stability.  To  avoid  feedback  into 
the  Nd:YAG  source  a Faraday  rotator  isolator  was  employed  prior  to 
the  1.5  - 1 beam  reducing  telescope  and  I..iNbO_  0P0.  The  generated 

O 
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OPO  output  was  transmitted  through  a polarizer  pair  for 

variable  attenuation  prior  to  the  LiNbO^  OP A. 

The  LiNbOg  OPA  was  pumped  by  350  mJ  8 nsec  1064  nm  pulses 

with  the  full  6.3  mm  diameter  near  field  unstable  resonator  beam. 

2 

The  smaller  area  0.13  cm  OPO  signal  beam  was  incident  directly 

into  the  OPA  to  avoid  extra  optical  components.  Thus  the  input 

2 

beam  did  not  utilize  the  full  0.30  cm  OPA  gain  region  in  this 
experiment.  A dichroic  beam  splitter  was  used  after  the  OPA  to 
divert  the  residual  1064  nm  pump  energy  prior  to  a power  meter 
which  measured  ihe  total  signal  and  idler  generated  output  energy. 
Figure  9a  shows  the  total  OPA  output  energy  available  as  a 

•• 

function  of  input  signal  wave  energy  for  pump  characteristics 

2 

intensities  of  52  and  69  MW/cm  . Figure  9 shows  OPA  gain  saturation 

by  plotting  the  total  signal  plus  idler  energy  output  divided  by 

signal  wave  energy  input  versus  the  signal  wave  energy  input.  The 
2 

52  MW/cm  pump  intensity  provides  a gain  of  55,  while  a gain  of 
136  is  realized  for  a 60  pm  input  in  both  cases.  The  LiNbO^  OPO 
followed  by  a LiNbO^  OPA  provides  the  capability  of  efficiently 
converting  input  pump  energy  to  tunable  output  energy  in  a manner 
that  is  analogous  to  dye  oscillator/amplifier  sources. 

A 1900  nm  signal  wave  source  generated  by  stimulated  Raman 
scattering  in  hydrogen  gas  was  also  used  an  an  input  to  the  OPA. 

The  experimental  arrangement  was  modified  as  shown  in  Fig.  7b  by 
substituting  a long  focal  length  telescope  and  the  one  meter, 
cell  for  the  Faraday  rotator,  beam  reducing  telescope  and  para- 
metric oscillator.  In  this  case  a 6 cm  long  LiNbO^  OPA  crystal  was 
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used.  Available  optics  permitted  a better  transverse  spatial 
overlap  between  the  input  Gaussian  signal  wave  with  spot  size 

= . 5 cm  and  the  unstable  resonator  pump  wave  than  in  the  OPO-OPA 
experiment.  The  one  meter  Raman  cell,  using  fused  silica  windows, 
operated  at  20  atmospheres  Hg  pressure.  The  threshold  for  the 
first  vibrational  Stokes  wave  was  18  mJ  for  a 1064  nm  beam  focused 
50  cm  into  the  cell. 

Figure  9b  shows  an  unsaturated  gain  value  of  140  for  a 
2 

73  MW/cm  pump  intensity.  This  gain  value  was  adjusted  to  an 

equivalent  gain  for  a 5.5  cm  long  OPA  crystal  used  in  the  0P0 

experiment.  The  change  in  transverse  beam  overlap  significantly 

improved  the  total  output  energy  available  at  higher  input  energies 

since  more  of  the  available  pumping  beam  area  is  now  utilized. 

Table  I shows  the  measured  parametric  amplifier  energy 

conversion  efficiencies  from  the  pump  to  signal  and  idler  waves 

for  various  input  pump  energies.  The  measurements  were  done  with 

2 

Aj  = 1900  nm  and  a pump  beam  area  of  0.45  cm  . The  conversion 

efficiency  increases  from  10  to  20%  for  pump  intensities  from 
2 

57.5  to  70  MW/cm  . The  effect  of  the  improved  transverse  spatial 
overlap  between  the  parametric  oscillator  driven  and  the  Raman 
coll  driven  parametric  amplifier  is  evident. 

The  possibility  of  extracting  additional  infrared  output 
energy  and  more  conversion  efficiency  from  the  pump  laser  beam 
was  confirmed  by  an  experiment  with  two  sequential  parametric 
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amplifier  crystals.  The  experimental  configuration  was  similar 

to  Fig. 7a,  but  with  the  polarizers  between  the  parametric 

oscillator  and  amplifier  source  removed.  A 4.1  cm  crystal  was 

located  after  the  pump  beam  and  input  beam  combining  beam  splitter 

followed  by  a second  5.5  cm  long  LiNbO^  crystal.  The  amplified 

signal  plus  idler  energy  was  measured  following  the  1064  nm  pump 

beam  diverter.  All  three  waves  traveled  directly  from  the  first 

crystal  to  the  second.  The  signal  wavelength  was  chosen  as 

1900  nm  with  the  idler  at  2420  nm . The  8 ns  pump  pulses  had  an 

2 

energy  of  248  mJ/pulsc  with  a pump  beam  area  of  0.45  cm  . The 
measurement  results  arc  shown  in  Table  II.  The  second  crystal 
does  provide  an  increase  in  conversion  efficiency  from  4 to  13% 
for  the  signal  wave  driven  case.  In  addition,  the  removal  of  the 
polarizer  permitted  the  excitation  of  the  first  LiNbO^  crystal 
by  both  the  signal  and  idler  waves  from  the  parametric  oscillator. 

A comp  ison  between  the  two  crystal  sequence  excited  by  the  signal 
wave  only  from  the  parametric  oscillator  and  signal  plus  idler 
shows  an  increase  from  32  mJ  output  to  43  mJ  with  a conversion  ■* 
efficiency  change  from  13%  to  17.4%. 

Figure  10  illustrates  the  input  and  output  spectra  for  the 
Raman  cell  driven  parametric  amplifier.  A 4.0  mJ  input  signal 
wave  at  1900  nm  was  used.  The  pump  energy  was  262  mJ  with  an 
output  signal  wave  energy  of  23  mJ.  The  output  linewidth  is 
predominately  controlled  by  the  input  bandwidth  of  0.4  cm  with 
some  broadening  due  to  OPA  gain  saturation. 
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V. 


SUMMARY  AND  CONCLUSIONS 


The  theoretical  derivations  of  Section  II  reviewed  and 
extended  the  coupled  wave  equation  analysis  for  parametric 

amplification.  The  connection  between  the  solutions  including  ( 

pump  depletion  and  earlier  derivations  assuming  non-depleted 
pump  waves  was  established.  A practical  method  of  generating 
the  Jacobian  elliptic  pump  depleted  solutions  with  a computer 

/ 

program  was  demonstrated.  As  an  example,  the  case  of  a parametric 

amplifier  with  Wj  , w2  and  w3  inputs  where  relative  phase  is 

important  in  the  conversion  process  was  considered.  A second  computer 

code  was  written  which  included  a Gaussian  pulse  time  dependence 

for  the  parametric  amplifier  inputs.  This  model  gave  more  realistic 

conversion  efficiency  values.  Small  area  gain  measurements  in  a 
* 

KD  P parametric  amplifier  vci-if ied  the  OPA  theoretical  calculations. 

The  KD*P  angle  tuning  curve  for  355  nm  pumping  was  also  carefully 

measured.  The  computer  calculations  wore  extended  to  include  beam 

transverse  profiles  which  included  a first  order  correction  for 

pump  beam  Poynting  vector  walkoff. 

Significant  values  of  gain  in  a LiNbO^  parametric  amplifier 

of  greater  than  50  for  both  signal  and  idler  waves  from  an  initial 

signal  wave  input  were  demonstrated.  The  importance  of  the 

analytical  derivations  were  underscored  by  thf>  saturation  behavior 

of  the  LiNbO,,  parametric  amplifier  which  generated  50  mJ  of  output 
o 

energy  for  a <1  mJ  input  energy  with  a conversion  efficiency  from 
the  pump  of  22 . 5r( . Finally,  additional  gain  and  conversion 
efficiency  was  demonstrated  by  a two  crystal  parametric  amplifier 


arra ngoment . 


* 

The  availability  of  high  optical  quality  KD  P and  LiNbO^ 
nonlinear  crystals  and  the  high  peak  power  near  diffraction 
limited  Nd:YAG  pump  source  has  led  to  a practical  demonstration 
of  optical  parametric  amplifier  capability.  The  OPA  advantages 
include:  a wide  tuning  range;  high  gain  in  the  forward  direction 
only,  thus  avoiding  the  need  for  isolation;  a relatively  narrow 
bandwidth  and  narrow  angle  gain  profile  which  avoids  super- 
fluorescence  problems;  and  high  conversion  efficiencies  in 
saturated  gain  operation  for  efficient  power  amplification.  These 
properties  make  parametric  amplifiers  useful  for  amplification  of 
dye  lasers,  F-center  lasers,  parametric  oscillators  and  tunable 
outputs  generated  by  mixing  or  summing  processes.  The  advantages 
of  parametric  amplifiers  which  are  well  known  in  the  microwave 
region  should  apply  equally  well  in  the  optical  region. 
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ENERGY  CONVERSION  EFFICIENCY  OF  LiNb0o  PARAMETRIC  AMPLIFIER 
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TABLE  II 


SEQUENTIAL  PARAMETRIC  AMPLIFIER  PERFORMANCE 


Crystals  at 
Phase- 
matching 
Angle 

Input  Fn 
from  0. 
Signal 

ergies 

P.0 

Idler 

Signal  and 
Idler  Output 
After  Second 
Crystal 

Energy  Conversion 
Efficiency 

1 Only 

2 . 5 m J 

0 

14  mJ 

4 % 

2 Only 

2.5  m J 

0 

14  mJ 

4 % 

1 , 2 

2 . 5 m J 

0 

32  mJ 

13  % 

1 , 2 

3 . 0 m J 

2.4  mJ 

43  mJ 

17.4  % 

Input  Pulse  Width  = 5 ns 

2 

Input  Pulse  Area  = 0.13  cm 


- 58  - 


* 


* 


APPENDIX  I 


CRYSTAL  ACCEPTANCE  ANGLE  AND  BANDWIDTH 


The  phase  velocity  mismatch  factor  is  defined  by 


(A.l) 


The  geometry  of  Fig.  A.l  illustrates  the  momentum  vectors  and 
associated  angles.  Assuming  small  angles  for  <p  and  if/  and 
taking  into  account  various  contributions  to  Alt  by  Taylor's 
series  expansions,  the  momentum  mismatch  expression  becomes 


Ak  = (u3/c) 


3n3(0) 

a© 


AG  + g 0 -bAu)^  , (A. 2) 


m 


where 


8n3(6)  ,,  2 .2  _ 2 2 Q 

= n n /[n  sm  9 ii  cos  9 ] 
e o'  o me  m 


a© 


g = k1k3/2k2 


and 


8k,  8kr 


b = 


aw,  a u)r 
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The  first  term  represents  an  beam  acceptance 

angle  in  the  horizontal  plane  with  the  input  beam  at 

angle  9^  from  the  crystal  optic  axis.  The  second  term  is 
either  the  beam  acceptance  angle  in  the  vertical  plane 

with  fixed  or  the  acceptance  angle  in  any  plane  with 


fixed.  The  final  term  relates  Ak  to  broad  linewidth 


inputs 


28,29 
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APPENDIX  II 


NONLINEAR  COEFFICIENT  VALUES  BY  MILLER’S  A SCALING 

The  value  of  d(-io0  , u>_.  , u>,_)  • ..  or  in  condensed  notation 

3 0 0 ljk 

30 

dA1  (see  BK)  is  estimated  by  utilizing  Miller's  A.  The 

nonlinear  tensor  element  is  related  to  Miller's  A and  linear 

. • • , 31 

susceptibilities  by 


“id* 


Eo  4ijk  Xkk"V 


(A2.1) 


where 


*ii  " (no  ~ 


at  a)., 


and 


" (no  - 15  at  "2 


Xkk  = <ne  - X)  at  “3 


Aijk  = L13  x 10  for  LiNb03  d31 


Aijk  = 5.99  x 10~3  for  LiNbOg  d22 


Aijk  = 4.04  • 10  2 for  KD*P  dJ4 


from  Choy  and  Byer 


30 


The  indicies  of  refraction  for  LiNbO„  are  available  from 

32  * 

previous  measurements,'  while  the  index  values  for  KD  P have  been 
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Y 


r)  r* 

curve  fit  to  the  following  Sellmeicr  equations  at  T = 25°C 
for  y in  vm  • 


1.23137 


0.2771624 


n^(KD*P)  = 1.0] 2233  + 


0 . 8381 R x 10 


-2 


1.00  - 


1.00  - 


15.00 

T2~ 


2 * 
n~(KD  P) 


1 .193722 

0.  933294  + — + 

0.7481954  x 10~2 

1.00 5 

r 


0.9423675 


1.00 


15.00 

~r 
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FIGURE  CAPTIONS 


Computer  solutions  of  Eqs.(26a-d)  for  P^(0)/Pg(0)  = 

1.89  x 10~3  , P2(0)  = 0 . 

a)  Plot  of  relative  photon  flux  magnitudes  versus  {;  where 
1 i = PiW(n3o)./niu)3) . 

b)  Plot  of  the  phase  angle  9 versus  £ 


Computer  solutions  of  Eqs.(26a-d)  for  P, (0)/P„(0)  = 

JL  O 

1.89  x 10~3  , P2(0)/P3(0)  = 1.59  x 10-3. 


a)  Plot  of  relative  photon  flux  magnitudes  versus  £ where 


1 i " ^(n^./n.^). 
b)  Plot  of  the  phase  angle 
initial  phase  of  +45°  . 


9 versus  £ starting  at  an 


Time  dependent  computer  solutions  for  Ti  = T2  ~ °-?5  t3, 

Q1(0)/Q3(0)  = 1.41  x 10-2  , Q2(0)  = 0 whore  the  Q values 

2 2 

are  the  time  integrated  values  of  u.(£)/u,,(0)  . The  output 
fluence  is  E.(  0 = Q . W( n„w . /n . ) . 

1 J-  1 1 «J 

* 

Experimental  configurations  for  the  KD  P OPA  investigations. 

* 

a)  Set  up  for  the  measurement  of  small  area  KD  P OPA  gain. 

* 

b)  Set  up  for  the  KD  P tuning  curve  verification. 

Signal  (w^)  and  idler  (w2)  intensity  gains  vs  pump  intensity 

* 

for  a KD  P OPA  pumped  at  354  nm.  The  estimate  is  the  para- 
metric approximation  given  by  Eq.(47).  The  calculation  is 
based  on  the  time  dependent  solutions  of  Eq.(26). 


- 63  - 


FIGURE  CAPTIONS 


coni  . 


afc 

6.  Angle  tuning  curve  for  a KD  P OPA  Type  I pha sematched 

with  355  nm  pump  at  T = 25°C.  The  solid  curve  represents 

26 

a computer  calculation  by  Herbst.  The  input  signal  wave- 
lengths are  shown  by  the  lines  and  the  observed  output  idler 
wavelengths  by  the  solid  dots. 

7.  Schematic  of  the  LiNbOg  OPA  experimental  configurations 

a)  LiNbOg  singly  resonate  OPO  source 

b)  A 1064  nm  pumped  H2  stimulated  Raman  source  at  = 1910  nm. 

8.  Angle  tuning  curve  for  the  LiNbOg  OPA  pumped  at  1064  nm 
at  T = 25°C. 

9.  a)  LiNbOg  OPA  energy  output  vs  input  signal  energy  for  the 
OPO  and  stimulated  Raman  1900  nm  sources. 

b)  LiNbOg  OPA  energy  ratio  vs  input  signal  energy  for  the 
LiNbOg  OPO  and  stimulated  Raman  source  inputs. 

10.  a)  4 mJ , 5 nsec  input  spectrum  at  1910  nm  to  the  LiNbOg  OPA 
generated  by  stimulated  Raman  scattering  in  pumped  by  1064  nm 
b)  LiNbO„  OPA  output  spectrum  at  23  mJ  output  energy. 

Al.  k-vector  diagram  showing  phasematching  angle  9m  and 

signal  (kj)  and  idler  (kg)  angles  with  respect  to  the  pump 
(kg)  vector. 
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RELATIVE  INTENSITY  RELATIVE  INTENSITY 
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INPUT  SPECTRUM 


Ht- 

SPECTROMETER 

RESOLUTION 


0.479  cm-* 


1909.8 


■Hh 


1910.0  1910.2  1910.4 

(Q)  WAVELENGTH  IN  nm 


1910.6 


OUTPUT  SPECTRUM 


SPECTROMETER 

RESOLUTION 


0.609  cm 


909.8  1910.0  1910.2  1910 

(b)  WAVELENGTH  IN  nm 


APPENDIX  III 

S.  J.  Brosnan  and  R.  L.  Byer,  " Optical  Parametric 
Oscillator  Threshold  and  Linewidth  Studies. 
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OPTICAL  PARAMETRIC  OSCILLATOR 
THRESHOLD  AND  LINEWIDTH  STUDIES 

S.J.  Brosnan  and  R.L.  Byer 
Applied  Physics  Department 


ABSTRACT 

This  paper  presents  a detailed  study  of  the  optimum 
design  parameters  for  the  LiNbO^  parametric  oscillator. 
Theoretical  and  experimental  studies  of  the  0P0  threshold 
parameters  and  of  linewidth  control  are  presented-  Consideration 
is  given  to  practical  factors  that  limit  0P0  performance  such  as 
laser  beam  quality  and  crystal  damage  mechanisms.  In  addition, 
stable  single  axial  mode  operation  is  reported. 
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OPTICAL  PARAMETRIC  OSCILLATOR 
THRESHOLD  AND  LINEWIDTH  STUDIES 

I.  INTRODUCTION 

The  optical  parametric  oscillator  (OPO)  has  been 

extensively  studied  and  developed  since  Giordmaine  and  Miller 

first  demonstrated  parametric  oscillation  in  LiNbO^  in  1965.* 

2 

Following  early  rapid  progress  reviewed  by  Harris  in  1969, 

work  in  parametric  oscillators  has  proceeded  at  a more  deliberate 

3 4 

pace  as  summarized  in  more  recent  reviews  by  Byer  and  by  Smith. 

q 

The  development  of  large,  good  quality  LiNbO^  crystals^ 
led  to  the  demonstration  of  the  Nd:YAG  laser  pumped  1.4  - 4.0  pm 
angle  tuned  LiNbO^  OPO  by  Herbst  et . al.,®  in  1974.  The  early 
results  were  soon  extended  to  higher  output  energies  and  narrower 

7 

linewidths  by  Byer  et . al . These  improvements  led  to  the 
application  of  the  LiNbO^  OPO  source  in  remote  air  pollution 

sensing, ^ laser  induced  chemistry, ***  and  laser  isotope  separation 

. , . 11 
studies . 

To  date  there  has  been  no  detailed  study  of  the  optimum 
design  parameters  for  the  LiNbO^  OPO  tunable  source.  This  paper 
fills  that  gap  by  presenting  theoretical  and  experimental  studies 
of  the  OPO  threshold  dependencies  and  linewidth  control.  Consider- 
ation is  also  given  to  practical  factors  that  limit  the  OPO 
performance  such  as  laser  beam  mode  quality  and  crystal  damage 
mechanisms  and  limits. 
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A model  for  describing  the  time  dependent  0P0  threshold 
pump  fluence  is  introduced  in  Section  II.  The  model  and 
computer  simulated  results  are  compared  with  detailed  experimental 
measurements  of  LiNbOg  0P0  threshold  as  a function  of  important 
parameters  such  as  pump  pulse  width,  cavity  length,  output 
coupling  and  pump  spatial  mode  parameters.  The  LiNbOg  0P0 
operating  limits  and  pump  laser  requirements  are  determined  by 
crystal  damage.  We,  therefore,  present  results  of  damage  studies 
in  Section  III.  In  Section  IV  a model  describing  the  LiNbOg  0P0 
linewidth  is  presented.  The  model  considers  the  effects  of  a 
grating,  beam  expander  and  etalons  within  the  0P0  cavity.  Using 
a grating  with  two  etalons  we  have  successfully  demonstrated 
stable  single  axial  mode  operation  of  the  LiNbOg  0P0  . 

The  theory  of  parametric  oscillation  requires  an  under- 
standing of  the  parametric  amplification  process  in  a nonlinear 

12 

medium.  In  an  accompanying  paper  Baumgartner  and  Byer  have 
treated  optical  parametric  amplification  (OPA)  in  detail.  We 
refer  to  the  OPA  results  in  our  theoretical  treatment  of  the 
optical  parametric  oscillator. 

II.  OPO  THRESHOLD 

A.  Theory 

In  this  section  a model  is  developed  describing  the 
threshold  pump  intensity  of  a pulsed  singly  resonant  oscillator 
(SRO)  including  the  effects  of  Poynting  vector  wallcoff.  Simplified 
expressions  for  OPO  threshold  are  derived  for  the  cases  of  a 
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simple  SRO  and  also  lor  the  double  pass  pumped  version  (DSRO). 

A schematic  of  the  simple  oscillator  is  shown  in  Pipe.  1.  A more 

complicated  resonator,  which  includes  line  narrowing  elements,  is 

easily  treated  in  this  model  by  generalizing  the  definition  of 

reflective  cavity  loss  R to  include  all  cavity  discrete  losses. 

The  tuning  curve  of  the  LiNbO^  0P0  is  shown  in  Fig.  2. 

Consider  a parametric  oscillator  with  a plane-parallel  two 

mirror  cavity.  Such  a configuration  in  a laser  is  unstable. 

Here,  however,  the  Gaussian  gain  profile  generated  by  the  Gaussian 

pump  quickly  confines  the  resonated  signal  wave  and  stabilizes 
13 

the  cavity.  ' A discussion  of  this  effect  and  a method  for 
calculating  signal  spot  size  are  given  in  Appendix  A. 

In  the  slowly  varying  envelope  approximation  the  equations 
describing  the  coll  inear  0P0  fields  in  terms  of  complex  phasor 
amplitudes,  assuming  no  pump  depletion,  are  given  by 


dE, 


dz 


+ a F = j k EE*  ejAkz 
s s ' s p 1 


dF.. 

l 

dz 


■n  . r r*  lAkz 
+ a . E . = .1  t ; . E e 
11  l p s 


(la) 


where  the  a’s  are  field  absorption  coefficients,  z is  the  spatial 

coordinate  in  the  propagation  direction  and  Ak  = k - k - k. 

* p S 1 

defines  wavevector  mismatch.  The  interaction  coefficients 
Krn  are  defined  by  the  following  equation. 
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“m  deff 


m 


(2) 


K 

m 


= s,i 


where  d is  the  effective  nonlinear  coefficient  and  n 

ef  f m 

is  the  refractive  index  of  the  nonlinear  crystal  at  angular 

frequency  wm  for  the  signal  or  idler  wave. 

In  Eq.(l)  the  signal  and  idler  fields  are  driven  bv  the 

polarization  waves  possessing  a spatial  profile  determined  by 

the  product  of  two  field  profiles.  We  assume  that  the  idler 

wave  is  not  resonated  and  is  free  to  accept  the  profile  of  its 

driving  polarization.  The  signal  wave,  however,  is  resonated. 

Its  driving  polarization  must  be  expanded  in  cavity  eigen  modes 

to  find  the  fraction  of  generated  signal  field  which  is  coupled 

into  the  TEM  mode  of  the  resonator.  Coupling  coefficients 
oo 

1 4 

of  this  type  have  been  calculated  by  Kogelnik  and  applied  to 

15  16 

OPO's  by  Boyd  et . al . and  Asby.  Here,  however,  it  is  convenient 
to  include  the  pump  Poynting  vector  walkoff  into  the  coupling 
coefficients  for  both  the  signal  and  idler  waves.  These  co- 
efficients and  a solution  to  the  coupled  equations  for  the  case 
ag  = = a and  Ak  = 0 are  given  in  Appendix  B. 

The  idler  field  is  assumed  to  be  zero  at  the  entrance  to 
the  crystal.  The  signal  field  at  the  end  of  a crystal  of 
length  1 is  given  by 

Eft)  = E (o)  e~atj  cosh  XX  (3) 

s s 
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where 


r ■ Avs'v2  <4) 

is  the  parametric  gain  coefficient. 

Here  is  the  signal  spatial  mode  coupling  coefficient 

defined  in  Appendix  B by 

2 

% 

Ss  + w2 

P s 

where  w and  w are  the  Gaussian  mode  electric  field  radii. 

P s 

An  effective  parametric  gain  length  is  introduced  where 


The  walkoff  length  £v,  is  given  by 


(5) 


SL 


w 


/fF 

2 


(6) 


where  p is  the  double  refraction  walkoff  angle.  The  walkoff 
length  is  closely  related  to  the  aperture  length  previously 
introduced  by  Boyd  et.al.15  For  large  arguments  the  error  function 
approaches  unity.  Therefore,  the  walkoff  length  is  the 

effective  gain  length  for  an  infinitely  long  crystal. 


f 
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The  single  pass  signal  power  gain  from  Eq.(3)  becomes 


|Esml2 

|ES(0)|2 


e 2aZ  cosh2  YCC 


(7) 


For  pulsed  0P0  operation,  the  buildup  to  threshold  may  be 
calculated  by  assuming  a Gaussian  time  profile  for  the  incident 
pump  intensity  which  yields  a time  dependent  gain  coefficient  V 
described  by 

r = r0  (8) 

2 

where  t is  the  1/e  intensity  halfwidth  of  the  pump  pulse. 

The  generated  signal  wave  is  amplified  as  it  makes  m cavity 

1 7 

transits  from  an  initial  parametric  noise  power.  During  a 

single  cavity  transit  the  pump  intensity  is  assumed  to  be  constant, 

or  equivalently,  t >>  c/n  £ . For  pass  m , the  signal  power  is 

P 

from  Eq . (7) , 

2 

P - P , [r  e~4aZ  cosh2[r  e“(tm/T)^]]  (9) 

m m-1  L o v ' 

The  factor  in  brackets  is  the  instantaneous  cavity  net  gain. 

Equation  (9)  can  be  easily  iterated  numerically  to  compute 

threshold,  incrementing  pump  intensity  until  a defined  threshold 

is  reached.  For  the  calculations  presented  here,  threshold  is 

defined  as  a signal  energy  of  100  y J » giving  a threshold  power. 

to  noise  power  ratio  of  £n(P  /P  ) = 33.  It  is  shown  in  the  next 

nr  o 
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section  that  this  numerical  model  is  in  excellent  agreement 
with  experiment. 

A simpler  model  is  desirable , however , for  optimizing  the 
0P0  parameters  that  affect  threshold.  To  allow  iterative 
calculations  we  introduce  a time  independent  gain  profile  of 
width  x and  magnitude  T . The  width  is  determined  by  the 
time  over  which  the  parametric  oscillator  has  net  instantaneous 
gain  greater  than  unity.  From  Eq.(9)  we  find  that 


(10) 


-1  /2 

where  the  relation  cosh  x = Jln(x  + /x  - 1)  has  been  used. 

Next,  we  relate  the  square  pulse  height  T to  the  Gaussian 
peak  T . For  consistency  the  total  gain  of  the  square  pulse  must 
equal  the  total  gain  of  the  Gaussian  case,  or 

Jin  cosh^  I '<£  - Jin  III  cosher  e ^q^T^  ]] 

2L  I q c ' 


= E Jin  cosh 


2 tr0*.-<V*>sl 


(11) 
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For  t >>  2L/c.,  this  summation  may  be  approximated  by  an 
integral . If  in  addition  the  hyperbolic  functions  are 
approximated  by  exponentials,  we  obtain 


Y9'  = r % 
o 


/? 

2 


erf ( t /2t  ) 
t /2x 


(12) 


Let  threshold  be  reached  after  the  p'th  round  trip,  where 
p = c7/2L  with  the  optical  cavity  length  L given  by 

L = L*  + ( n-1 ) H . (13) 

Here  h'  and  it  are  the  cavity  and  crystal  physical  lengths.  Re- 
writing Eq.(9)  in  the  form 

S,n(Pn/PQ)  = p Hn  R - 4 pa£  + 2 p in  cosh  F SC 

and  making  the  exponential  approximation  that  cosh  a = \ exp  2a 
we  find 

T£  = L/ct  Jin  Pn/PG  + 2a Jl  + Jln(l//Rj+  Jin  2 (14) 


This  equation  has  the  clear  physical  interpretation  of  a gain 
equals  loss  threshold  condition.  The  first  term  is  an  effective 
loss  due  to  the  buildup  time  during  pulsed  operation.  The 
following  terms  describe  distributed  crystal  absorption  losses 
and  the  cavity  output  coupling  loss.  The  final  term  Jln2  is  due 
to  the  SRO  operation. 


Substituting 

r y - 2//fF 

o 


Kq .(14)  into 

7 /2t 
erf(T  /2t ) 


Eq.(12)  we  have  finally 


L/ct  tn(Pn/Po)+  2a Jl  + Hn  (l/ Jin,  2 

(15) 
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Now  Eq.(10)  and  Fq.(15)  are  iteratively  solved  for  x 


and  T 

o 

The  threshold  pump  peak  intensity  using  Eqs.(2)  and  (4) 

2 

and  the  relation  L = I nc  e |E  I is 

p o 1 p 1 

"o  - <V2/™f2  <16> 


where 


K 


2 to  to . d 
s 1 


n n ^ n e 
s x p o 


2 

eff 


(17) 


For  LiNbO„ 
giving 


we  calculate  in  Appendix  C that 


■(9- 


k = (9.33  x 10 


-10 


2 \ , 
cm  \ 1 


MW  ' \\ 


deff=5.72  x 10~12  m/V, 
(18) 


with  Vs  in  cm. 

The  energy  fluence  (energy/area)  threshold  is  calculated 
by  integration  over  the  temporal  intensity  profile  and  gives, 


J *=  /|  T I = /| 

o V 2 o V 2 


(I’d2’) 


7 K 


g # 


(19) 


where  rQ  is  found  from  an  iterative  solution  of  Eqs.(10) 
and  (15)  and  ff  is  given  in  Eq.(5). 

For  threshold  fluence  estimates  one  can  simplify  the  model 
further  by  assuming  7 is  fixed  at  2 x,  which  is  approximately 
true  for  most  0P0  configurations.  With  this  assumption  from  Eq.(12) 
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and  Eq.(14)  we  find  that  t 2*  = 1.34  Y'/'  so  that  the  energy 
fluence  is  directly  given  by 


Jq(t  = 2x) 


2.25 

’ o T 


L 

2tc 


P 1 

Jin  — + 2a l + Jin  — + In  2 

P /R 

o 


2 


(20) 


These  are  the  SRO  threshold  fluence  values  for  a Gaussian  profile 
time  dependent  pump  pulse. 

A reduction  in  pump  fluence  threshold  and  some  enhancement 

1 8 

in  0P0  energy  efficiency  may  be  obtained  by  back  reflecting 
the  pump  radiation  to  double  pass  the  0P0  thus  creating  signal 
gain  on  both  the  forward  and  backward  transits  of  the  crystal. 

If  we  let  y be  the  ratio  of  backward  to  forward  pump  field 
amplitude  inside  the  crystal,  then  the  DSRO  relation  analogous 
to  Eq . (9)  is 

P = P R e~4a£  cosh2[r  SP  e“(tn/T)2]COsh2[r  SPy  e"(tm/T)2] 
m m-1  o o 


(21) 

If  the  square  gain  model  is  applied  to  the  DSRO,  the  result  is 


r«z? 


1 

1 + Y 


/ p \ 

1 

£n  [ — ) 

+ In  — + 2al  + £n  t 

\P  / 

/R 

' o' 

m 

(22) 


where  t is  determined  from 

R cosh2  hye-'P2^2]  cosh2' [ro*  y e-Cp2’>2].  l 


(23) 
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This  equation  is  satisfied  only  on  the  wings  of  the  gain 
envelope,  so  arguments  of  the  hyperbolic  cosines  are  small. 

In  this  approximation  it  can  be  shown  that  cosh  yx  = (cosh  x )^2 
We  then  obtain 

.2 


cosh  [ro*e-<T/-)2] 


1+Y 


1 

/IT 


(24) 


Thus  t is  given  by  Eq.(10)  with  R replaced  by  R'  , 


R = (R) 


1+Y 


(25) 


A calculation  of  the  DSRO  relation  between  T and  r shows 

o 

that  Eq.(12)  is  valid  for  a DSRO  also.  We  obtain  for  the  threshold 
gain , 


r .«/’  = — 

° Si 


r t / 2t 


orf(x  / 2t) 


J + Y/ 


— in 
CT  \P 


V 1 


+ 2a  i + Jin  — + in  4 

Sr 

(26) 


Now  Kqs.(10  and  (26)  are  solved  iteratively  to  find  T SC  . and 

v o 

threshold  fluence  JQ  is  found. from  Fq.(19). 

The  estimated  threshold  flucnce  for  7 = 2t  is  given  by 


2.25 


J = 
o 


(1  + Y)2 
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2ct 
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+ 2ai  + Jin  — + in  4 
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In  the  model  described  here,  an  implicit  assumption  has 
been  that  the  pump  beam  is  collimated  and  has  a uniform  phase- 
front  of  high  optical  quality.  For  an  SRO  the  requirement  for 
pump  collimation  is  somewhat  relaxed  since  the  idler  wave  is 
free  to  propagate  noncollinearly  with  phasematching  being 
preserved.  The  0P0  threshold  is  not  strongly  dependent  upon 
pump  collimation  if  the  divergence  cone  angle  is  roughly  less 
than  the  Poynting  vector  walkoff  angle. 

However,  threshold  can  be  markedly  increased  if  the  pump 
beam  optical  quality  is  nonuniform.  If  a transverse  phase 
variation  d^/dx  is  introduced,  then  as  the  pump  traverses  the 
crystal  with  Poynting  vector  walkoff  in  the  x direction,  the 
generated  signal  wave  sees  a pump  with  a changing  phase  \p  where 

Vz)"(Pfr)Z  (2S> 

From  Eq.(l)  we  see  that  the  transverse  phase  variation  acts  as 
an  effective  phase  mismatch  and  reduces  0P0  power  gain  by 
sine  <J*p(£).  The  critical  pha.se  variation  value  is  set  by 


djjj  _ jt_ 
dx  Pfc 


(29) 


In  a 5 cm  LiNbO^  0P0  crystal  the  critical  phase  variation  is 
approximately  a quarter  wave  phase  change  per  millimeter  of 
transverse  dimension.  High  quality  lenses  can  easily  meet  this 
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requirement,  but  commercial  grade  elements  may  not.  In  any 
1 case,  near  diffraction  limited  pump  field  and  imaging  optical 

components  are  recommended. 

I B.  Experimental  Results 

L 

In  this  section  the  SRO  theory  and  simplified  model  results 

are  compared  to  experimental  measurements  for  a L0.6  pm  pumped 

angle  tuned  LiNbO^  SRO.  The  pump  source  is  an  unstable  resonator 

19 

Nd:VAG  oscillator  operated  in  the  far  field  where  the  central 
lobe  of  the  modified  Airy  disk  function  closely  approximates  a 
Gaussian  beam.  The  SRO  threshold  energy  fluence  is  determined 
as  a function  of  cavity  length,  crystal  length,  signal  wave 
output  coupling,  pump  spot  size,  and  pump  pulse  width.  Agreement 
with  theory  is  very  good  for  the  highest  quality  LiNbO^  crystals 
and  input  pump  beam  profiles. 

In  making  these  measurements  a total  of  26  LiNbO^  crystals 

20 

were  tested  over  a nine  month  period.  The  threshold  energy 
fluence  of  the  best  crystal  tested  agreed  very  well  with  theory 
and  was  used  to  scale  the  experimental  data  to  test  the  theory. 

The  threshold  energy  fluence  for  an  average  crystal  was  1.7  times 
, greater  than  that  of  the  best  crystal.  The  standard  deviation 

in  threshold  fluence  for  the  crystals  tested  was  0.45.  Thus  the 
best  crystal  was  exceptional  and  as  yet  not  generally  reproducible. 
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The  Nd:YAG  laser  pump  beam  phasefront  quality  was 
another  critical  factor  for  an  accurate  experimental  to 
theoretical  comparison  of  threshold.  Beam  distortion  in  the 
output  of  some  Nd:YAG  lasers  is  severe  enough  to  prevent  0P0 
oscillation.  In  such  cases  propagation  to  the  far  field  to 
allow  diffraction  to  clean  up  phase  perturbations  may  help. 

Beam  reduction  or  expansion  telescopes  using  commercial  grade 
lenses  induced  enough  wavefront  distortion  to  increase  the  0P0 
threshold  by  60%.  It  should  be  noted  that  the  near  field  phase- 
front  of  the  Nd:YAG  laser  source  used  in  these  experiments  was 
of  good  quality  as  determined  by  interferometry  and  by  the 
fractional  conversion  to  the  central  lobe  of  the  far  field  Airy 
disk  profile. 

The  simple  SRO  cavity  configurations  are  shown  in  Fig.  1. 
Figure  2 shows  the  LiNbOg  angle  tuning  curve  and  the  crystal 
gain  bandwidth  determined  by  crystal  dispersion.  For  the 
theoretical-experimental  comparisons  discussed  in  this  section, 
the  simple  SRO  configuration  was  used. 

Figure  3 shows  the  measured  SRO  threshold  fluence  vs 
cavity  length  for  a 4 mm  diameter,  20  nsec  duration  pump  beam. 

The  solid  line  is  derived  from  the  numerical  solution  of  the 
SRO  threshold  condition  given  by  Fq.(9)  and  is  in  excellent 
agreement  with  experiment.  The  nonlinear  coefficient  used  for 
the  theoretical  calculation  is  discussed  in  Appendix  C and  yields 
the  gain  constant  given  by  Eq.(17).  The  dashed  line  showing  the 
simplified  square  pump  pulse  model  results  given  by  Eq.(19)  is 
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within  of  the  numerical  results.  Figure  3 shows  that  the 

simplified  threshold  model  is  of  adequate  accuracy  for  design 
purposes.  Its  advantage  lies  in  the  ease  with  which  it  can  be 
applied  to  0P0  design  calculations  where  a number  of  variables 
affecting  threshold  are  considered. 

The  FRO  threshold  fluence  vs  LiNbOg  crystal  length  is 
shown  in  Fig.  4.  For  these  measurements  a number  of  different 
crystals  of  comparable  quality  were  used  to  minimize  scatter  in 
the  measurements.  The  threshold  values  were  then  scaled  to  the 
value  in  agreement  with  the  best  5 cm  long  crystal  tested.  The 
effect  of  pump  beam  Poynting  vector  walkoff  is  evident  for  the 
longest  crystals  tested.  From  the  model  and  these  measurements 
the  optimum  crystal  length  is  approximately  5 cm. 

The  SRO  threshold  fluence  dependence  on  signal  reflectance 
is  shown  in  Fig.  f>.  The  data  was  taken  by  using  a series  of 
mirrors  coated  for  a range  of  signal  reflectances.  However,  each 
mirror  also  had  a different  reflectance  at  the  pump  and  idler 
wavelengths.  The  data  points  have  been  individually  scaled 
according  to  the  DSRO  to  SRO  ratio  to  reduce  the  variation  on 
threshold  due  to  pump  feedback..  The  present  model  does  not 
include  the  effect  on  SRO  threshold  of  varying  amounts  of  feed- 
back at  the  idler  field.  The  problem  of  parametric  amplification 
with  both  signal  and  idler  present  was  discussed  in  the  previous 
paper.  It  was  shown  that  if  the  idler  is  present,  but  of  the 
incorrect  phase,  the  parametric  gain  in  the  crystal  is  effectively 
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delayed  until  the  phasing  is  corrected.  The  re-phasing  occurs 
by  idler  depletion  to  noise  followed  by  re-growth  at  the  proper 
phase  and  thus  reduces  the  effective  crystal  gain  length.  The 
error  brackets  in  Fig.  5 are  an  estimate  of  this  effect. 

The  pump  spot  size  dependence  of  parametric  gain  for 

15 

critically  phasematched  interactions  is  well  known.  Figure  6 
shows  the  first  experimental  verification  of  the  decrease  in  gain 
for  small  pump  spot  sizes  as  a result  of  Poynting  vector  walkoff. 

The  model  and  experimental  results  are  in  good  agreement  and 
show  that  for  a 5 cm  long  LiNbO^  SRO  pump,  spot  sizes  of  A rmn 
diameter  or  greater  are  preferred  to  minimize  the  threshold  fluence. 

The  model  can  also  be  applied  to  the  0P0  pumped  by  the 
near  field  unstable  resonator  beam  profile.  The  measured  and 
calculated  threshold  fluence  for  the  near  field  pump  beam  is  50% 
higher  than  that  for  the  far  field.  Therefore  to  obtain  the 
lowest  operating  0P0  threshold  far  field  conversion  of  the  unstable 
resonator  mode  is  preferred.  In  the  present  experiments  the  far 
field  conversion  was  accomplished  by  propagating  the  beam  for  a 
distance  of  20  m by  multiple  reflections  between  plane  mirrors. 

The  pump  pulsowidths  dependence  of  threshold  fluence  is  shown 
in  Fig.  7.  The  pump  pulsewidth  was  adjusted  by  operating  the 
Nd:YAG  oscillator  at  varying  degrees  above  threshold  and  using 
a following  Nd:YAG  amplifier  to  provide:  the  required  energy.  'Die 
experimental  values  shown  in  Fig.  7 are  somewhat  higher  than  theory 
in  this  case  because  of  the  pump  phase  front  distortion  caused  by 
the  Nd:YAG  amplifier  crystal  and  also  because  the  measurements 
were  made  in  the  near  field.  It  is  clear  from  the  figure  that  the 
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data,  when  scaled  by  .37,  agrees  well  with  theory. 

The  results  shown  in  Fig.  7 are  for  a short  0P0  cavity. 

For  longer  cavity  lengths  near  15-20  cm  required  to  accommodate 
the  linewidth  control  elements,  there  exists  a broad  threshold 
fluence  minimum  resulting  in  nearly  constant  threshold  values 
for  pulsewidths  in  the  10-30  ysec  range.  For  shorter  pulse 
widths  the  threshold  increases  due  to  buildup  time  effects. 

For  longer  pulse  lengths  the  threshold  approaches  steady  state 
values  and  is  constant  in  pump  intensity. 

A significant  improvement  in  SRO  performance  is  possible 

by  double  passing  the  pump  beam  (see  Fig.  1).  In  addition  to  a 

1 8 

threshold  reduction  first  noted  by  Bjorkholm  and  described  by 
Fq.(26),  this  arrangement  improves  the  oscillator  stability,  increase 
the  conversion  efficiency  and  provides  for  separation  of  the  tunable 
0P0  output  from  the  incident  pump  beam.  With  a pump  high  reflector 
mirror  after  the  0P0  output  mirror,  typically  GO % of  the  forward 
propagating  pump  beam  is  reflected  back  into  the  crystal.  The  rest 
is  lost  to  mirror  and  crystal  surface  reflections.  For  the  present 
case  the  theoretically  predicted  DSRO  to  SKO  threshold  reduction  is 
0.5'i  which  compares  to  a measured  threshold  reduction  of  0.68  i .15. 

The  DSRO  configuration  can  be  simplified  by  coating  the 
output  mirror  to  reflect  both  the  resonated  signal  wave  and  the 
pump  wave.  Coatings  of  this  type  have  been  designed  and  used  in 
DSRO  operation  in  spite  of  the  more  complex  coating  requirements. 
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The  advantages  of  the  DSRO  conf iguration  are  not  as 
clear  cut  as  it  appears,  however.  The  net  threshold  fluence 
reduction  is  real,  but  the  fluence  at  the  crystal  surface, 
which  governs  crystal  damage  is  the  sum  of  the  forward  and 
backward  travelling  pump  waves  and  in  fact  may  be  larger  for  the 
DSRO  at  threshold  than  for  the  SRO.  The  double  pass  pumping  of 
the  LiNb03  crystal  also  leads  to  a reduction  in  the  intensity 
at  which  small  scale  self  focusing  occurs, thus  increasing  the 
risk  of  this  type  of  damage.  Finally,  the  reflected  pump  beam 
must  be  isolated  from  the  pump  laser  source  to  prevent  the 
possibility  of  damage  to  the  laser  oscillator  from  excessive 
feedback.  A pulsed  Faraday  rotator  isolator  is  usually  used 

2i 

for  near  field  operation  of  the  SRO  and  DSRO.  For  operation 
in  the  far  field,  the  20  m one  way  propagation  distance  to  the 
0P0  provides  a round  trip  time  delay  of  120  nsec  which  is 
sufficient  temporal  isolation.  For  SRO  operation  in  the  far 
field  no  feedback  is  observed.  For  DSRO  operation  some  feedback 
influence  is  evident,  but  is  not  strong  enough  to  affect  the  pump 
laser  oscillato'  or  the  DSRO  performance. 

The  experimental  or  useful  conversion  efficiency  n , 
defined  as  the  ratio  of  the  external  signal  plus  idler  output 
energy  to  incident  pump  energy,  is  shown  in  Fig.  8 for  the  SRO. 

At  three  times  threshold  the  observed  energy  conversion  efficiency 
is  18%.  Bjorkholm  has  considered  the  power  conversion  efficiency 
for  the  SRO  pumped  by  a plane  wave  and  Gaussian  intensity  beam. 
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‘ In  his  analysis  an  integral  over  the  spatial  mode  profile  of 

| the  pump  wave  was  performed  yielding  a fractional  transmitted 

pump  power  vs  the  number  of  times  above  tb  -eshold  N.  For 

pulsed  0i>0  operation  considered  here,  a temporal  integral  must 

also  be  performed  to  obtain  the  energy  conversion  efficiency  from 

the  peak  power  efficiency.  The  result  is  that  the  peak  power 

conversion  efficiency  must  be  scaled  by  the  ratio  of  the  signal- 

to-pump  pulsewidth,  which  is  approximately  given  by  1-1/N. 

Finally,  to  compare  Bjorkholm's  results  to  the  measured  conversion 

efficiency,  the  internal  conversion  efficiency  must  be  scaled  by 

the  crystal  surface  and  mirror  transmittance  losses.  For  the 

simple  SRO  cavity  considered  here  the  scale  factor  is  0.6.  The 

agreement  between  the  theoretically  expected  conversion  efficiency 

and  the  measured  efficiency  is  quite  good  as  shown  in  Fig.  8.  The 

conversion  efficiency  of  the  SRO  continues  to  increase  monotonical 1 y 

7 

with  N as  shown  earlier  by  Ryer  et.al.  However,  reliable 

operation  at  values  of  N much  greater  than  three  is  difficult 

due  to  the  limitations  imposed  by  crystal  damage.  A more  appropriate 

means  of  increasing  the  net  0P0  output  energy  is  to  increase  the 

pump  beam  area  but  maintain  operation  at  three  times  threshold  for 

the  maximum  available  laser  energy.  An  alternative  is  to  utilize 

an  optical  parametric  amplifier  (OPA)  after  the  0P0  as  an  efficient 

energy  converter.  The  characteristics  of  the  OPA  as  a power  amplifier 

1 2 

were  discussed  in  the  previous  paper. 


- 101 


III.  DAMAGE  LIMITATIONS 

A good  understanding  of  the  factors  affecting  the  0P0 

threshold  coupled  with  a detailed  knowledge  of  LiNbOg  damage 

mechanisms  and  limits  is  required  for  the  design  of  a reliable 

oscillator.  In  this  section  experimental  and  theoretical  results 

are  presented  which  describe  the  observed  LiNbOg  damage  mechanisms. 

Surface  damage  is  discussed  first  followed  by  a discussion  of  small- 

scale  self -focusing  which  leads  to  internal  tracking  or  "angel 

hair"  damage  of  LiNbOg  crystals. 

The  surface  damage  of  LiNbOg  has  been  extensively  considered 

22  23  2d 

by  a number  of  workers.  In  a series  of  papers  Zverev  et.a}.,  ’ ’ 

proposed  that  LiNbOg  surface  damage  is  caused  by  thermal  heating 
of  an  oxygen  depleted  niobium  surface  layer.  In  addition,  the 
reduced  niobium  is  assumed  to  be  concentrated  at  crystal  imperfect- 
ions, inclusions,  or  microcracks.  Damage  measurements,  initially 
made  by  Zverev  et.al.,  and  recently  repeated  in  this  laboratory, 
are  in  agreement  with  this  basic  model . Measurements  of  surface 
energy  damage  as  a function  of  spot  size  show  that  the  damage 

threshold  fluence  rises  sharply  for  pump  spot  sizes  less  than 

2 

60-70  pm  . The  damage  fluence  varies  from  2.7  J/cm  for  large  spot 

2 25 

sizes  to  greater  than  11  J/cm  for  tightly  focused  spots.  This 

indicates  a mean  separation  of  surface  defects  on  this  order  and 

is  not  inconsistent  with  the  recent  observation  of  microdomains 

26  27 

and  antidomains  in  LiNbOg  and  LiTaOg,  ’ J at  this  density. 
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The  pulsewidth  dependence  of  the  damage  fluenoe  can  be 

24 

found  by  modeling  the  trap  sites  as  an  absorbing  layer. 

28 

Using  thermal  diffusion  arguments  as  outlined  by  Bliss  we 
find  that  for  short  pulses  during  which  the  layer  is  not  completely 
heated  the  fluence  limit  scales  as  /t  . When  complete  heating 
occurs,  but  pulsewidths  are  short  compared  to  the  thermal 
diffusion  time  into  the  bulk,  the  damage  fluence  is  independent 
of  x . For  longer  pulses,  diffusion  occurs  during  the  pulse 
and  the  damage  fluence  scales  as  t , thereby  becoming  intensity 
limited . 

24 

Zverev  reports  a constant  fluence  damage  threshold  of 
2 

4.6  J/cm  for  pulsewidths  less  than  30  nsec  and  a constant 

2 

intensity  limit  of  170  MW/cm  for  longer  pulses.  In  an  earlier 

2 

paper  he  quotes  a damage  limit  of  120  MW/cm  for  a 30  nsec  pulse. 

The  result's  of  damage  measurements  made  in  our  laboratory 

at  1.06  pm  are  shown  in  Fig.  9a.  We  measure  a constant  fluence 

2 

damage  limit  of  2.7  ± .4  J/cm  for  uncoated  LiNbO^  in  the  10-30  nsec 
region.  This  value  is  somewhat  lower  than  Zverev's  results  . For 
reference,  Fig. 9 also  shows  the  energy  fluence  just  outside  the 
surface  of  an  average  LiNbO^  crystal  operating  at  three  times 
threshold  for  a short  8 cm  cavity  SRO.  Clearly  some  means  of 
raising  the  surface  damage  limit  is  desirable  for  long  term 
damage  free  OFO  operation. 

The  niobium  ion  reduction  model,  proposed  by  Zverev, 

provided  the  incentive  to  test  the  effect  of  oxygen  on  the  damage 

29 

fluence  level  of  LiNbO^.  The  measurements  were  performed  by 
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mounting  LiNbOg  test  plates  in  an  oxygen  atmosphere  in  a 
temperature  controlled  environment.  Surface  damage  was  defined 
as  the  maximum  fluence  at  which  the  crystal  surface  did  not 
damage  when  irradicated  for  one  minute  at  10  pps  by  Q-switched 
laser  pulses.  All  measurements  were  carefully  normalized  to  a 
standard  room  temperature  uncoated  LiNbOg  test  plate.  It  has 
been  shown  that  the  damage  fluence  is  strongly  dependent  upon 
the  number  of  integrated  pulses  as  well  as  pulse  spatial  diameter 
and  spatial  mode  quality.  The  1 minute  survival  time  was  chosen 
to  eliminate  the  dependence  on  integrated  pulses  and  to  be 
representative  of  actual  0P0  operating  damage  levels.  The  spot 
size  was  set  significantly  larger  than  the  60  - 70  pm  defect 
spacing . 

When  oxygen  was  supplied  to  the  crystal,  even  at  room 
temperature,  a significant  increase  in  damage  level  resulted  as 
shown  in  Fig.  9a.  This  result  was  not  totally  unexpected  since 
it  has  been  known  for  some  time  that  oxygen  diffuses  rapidly  into 
and  out  of  LiNbOg  crystals  and  that  SiOg  anti-ref lection  coatings 
significantly  increase  the  LiNbOg  surface  damage  level  as  shown 
in  Fig.  9a.  It  is  not  clear  whether  the  RiOg  coating  itself 
or  the  coating  procedure  acts  to  replace  the  LiNbOg  surface  oxygen 
and  yield  the  increased  damage  threshold.  It  is  known  that  LiNbOg 
anti-reflection  coated  with  ThF^  shows  the  same  surface  damage 
threshold  as  the  uncoated  crystal.  The  convenience  of  the  SiOg 
coating  and  its  anti-reflection  reduction  of  surface  Fresnel 
reflection  losses  to  1-2%  clearly  make  it  the  choice  for  LiNbOg 
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OPO  operation.  However,  it  should  be  kept  in  mind  that  the 
Si0o  coating  must  be  flawless  since  imperfections  lead  to  phase 
perturbations  of  the  pump  beam  that  can  aggravate  the  onset  of 
self  focusing. 

Small-scale  self-focusing  using  an  instability  theory 

11  12 

approach  is  given  by  Suydam.*  ’ Small  spatial  beam  perturbations 
in  amplitude  and  phase  are  amplified  through  the  nonlinear  intensity 
dependent  index  of  refraction  ng  resulting  in  beam  breakup  and 
focusing . 

For  a Gaussian  beam  the  critical  power  Pc  for  whole  beam 
self -focusing  is  given  by 

,2 

AG  ri 

P„  = % x 10"'  (30) 

c 32  tt  rig 


where  Pc  is  in  Watts  and  the  other  parameters  are  in  cgs  units. 
For  small-scale  self-focusing  the  intensity  damage  limit  given 
by  Suydam  is 


I 


d 


X cn 
o 


48  tt2  n^ 


*n(3/C)  n 
(10  ') 


(31) 


where  is  the  focusing  length  and  6 is  the  depth  of  the 

small  scale  perturbation.  Suydam  also  shows  that  a.n  optimum 

transverse  spatial  frequency  perturbation  for  small-scale  self- 

focusing  exists.  The  spatial  modulation  dimension  > that 

m 
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produces  the  greatest  self-focusing  is 


10 


-7 


12  tt  n2  I 


(32) 


To  apply  small-scale  self-focusing  theory  to  LiNbOg  we 
need  an  estimate  for  the  magnitude  of  n2  . We  are  unaware  of 
any  measured  value  for  LiNbOg  n2  and,  therefore,  choose  to  use 
the  approximate  expression 


n2(10  13  esu)  = 391 


34 

based  on  the  empirical  relation  of  Boling  et.al. 

Here  is  the  Abbe  number,  = (n^  - l)/n^  - n£)  where 

n^  and  nc  are  the  crystal  indices  of  refraction  at  the  wavelength 

0,486,  0.588  and  0.656  pm  . Using  the  above  relation  we  find 

-13 

n2(LiNb02)  = 11.1  x 10  esu.  This  value  may  be  an  underestimate 
of  the  true  value  due  to  the  effect  of  the  d states  of  the 
niobium  ion  and  to  the  application  of  the  empirical  model  developed 
for  glasses  to  a high  index  of  refraction  crystal. 

For  LiNbOg  we  find  from  Eq.(30)  that  the  critical  power 

for  whole  beam  self  focusing  is  P = 100  kW.  This  is  comparable 

c 

to  critical  powers  for  materials  with  large  values  of  n2  , such 
as  CSg  . 
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If  we  set  7,^  equal  to  the  crystal  length  of  5 cm, 

(10  cm  for  the  DSRO) , then  for  a perturbation  of  order  6=1 

we  find  from  Eq.(31)  that  Id  = 290  MW/cm2  (Id  = 145  MW/cm2  DSRO) 

for  the  SRO.  If  we  take  6 = .15,  as  for  a defect  in  an  anti- 

o 

reflection  coated  LiNbO^  surface,  then  Id  = 800  MW/cm  for  an 

2 

SRO  and  400  MW/cm  for  a DSRO.  The  optimum  transverse  modulation 

dimension  is  found  from  Eq.(32)  to  be  ) = 0.5  mm  for 

m 

2 

I = 300  MW/cm  . This  dimension  is  approximately  the  spacing  of 

crystal  Schlieren  defects  and  the  Fresnel  intensity  modulation 

of  a near  field  unstable  resonator  beam. 

Since  n^  for  LiNbO^  and  the  detailed  character  of  the 

perturbations  are  not  known,  the  damage  values  presented  here  are 

only  approximate.  One  measured  data  point  gave  a DSRO  damage 

2 

track  at  Id  = 300  MW/cm  . This  value  is  shown  in  Fig.  9b.  Also 
shown  in  Fig.  9b  is  the  intensity  inside  an  average  LiNbOg  crystal 
operating  at  three  times  threshold  for  an  8 cm  cavity  SRO.  Clearly, 
for  pump  pulses  less  than  30  nsec  duration  the  danger  of  self- 
focusing  is  critical.  This  is  particularly  the  case  for  crystals 
with  Schlieren  defects  or  for  beam  perturbations  produced  by 
imperfect  coatings  or  damaged  optical  components  prior  to  the  crystal 
An  additional  factor  may  be  introduced  by  partial  mode  locking  or 
axial  mode  modulation  in  the  pump  beam.  The  resultant  instantaneous 
intensity  values  may  easily  exceed  twice  the  average  peak  intensity. 

Figure  9 clearly  shows  that  longer  pump  pulsewidths  reduce 
the  0P0  intensity  threshold  and  thus  reduce  the  probability  of  self- 
focusing.  Experimentally  we  have  found  that  for  pulsewidths  longer 
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than  15  nsec  self-focusing  is  avoided  even  in  the  case  of 
severe  surface  damage  to  the  crystal.  The  longer  O-switched 
pulse  lengths  were  obtained  by  increasing  the  Nd:YAG  optical 
cavity  length.  Added  benefits  of  the  longer  pump  pulse  are  a 
decrease  in  the  operating  linewidth  of  the  0P0  due  to  the 
increased  number  of  passes  during  buildup  and  a greater  in- 
sensitivity of  intensity  threshold  on  0P0  cavity  length. 

In  summary,  Fig.  9 shows  that  any  increase  in  0P0  threshold 
brings  the  device  closer  to  the  surface  damage  or  small-scale 
self-focusing  damage  limit.  One  may  either  rely  on  the  RiOg  anti- 
reflection coating  for  increased  surface  damage  fluence  or  operate 

2 

at  a fluence  less  than  the  2.7  J/cm  set  by  the  bare  LiNbO^ 
crystal  surface.  In  any  case,  care  should  be  taken  to  minimize 
the  energy  fluence  threshold  by  the  proper  choice  of  0P0  para- 
meters and  by  the  use  of  a high  quality  pump  beam  of  proper  temporal 
durat ion . 

IV.  LINF.WIDTII  CONTROL 

A.  Introduction 

The  gain  linewidth  of  the  parametric  oscillator  is  set  by 
2 3 

crystal  dispersion.  ’ For  the  angle  tuned  LiNbO^  SRO  the  crystal 
bandwidth  varies  from  10  cm  * to  approximately  200  cm  * at 
degeneracy  as  shown  in  Fig.  2.  Therefore,  for  most  applications 
some  linewidth  reduction  is  desired.  Theory  shows  that  the  SRO 
is  capable  of  operating  in  a single  axial  mode  without  reduction 
of  conversion  efficiency  at  the  resonated  wave.  The  non-resonated  wave 
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reflects  the  axial  mode  spectrum  of  the  pump  field.  Single 
axial  mode  operation  has  been  achieved  earlier  for  the  Nd.YAG 
pumped  non-critically  phasematched  LiNbO^  SRO  using  a single 
high  resolution  etalon  line  narrowing  element. 

The  present  angle  tuned  LiNbO^  SRO  requires  a primary  line- 

width  control  element  and  an  etalon  for  narrower  linewidth  operation. 

Primary  linewidth  control  elements  that  have  been  tried  include  a 

thin  tilted  etalon,  a birefringent  filter,  and  a diffraction 
7 37 

grating.  ’ The  most  attractive  method  of  primary  linewidth 
control  is  a grating  proceeded  by  beam  expansion  prisms  as  shown 
in  Fig.  10.  An  analysis  of  the  prism  expander-grating  linewidth 
control  of  the  0P0  is  presented  next.  In  Section  C experimental 
results  are  described  including  stable  single  axial  mode  operation 
of  the  LiNb03  SRO. 

B.  Linewidth  Model 

In  this  section  a model  is  developed  that  describes  the 
SRO  linewidth  when  operated  with  a resonant  cavity  consisting 
of  beam  expansion  prisms  and  a grating.  This  method  of  primary 
linewidth  control  was  selected  because  of  linear  tuning,  a wide 
tuning  range  and  ease  of  implementation  at  a low  cost. 

The  0P0  linewidth  is  the  result  of  the  crystal  and  grating 
component  lineshape  functions  after  consideration  for  multiple 
passes  in  the  optical  cavity.  The  resultant  linewidth  Av  is 
calculated  from  the  product  of  the  individual  multiple-pass 
lineshapes  of  the  crystal  and  of  the  grating.  Since  the  lineshapes 
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are  approximately  Gaussian,  we  may  write  the  linewidth  as 


- - / 1 , 1 

Av  Y fAvc(p)  ]2  [Av^(p)]2 


(33) 


where  Avc  and  Av^  are  the  crystal  and  grating  linewidths  and 

p is  the  number  of  cavity  transits. 

The  lineshape  function  for  the  crystal  bandwidth  Avfi  is 

2 

given  by  evaluating  the  sine  (Ak£/2)  function  for  a phase  mis- 
match of  Ak  and  a crystal  length  of  l . An  analysis  similar 
to  that  for  the  diffraction  grating  presented  in  Appendix  D shows 
that  for  p passes  the  crystal  linewidth  reduces  to 

1 

Av  (p)  = — Av  (34) 

C /p  C 


where  Av  is  determined  by  the  crystal  dispersion  between 
signal  and  idler  wavelengths  and  in  cm"'*  is  approximately 


0 . is  defined  by 
si 


si 


A v = 
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(35) 


(30) 


and  is  approximately  equal  to  2Ang^ 
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This  basic  crystal  linewidth  assumes  a collimated 

pump.  For  divergent  pump  beams  the  linewidth  increases 
significantly.  An  expansion  of  the  wave  vector  mismatch  for 

12 

small  deviations  of  the  pump  input  angle  and  signal  frequency 
gives  the  linewidth  increase  due  to  pump  divergence  6(Av)  of 


approximately 


(37) 


where  a is  the  pump  beam  cone  angle.  For  LiNbO^  0P0  signal 
and  idler  frequencies  away  from  degeneracy,  B . is  approximately 

S X 

0.1.  For  p = 0.03  rad  and  } =1.06  ym  , the  linewidth 

p 

broadens  nearly  6 cm  * per  milliradian  of  pump  divergence.  The 

crystal  bandwidth  is  then  significantly  broadened  by  a diverging 

incident  pump.  If  line  narrowing  elements  are  used,  an  increase 

in  threshold  may  be  observed  because  part  of  the  generated  signal 

wave  is  removed  from  the  gain  medium  by  cavity  frequency  selection. 

Consider  the  0P0  cavity  shown  in  Fig.  10  which  uses  a 

Littrow  configuration  grating  at  resonant  wavelength  > and 

s 

grating  angle  0 relative  to  normal.  For  added  resolution  a 

beam  expansion  element  of  linear  magnification  M and  angular 

magnification  M is  included. 

The  prism  beam  expander  has  been  described  earlier  as  a 

means  of  beam  expansion  within  a dye  laser  cavity.*  Klauminzer 

has  discussed  its  advantages .which  include  ease  of  construction 

39 

and  insertion  into  the  optical  cavity. 
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Using  LiNbOg  as  the  prism  material,  we  designed  and 
constructed  a two  element  prism  expander  as  discussed  in 
Appendix  E.  The  magnification  of  a single  prism,  taken  as 
the  output  to  input  beam  width  ratio  W/w  is  given  by  Eq.(E5) 

M = tan  G (38) 

For  N prisms  in  sequence  the  net  magnification  is  the  product 

M = Mn  as  expected. 

The  angular  magnification  of  a Brewster  angle  prism  is 
shown  in  Appendix  E to  be 


where  n is  the  prism  index  of  refraction.  For  comparison,  a 
telescope  expander  has  an  angular  magnification  of 


/ 

I 


- 112  - 


with  respect  to  0 


which  yields,  for  the  Littrow  case 


out 


d dv  d0 

s _ s _ 

>s  vs  2 tan0 


(42) 


Here,  d is  the  groove  spacing  and  m is  the  order  number.  The 
angular  width  A0  is  calculated  in  Appendix  D for  the  case  of 
a multiple  passed  grating  with  an  optical  length  2L  between 
passes.  The  calculation  shows  that  diffraction  not  only  sets 
the  single  pass  angular  width  but  also  puts  a lower  bound  on  the 
linewidth  attainable  by  multiple  passing.  Using  Eq.(D.12)  we 
write  the  diffraction  limited  full  width  at  half  maximum  (FWHM) 
grating  linewidth  as 
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(43) 


s 


: g 

where  the  ,f  actor  /2£n2  converts  the  halfwidth  at  l/eJ  line- 

width  to  the  FV. I iM  value.. 


The  grating  resolution  may  be  broader  than  that  specified 
in  F.q.(43)  when  the  finite  aperture  set  by  the  gain  medium  size 
is  considered.  Wave  packets  outside  the  diffraction  limited  band- 
width may  oscillate  by  propagating  off-angle  in  the  cavity.  These 

waves  may  experience  gain  over  the  full  pump  beam  area.  Therefore, 

2 

to  be  consistent  with  the  1 /o  angular  width  used  to  obtain 
Fq (43)  we  take  the  proper  pump  aperture  diameter  to  be  2 w^  . 
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+ A X which  is 


Consider  a light  ray  of  wavelength  Xq 
initially  propagating  parallel  to  the  cavity  axis.  It  is 
shown  in  Appendix  D that  after  its  q'th  reflection  from  the 
grating,  its  angle  to  the  axis  is 


A 


a 


MAX 


d cos  0 


(44) 


After  a beam  expansion  device,  this  angle  is  M A . If  we 

ci  M 

note  that  the  gain  region  establishes  an  effective  aperture  of 
width  equal  to  the  pump  spot  diameter  2w^,  then  the  number  of 
passes  p in  which  the  ray  travels  across  the  aperture  is 
determined  by 

2w  = 2DM  A,  £ Q (45) 

P a Q=1 


where  the  sum  has  the  value  p(p+l)/2  and  D is  an  effective 
cavity  length.  In  air,  D would  simply  be  the  physical  cavity 
length.  Here,  however,  the  presence  of  a crystal,  of  index  n 

o 

- 2 
and  a beam  expansion  device  of  angular  expansion  M (M  for  prisms, 

a 

M for  a telescope)  gives 

D - Vs  + 'a  + WMa  <46> 

where  the  crystal,  air  and  expander  lengths  are  given  by  i , 

H and  respectively.  In  practice,  p is  fixed  by  the 

a 

cavity  buildup  time.  The  contribution  to  grating  linewidth  due 
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to  cavity  geometry  is  then 


Av  = /2£n2 

cav 


2 w cos  0 
P 


m p(p  + 1)  D Aq  Ma 


(47) 


The  full  grating 
aperture  width  t 

Av  (p)  = /2£n2 

g 


1 inewidth 
erms  given 

1 

2 tan  0 


is  thus  a sum  of  the  diffraction 
by 
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m p(p  + 1)  MaD  Ag 


(48) 


We  now  consider  the  relative  magnitudes  of  the  diffraction 

and  aperture  width  terms  in  Fq.(48)  for  typical  cavity  configurations, 

in  an  effort  to  understand  which  parameters  are  important  to  the 

design  of  a narrowband  QPO.  If  we  let  w = 0.2  cm,  w = 0.06  cm, 

P s 

Ag  = 1.7  pm  , i = 6000  lines/cm,  D = 8 cm  and  m = 1 , we  find 
that  for  a 10  ns  pump  pulse  and  no  expansion,  the  inclusion  of 
the  aperture  term  increases  the  linewidth  by  30".  For  a 20  ns 
long  pump  pulse  the  aperture  term  contribution  to  the  linewidth 
is  reduced  by  multiple  passing  to  8"  . If  beam  expansion  prisms 
are  used,  the  effect  of  the  aperture  term  becomes  negligible. 
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From  Eq.(48)  it  is  clear  that  to  improve  linewiclth  performance 
we  may  increase  grating  0 and/or  increase  the  beam  diameter 
on  the  grating,  2 M w . However,  each  approach  has  limits. 

b 

As  0 increases,  the  free  spectral  range  of  the  grating  decreases 
and  does  not  allow  0P0  operation  over  the  full  tuning  range  from 
1.4  ym  to  2.1  pm  . Increasing  M requires  a 4 prism  beam 
expander  and  a larger  grating  and  has  the  disadvantage  of  increasing 
threshold  due  to  the  added  cavity  length.  However,  for  sufficiently 
long  pump  pulsewidths  near  20  ns  the  increase  in  threshold  may 
be  tolerable. 

Various  techniques  exist  for  additional  linewidth  control. 

Single  frequency  operation  using  interferometer  arrangements 

have  been  demonstrated  for  lasers^  and  for  parametric  oscillators.2"* 

41 

Internal  tilted  etalons  have  also  been  used  successfully,  but 
have  suffered  from  insertion  losses,  a problem  which  is  alleviated 
by  the  use  of  beam  expansion  prisms.  Due  to  their  ease  of 
alignment  we  have  chosen  to  use  tilted  etalons  to  attain  narrow- 
band  operation . 

The  lineshapc  function  fe(Av)  for  a transmitted  beam  of 
a tilted  etalon  is  Lorentzian  in  the  frequency  deviation  from 
resonance,  Av  = v - v , for  Av  <<  v , 

D O 

fE(Av)  = [1  + 4 .^2(Av)2/(AvpsR)2]“1  (49) 
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is  the  free  spectral 


where  $ is  the  etalon  finesse  and  Av^^ 
range  2nt  for  a material  of  thickness  t and  index  n . For 
p passes  through  the  etalon,  the  linewidth  reduces  to 


Av„(l) 

Av  (p)  = — * (50) 

/p 

-2 

Clearly  Eq.(33)  could  be  extended  by  adding  (Av,-,(p)  to  the 

h 

sum  of  linewidth  factors  under  the  square  root. 

For  the  present  application  as  an  internal  cavity  element 

with  short  optical  pulses,  various  loss  mechanisms  arise.  The 

insertion  loss  y for  a tilted  etalon  in  a cw  laser  cavity 

42 

has  been  treated  by  Leeb  and  found  to  be 


2R 


2t  9 


YF,= 


E 


(1  - R)‘ 


n w 


o 


(51) 


where  R is  the  power  reflectivity,  0^  is  the  tilt  angle  and 

w0  is  the  spot  size  in  the  direction  perpendicular  to  the 

rotation  axis.  . v,'hen  the  prism  expander  is  used,  w is  large 

and  Yp.<  10%  even  for  a & = 20  etalon. 

A more  important  loss  mechanism  for  short  pulse  applications 

is  due  to  the  delay  propagation  time  through  the  etalon  as 

43 

considered  by  Daussy  et . al . They  find  that  a delay  near  1 ns 
may  occur  for  high  finesse  etalons.  This  can  amount  to  a considerable 
loss  during  the  buildup  time,  increasing  0P0  threshold  prohibitively. 
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Another  finding  is  that  the  effective  pulsed  finesse  may 
be  significantly  (~  1/3)  lower  than  that  calculated  for  a 
steady  state  condition.  These  results  are  consistent  with 
our  experimental  linewidth  results  for  the  pulsed  LiNbO^  OPO 
operating  with  a tilted  etalon  line-narrowing  element. 

C.  Experimental  Results 

The  linewidths  for  a number  of  OPO  cavity  configurations 
(see  Appendix  F) , have  been  measured  and  compared  with  theory. 

In  general  the  agreement  is  good.  The  observed  discrepancies 
can  often  be  attributed  to  factors  such  as  crystal  quality  or 
optical  alignment.  As  in  the  threshold  measurements  described 
in  Part  II,  we  have  found  that  crystal  optical  uniformity  is  an 
important  factor  and  that  the  presence  of  Schlieren  variations 
in  the  index  of  refraction  of  a crystal  increase  the  measured 
linewidth  up  to  a factor  of  two. 

Optical  alignment  of  the  cavity  and  of  pump  beam  collimation 
telescopes  is  also  more  important  for  narrow  linewidths  than  for 
threshold  minimization.  Some  care  must  be  exercised  in  the 
centering  of  the  grating  resonance  upon  the  peak  of  the  crystal 
bandwidth.  A lineshape  asymmetry  caused  by  a broadening  of  the 
high  frequency  side  of  the  line  occurs  if  this  condition  is 
violated.  If  the  grating  is  tuned  off  line  center,  the  crystal 
phasematching  conditon  is  broken  more  slowly  on  the  high  frequency 
side,  since  the  OPO  waves  may  then  angle  phasematch  and  propagate 
noncol linearly . On  the  low  frequency  side  Ak  >0  so  that  the 
OPO  waves  remain  col  linear  but  are  reduced  in  magnitude.  The 
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amount  of  broadening  may  be  calculated  by  expanding  the  wave- 
vector  mismatch  allowing  for  angle  phasematching  and  is  found 
to  be  quadratic  in  the  angle  between  signal  and  pump  beams. 

For  minimum  linewidths  and  col linear  beams  it  is  then  recommended 
that  the  grating  be  tuned  slightly  to  the  low  frequency  side  of 
the  crystal  bandwidth. 

The  measured  and  calculated  linewidths  for  OPO  cavities 
with  increasing  spectral  selectivity  are  given  in  Table  I.  The 
first  block  of  entries  shows  that  a prism  power  of  10  is  required 
to  give  wavelength  control  to  the  grating.  The  resulting  grating 
linewidth  is  0.53  cm  * . A typical  linewidth  trace  taken  using 
a scanning  1 meter  spectrometer  is  shown  in  Fig.  11. 

It  is  experimentally  observed  that  the  OPO  output  beam  at 

2-3  times  threshold  has  the  same  spot  size  as  the  incident  pump 

beam.  This  is  explained  by  noting  that  the  confinement  of  the 

signal  beam  is  relaxed  as  the  peak  of  the  pump  is  depleted,  and 

the  signal  spot  size  may  increase  to  approach  the  incident  pump 

spot  size.  It  might  be  argued  that  the  proper  w to  use  in 

s 

Eq.(48)  is  this  larger  spot.  However,  the  linewidth  is  determined 
during  buildup  before  depiction  occurs  and  is,  therefore,  set  by 
the  smaller  confined  signal  spot  size  as  calculated  in  Appendix  A. 

If  a tilted  etalon  is  inserted  into  the  cavity,  the  OPO 
linewidth  is  reduced  as  shown  in  the  second  group  of  entries  in 
Table  I.  As  a rule  of  thumb,  the  etalon  free  spectral  range  is 
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TABLE  I 


OPO  LINEWIDTHS  : THEORY  AND  EXPERIMENT 

PARAMETERS  : } * 1.6  ».  w = .2  cm,  t = 10  ns, 

S r ' p r p 

L = 20  cm,  D = 8 cm,  N = 6000  lines/cm 


Configuration  Linewidth  (cm  ^ ) Comments 


Theory  Expt 


Grating  only 

2.2 

2.2 

w = . 1 5 cm 

Crystal  has  wavelength  contro 

Grating,  x5  prisms 

.86 

.73-1.2 

Grating  & crystal  both  contro 

Grating,  xlO  prisms 

.40 

.53 

Grating  has  X control 

Grating, x5  prisms 
+ 2 mm,  F = 7 etalons 

.059 

< .1 

2 etalon  modes 

+ 1 mm,  F = 7 etalons 

.15 

.14 

1 etalon  mode 

+ .5  mm,F  = 7 etalons 

.3 

.63 

1 etalon  mode 

Grating,  xlO  prisms 

.02 

.075 

20  ns  pump  pulse 

2 mm,  F = 7 etalons 

3 axial  modes  FWHM 

Grating,  xlO  prisms 

2 mm , F = 7 etalons 

< .02 

Single  axial  mode  operation 

Reflective  etalon 
(R  = 507.,  60%) 

- very  stable 
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chosen  to  be  twice  the  grating  linewidth.  A single  etalon 
mode  is  then  produced  with  the  linewidth  determined  by  the 
effective  multipass  finesse.  As  discussed  above,  the  insertion 
loss  of  the  etalon  increases  with  increasing  & imposing  a 
limit  on  the  linewidth  attainable  with  a single  etalon.  We 
have  found  experimentally  that  for  a finesse  greater  than  10 
the  insertion  loss  raises  the  0P0  threshold  excessively,  even 
when  the  etalon  is  used  in  the  prism  expanded  beam. 

In  a final  series  of  measurements  we  have  improved  the 
line-narrowed  0P0  performance  and  have  successfully  demonstrated 
stable  single  axial  mode  operation.  An  elliptical  20  ns  pump 
beam  was  incident  with  a 6 mm  major  diameter  oriented  in  the 
W’alkoff  plane  to  minimize  threshold  and  a 2 mm  minor  diameter 
set  as  the  aperture  for  grating  narrowing.  A reflective  etalon 
was  added  to  the  cavity  in  the  form  of  a pair  of  output  coupler 
mirrors  spaced  5 cm  apart.  In  this  experiment  the  non-optimum 
choice  of  50%  and  60%  reflecting  mirrors  was  used.  The  output 
linewidth  was  analyzed  by  diverging  the  0P0  output  through  an 
air  spaced  plane  parallel  Fabry-Perot  etalon  with  a finesse  of 
20  in  the  visible.  The  resulting  rings  wore  observed  with  an 
infrared  vidicon  and  photographed  from  the  monitor  display. 

With  the  output  reflective  etalon  mirror  misaligned  the  etalon 
ring  pattern  shown  in  Fig.  12a  was  obtained.  The  pattern  shows 
that  the  0P0  operated  at  a FWIIM  linewidth  of  0.075  cm  ^ , or  three 
axial  modes  spaced  0.025  cm  ^ apart.  Whenthe  outside  mirror  was 
aligned  to  form  a reflecting  etalon  the  linewidth  collapsed 


i 
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to  a single  axial  mode,  as  shown  in  Fig.  ]2b.  The  etalon 
effect  of  the  output  mirror  pair  was  verified  by  observing 
the  frequency  jumps  of  0.1  cm  ^ by  tilting  the  lower  resolution 
etalon  angle.  The  single  axial  mode  frequency  was  quite  stable 
and  repeatable  which  is  noteworthy  since  no  special  experimental 
precautions  were  taken  to  establish  an  interferometrically  stable 
0P0  cavity . Also , due  to  the  1 sec .persistence  time  of  the  vidicon, 
observation  of  clear  axial  modes  supports  the  claim  of  their 
frequency  stability.  The  0P0  is  a reactive  device  with  small 
thermal  loading  of  the  LiNbOg  crystal  and  optical  cavity  thus 
avoiding  resonator  instabilities.  To  the  resolution  limit  of 
the  analyzer  etalon  the  FWHM  of  a single  axial  mode  appears  to 
be  less  than  0.01  cm-^  or  300  MHz . The  expected  Fourier  transform 
limited  bandwidth  is  near  100  MHz.  The  measured  output  pulse 
energy  in  single  mode  operation  was  10  mJ,  which  is  approximately 
80%  of  the  un-narrowed  0P0  output  energy  under  similar  pumping 
conditions. 

V.  CONCLUSION 

This  paper  summarizes  the  parameters  relevant  to  the  design 
of  a reliable  LiNbOg  parametric  oscillator  capable  of  stable 
narrowband  operation.  The  prime  requirements  involve  the  quality 
of  both  the  LiNbOg  crystal  and  the  pump  laser  beam. 
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Imperfect  crystal  quality  is  a major  cause  of  increases  jn 
both  the  0P0  threshold  and  linewidth.  The  oscillation  threshold 
for  an  average  LiNbOg  crystal  is  70%  higher  than  that  predicted 
for  an  ideal  crystal , Linewidths  may  be  increased  by  a factor 
of  two  in  crystals  with  Schlieren  defects.  Pump  laser  phase- 
front  quality  is  also  important  and  high  quality  optical  components 
are  essential.  For  example,  the  use  of  commerical  grade  lenses 
in  a telescope  may  increase  the  threshold  by  60%. 

The  operating  limits  of  the  0P0  are  set  by  the  onset  of 

crystal  damage.  There  are  two  types  of  damage,  and  each  imposes 

a constraint  upon  the  configuration  of  the  pump  beam.  Internal 

crystal  tracking  due  to  small  scale  self-focusing  scales  as  the 

peak  power  of  the  pump  pulse.  Pulse  lengths  greater  than  10  ns 

are  essential  to  avoid  self-focusing,  and  pulses  longer  than  15  ns 

are  highly  recommended.  In  addition,  the  surface  of  LiFbOg  may 

2 

be  damaged  by  a peak  energy  fluence  of  2.7  J/cm  for  uncoated 
LiNbOg.  The  damage  limit  is  somewhat  higher  for  an  Si02  anti- 
reflection coated  surface. 

We  have  designed  the  LiNbOg  0P0  for  a minimum  threshold 
fluence  to  provide  for  reliable  damage  free  operation.  The 
threshold  increase  due  to  pump  beam  walkoff  leads  to  a minimum 
pump  spot  diameter  of  4 mm  and  crystal  length  of  5 - 6 cm. 

In  addition  the  cavity  length  should  be  minimized.  A complete 
cavi ty , includi ng  grating,  prism  ten-power  expander  and  tilted 
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etalon  may  be  constructed  with  a 16  cm  physical  length.  Since 
the  effective  loss  due  to  buildup  time  dominates  other  loss 
mechanisms,  the  threshold  is  not  a strong  function  of  output 
coupler  reflectivity  R.  For  greater  useful  energy  conversion 
it  is  necessary  that  the  output  coupler  provide  the  dominant 
cavity  loss.  Therefore,  we  choose  to  operate  with  an  output 
reflectivity  of  less  than  50%. 

The  required  pump  energy  to  operate  the  LiNbO^  SRO  varies 
with  crystal  and  pump  laser,  but  lies  in  the  100-150  mJ  range. 

As  an  added  note,  we  have  operated  a considerable  number  of 
0P0  cavity  configurations  as  summarized  in  Appendix  F.  However, 
the  best  results  to  date  have  been  obtained  with  the  "L"  cavity 
which  utilizes  the  input  beam  splitter  as  shown  in  Fig.  10. 

Narrow  linewidths  have  been  achieved  in  the  LiNbO  0P0 

o 

with  the  use  of  a 600  l /mm  diffraction  grating  in  first  order 
with  the  addition  of  alOx  prism  beam  expander.  The  measured 
linewidths  were  0.5  cm-^ . A 2 mm  tilted  etalon  with  & = 7 
reduces  this  linewidth  to  0.08  cm”\  corresponding  to  operation 
with  3 axial  modes.  A reflective  etalon  output  coupler  con- 
figuration further  collapses  t)  • output  to  a single  axial  mode 
which  is  remarkably  stable  in  frequency. 

The  angle  tuned  LiNbO^  parametric  oscillator  is  a useful 
laboratory  tool  for  studies  which  require  a broadly  tunable 
source  of  narrowband  infrared  radiation. 
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APPENDIX  A 


RESONATED  WAVE  SPOT  SIZE 


For  a Gaussian  profile  Rain  medium,  the  signal  spot 
size  evolves  during  buildup  to  threshold  to  a radius  determined 
by  the  pump  spot  size  and  diffraction.  We  calculate  in  this 
appendix  the  resonated  signal  wave  steady  state  spot  size  and 
determine  whether  this  condition  is  indeed  reached  during  the 
0P0  buildup  time.  We  proceed  by  considering  the  signal  spot 
size  as  it  is  iteratively  narrowed  by  the  gain  polarization 
and  broadened  by  diffraction.  The  balance  determines  the  final 
signal  wave  spot  size. 

Let  the  fields  have  Gaussian  profiles, 


2.2 

E.  = E.  e~r  /wj 
3 3 o 


(A.l) 


From  Eq.(l)  we  see  that  the  driving  polarization  for  the  0P0 
waves  have  radii  w given  by 


W 


2 

wi 


w 


P 


-2 

wi 


+ -o  (A. 2) 


w. 


The  idler  wave  is  free  and  therefore,  assumes  the  profile  of 
its  driving  polariztion,  or  w = w . The  signal  polarization 
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radius  may  now  be  written  as 


(A. 3) 


The  radius  of  this  wave  when  allowed  to  propagate  a cavity 
round  trip  is 

- 2 -2  / 2L  > 
ws  " ws  1 + [~ 

We  let  w'  ~ w and  iterate  to  see  how  soon  a stable  signal  mode 
s s 

radius  is  established.  Numerical  iteration  shows  that  a spot 
size  within  30%  of  steady-state  is  established  in  only  2-3  cavity 
round  trips.  This  verifies  the  assumption  made  in  the  model. 

The  steady  state  signal  spot  size  is  found  directly  by 
again  letting  w'  = w and  substituting  (A.d)  into  (A. 3).  The 

o S 

—2 

resulting  condition  is  a cubic  equation  in  w , 


2. 


2i 

s 


(A. 4) 


= 0 


(A. 5) 


From  this  calculation  we  find  that  the  signal  spot  size 
is  significantly  less  than  that  of  the  pump.  For  example,  a 
pump  spot  of  2 mm  radius  incident  on  a 15  cm  long  0P0  cavity 
operating  at  1.9  pm  produces  a signal  spot  size  of  0.62  mm. 
This  spot  size  is  important  for  later  linewidth  calculations. 
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Some  experimental  evidence  exists  to  support  the  prediction 
for  a small  confined  signal  spot  size.  In  an  initial  experiment 
to  measure  0P0  threshold,  a 6 mm  diameter  pump  beam  was  apertured 
to  a diameter  of  2 mm  with  no  change  in  threshold.  This  showed 
that  near  threshold  only  the  central  portion  of  the  pump  is 
useful  in  providing  signal  gain.  Next,  the  flat  0P0  cavity 
output  mirror  was  replaced  by  a curved  mirror  to  produce  a 
Gaiissian  mode  stable  cavity  with  a signal  spot  size  equal  to 
that  predicted  by  the  confined  signal  calculation  above.  No 
change  in  0P0  threshold  was  observed  lending  support  to  the  model 
presented  here. 

It  should  be  noted  that  for  higher  pump  powers  the  Gaussian 
profile  assumed  in  this  calculation  becomes  flattened  and 
distorted  due  to  depletion.  This  relaxes  the  constraint  on  the 
signal  spot  size  which  may  expand  to  match  that  of  the  pump. 
However,  the  spot  size;  calculated  here  should  be  valid  during 
the  0P0  buildup  time  and  at  operation  near  threshold. 
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APPENDIX  B 

SRO  - Single  Pass  Gain  Solution  Including  Poynting  Vector 

Walkof f 

To  calculate  the  signal  and  idler  field  buildup  in  an  0P0 
resonator  some  care  must  be  taken  to  consider  only  that  fraction 
of  the  generated  signal  wave  which  couples  into  the  resonator. 

In  addition,  the  coupling  of  the  pump  and  generated  waves  is 
gradually  broken  due  to  pump  wave  Poynting  vector  walkof f.  It 

is  convenient  to  treat  both  effects  in  the  coupling  coefficient 

14 

formalism  of  Kogelnik. 

The  radii  of  the  polarization  fields  for  the  signal  and 
idler  waves,  w and  w^  , are  found  in  Appendix  A to  be 

w.  = Wi  \ = IX/ vg  + 2 / \Vp ] “ (B.l) 


Because  the  signal  wave  is  resonant  its  polarization  profile 

must  be  expanded  in  terms  of  cavity  modes. 

Only  the  TEM  mode  is  of  interest  here.  Let  P^  represent 
oo  .1 

the  Gaussian  profile  of  a polarization  which  drives  the  signal 
or  idler  wave.  If  we  describe  the  fields  in  transverse  coordinates 
x and  y the  coupling  coefficient  can  be  written  as 


c 


oo 


00  * 

oo 

/ P.Cx)  E . ( x)  dx 

J B-j(y)  Ei(y)  dy 

OO 
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(B.7) 
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The  polarization  terms  of  Eq . ( 3 ) can  now  be  re-written  in  terms 
of  cavity  mode  field, 


dE, 


dz 


+ “sEs  = J KscsEpEi  exP 


T(  Z 


JAkz 


4 £. 


w 


(B.8) 


dE. 


— - + a.E.  = j k.E  E exp 

dz  11  1 P s 


ir  z 


Z 2 
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JAkz 


(B . 9) 


Let 


Then 


E_  = E;  e“asz  , E.  = E'  e"aiz 

O 11 


(B . 10) 
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where  ag  = ai  = a is  assumed.  Let  Ak  = 0,  since  we  are 
interested  in  a phasematched  0P0.  Then  substitution  yields, 
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dz 


exp 
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( B . 1 2 ) 
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Changing  variables  by  letting 


yields 


d E. 


K 


K . 
S 3 


The  solution  satisfying  the  boundary  condition  that  Y' 

s 

and  EC  =0  at  z'  = 0 is, 


E (z)  = E (o)  e asz  cosli  Y 7,' 

b o 


where 


r = •• 


with  vc  and  I defined  by 
P J 


2 u>  a),  d"  „ T . ip  1 2 

S3,  off  1 = J nee  E 

3-  , p o'  p1 

n n . n e c 
s 3 p o 


The  idler  wave  has  the  solution, 


M*>  - K.<o)  e-v-  fh-  ♦ 3*)  slnh  rz- 


g . K 
f 3 S 


where  ij>  is  the  plia.se  of  the  pump  field. 


(13.13) 


(B.14) 


= Es(o) 


(B . 15) 


( B . 1 6 ) 


(B.37) 
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APPENDIX  C 


Effective  Nonlinear  Coefficient  for  I.iNbO^ 


For  LiNbOg  with  a point  group  symmetry  of  3 m,  the  effective 

44 

nonlinear  coefficient  is 

deff  = d31  sin  9 + d22  cos  9 (C.l) 

where  © is  the  angle  between  the  optic  axis  and  pump  propagation 
direction.  Also,  d22  = .53  d31.45  So  for  © = 47°  , in  the  middle 
of  the  0P0  tuning  range,  we  have 

deff  = 1,09  d31  (C.2) 

From  1.06  pm  second  harmonic  generation  measurements  the 

recommended  value  of  d^^  given  by  Choy  is  5.95  x 10  ' m/V. 

We  scale  this  to  the  0P0  wavelength  range,  by  using  the  dispersion! c.c 

Miller's  Delta  A.  ..  value  of  1.13  x 10  ^ , where 

i.l  k 


APPENDIX  D 


V 

I 


Multiple  Pass  Diffraction  (Irating  Resolution 

In  this  Appendix  we  consider  the  multiple  passed  reflection 
of  monochromatic  light  from  a detraction  grating  and  calculate 
the  resulting  angular  halfwidth  for  the  case  in  which  a propagation 
length  2L  is  allowed  between  grating  reflection.  This  distance 
permits  diffraction  to  again  broaden  the  grating-narrowed  diverg- 
ence angle  and  to  establish  a steady-state  ABp  independent  of 
further  grating  reflections. 

Neglect  for  the  present  the  effect  of  L and  consider  the 

normalized  grating  intensity  function  for  infinitely  narrow 

. 48 

grooves  as  given  by 

. 2 * 
sm  .A'  Y 

1 = “5 p (D.l) 

N sin  y 

where  N is  the  number  of  illuminated  grating  grooves  and  y 
is  defined  by 

ltd 

y = — [sin  0q  + sin(©o  + A9)  ] (D.2) 

Here  d is  the  groove  spacing,  0 is  the  incident  angle 
relative  to  grating  normal,  and  AB  is  the  deviation  of  a reflected 
light  ray  from  Lit trow  configuration.  The  grating  resonance  is 
defined  by 

-y-  sin  0Q  = nir  (D.3) 
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whore  n is  an  integer.  If  A0  <<  1,  wo  then  find  that 


Y 


nd 

X 


cos  0 A 0 
o 


(D.4) 


After  p reflections  the  grating  intensity  function  becomes 


I(P) 


• 2 , 
sin  .yt'y 


Ur  Y 


(D.5) 


We  define  the  angular  halfwidth  of  the  grating  to  be  the  angular 

deviation  from  resonance  A0^  which  reduces  the  grating  intensity 
2 

function  to  1/e  of  its  maximum.  A numerical  calculation 
using  Fq.(D.5)  shows  that  to  a good  approximation 

2 

= — CD.  6) 

S 

for  p >>  1 . The  angular  halfwidth  of  a grating  narrowed  light 
packet  is  then,  from  Eq.(D.d)  and  (D.6), 


A©d  = ( X/tt\V)  l//p 


( D . 7 ) 


2 

where  W is  the  1/e  spot  radius  and  ~ 2W/d  has  been  used. 
The  far  fieid  diffraction  limited  divergence  angle  is  then  reduced 
by  l//p"  for  a multi- passed  grating. 
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The  propagation  length  2L  is  now  included  to  calculate 
the  effect  of  diffraction.  For  a Gaussian  beam,  the  divergence 
angle  0.  is  W(z)/R(z)  where  w(z)  is  the  soot  size  and  R(z) 

V. 

a q 

is  the  beam  curvature.*"  we  have  then 


(D.8) 


2 

where  zR  is  the  Rayleigh  range  tt  Wq/  / . If  beam  expansion 

prisms  are  used,  this  angle  is  reduced  by  1/M  . Since  w = ^/M 

iX  s s 

we  obtain 


(D.  9) 


We  are  interested  in  the  steady  state  conditon  which  must 
occur  for  some  p'  . We  may  approximately  write, 


TUtf 


M/M, 


TFT 


1 + /l  + (zr/2L): 


X 1 

— (D.IO) 


u Ws  /p' 


where  L is  the  cavity  optical  length.  Assuming  p'  > 1 yields 


4 p'  (p'  + 1) 


r*V| 

v M i 


1 + 


'ZT 


,2L> 


2 n 
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(D.ll) 


Therefore 


/M  z ' 

p-  * (.A  _R 

\ M 4Li 
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If  Eq.(D.ll)  is  less  than  unity  we  set  p'  = 1 . For  w = 0.6  mm, 

s 

?s  =1.9  gm  , L = 15  cm  and  5 x prisms,  p'  = 3.  This  is  much  less  thr 

the  transits  possible  during  the  0P0  buildup  time , ct /2L-.  Thus  the 
same  diffractive  process  which  stabilizes  the  spot  size  evolution 
in  Appendix  A is  also  responsible  for  the  rapid  stabilization  of 
the  diffraction  limited  halfwidth  to  the  steady  state  value 


agd  - 
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We  also  need  to  understand  the  mechanism  of  the  ejection 

of  off-resonant  wavelengths  from  a grating  mirror  cavity  of 

center  wavelength  ) and  grating  angle  G . If  1 = 1 + A 1 

s o o 

the  grating  equation  becomes 


sin  9"  + sin  0"  = ml 

T 

where  d is  the  groove  spacing,  and  for  pass  n 


(9.14) 


e ' = e - a . , G~=e  + a 

o n-1  ’ o n 


where  9^  and  9"  are  the  incident  and  exit  angles  relative 
to  the  grating  normal.  Expansion  of  Eq.(I1.13)  yields 

mN ( A 1 ) 


A = A , + 
n n-1 


cos  9^ 


(9.14) 
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For  the  resonant  wavelength  Xq  , the  grating  acts  as  a simple 
mirror.  For  a wavelength  diviation  A 7 , the  grating  reflection 
angle,  relative  to  the  cavity  axis,  is  incremented  by  7/d  cos  9 

on  each  pass.  If  we  take  Aq  = 0 , 


An  = 


m./f  A 7 


d cos  9 


o 


(D.  15) 
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APPENDIX  E 


Brewster  Angle  Beam  Expansion  Prisms 


Consider  the  two  prism  beam  expanders  shown  in  Fig.jo 
The  linear  magnification  of  a single  prism  is  given  by 


cos  9 

M = (E.l) 

cos  9^ 


where  9^  and  9g  are  the  angles  of  incidence  and  exit  from 

the  prism  first  surface.  The  maximum  magnification  occurs  when 

the  beam  exits  normal  to  the  prism  second  surface.  Then  9 = a 

o ’ 

where  a is  the  prism  apex  angle  and  Snell's  law  results  in 


1 tan  9a 

- • (E.2) 

n tan  a 


For  low  insertion  loss  the  incident  angle  is  set  at  Brewster's 
angle  so  that 


n = tan  9..  (E.3) 

and 

1 

M = = n (E.d) 

tan  a 


Therefore,  near  Brewster's  angle 


M = tan  9i 


(E.5) 


- 1 30  - 


The  angular  magn i fi.cat ion  for  a prism  beam  expander  is 


found  by  differentiating  Snell's  law 


6 0.  cos  0.  = n cos  0„  6 0 
li  e e 


(E.C) 


so  that  . _ 

e 

6 0. 
l 

Using  Eqns.(E.l)  and  (E.d) 


1 cos  0^ 
n cos  0 

e 

we  have  for  a Brewster  angle 


(E.7) 

prism , 


Therefore , 


For  comparison, 


1 6 0 1 

_ e _ 

M 60.  nM 

a l 

M = Mn  = M2  ( E , 9 ) 

a v ' 

a telescope  has  an  angular  mangif ication  of 


(E.8) 


M = M (E. 10) 

The  minimum  optical  path  length  £q  through  a prism  is  given  by 

£ = n W tan  a = W (E.ll) 

The  two  jirisni  expander  is  nearly  achromatic.  The  angular 
deviation  AO  for  an  index  dispersion  An  is  to  first  order 

A0  = (l/n()  - 1/n^)  An  (K.12) 
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where  n is  the  mean  index  of  refraction  of  the  medium.  For 
o 

glass  prisms  the  angular  deviation  is  only  2.7  mrad  over  the 
entire  visible  spectrum. 

The  insertion  loss  of  the  prisms  is  given  by  the  P'resnel 
reflection  loss  at  the  first  surface  plus  the  anti-reflection 
coating  loss  at  the  near  normal  surface.  For  the  LiNbO^  prisms 
this  loss  is  less  than  1%  per  prism,  which  is  insignificant 
compared  to  output  coupling  losses.  Thus  the  LiNbO^  prisms 
beam  expander  is  readily  inserted  into  the  0P0  cavity  without 
significantly  increasing  threshold. 
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APPENDIX  F 


Experimental  Cavity  Conf igurations 

In  the  course  of  this  work,  a number  of  0P0  configurations 
have  been  constructed  in  an  effort  to  reduce  threshold  and  to 
improve  linewidth  performance.  We  have  found  that  the  Gaussian 
pump  profile  is  necessary  for  minimum  threshold  due  to  the 
longer  walkoff  length.  Attempts  to  improve  performance  have 
concentrated  on  a reduction  of  cavity  length,  as  shown  by  the  in- 
line arrangements  using  Fresnel  pump  coupling  (Fig.  14a)  and 
using  a dichroic  input/output  separation  mirror  (Fig.  14b). 

The  latter  is  necessarily  a DSRO  but  has  the  disadvantage  that 
for  moderate  pump  levels  the  signal  wave  feedback  from  the  pump 
high  reflector  mirror  is  sufficient  for  oscillation.  In  the 
Fresnel  coupled  cavity,  the  pump  is  reflected  off  the  first  prism 
surface  This  arrangement  works  well  but  is  somewhat  inconvenient 
to  alig.,  , id  does  not  efficiently  couple  the  pump  beam  into  the 
cavity.  The  L shaped  cavity  using  a dichroic  input  mirror 
shown  in  Fig.  10,  performs  bettor  than  any  other  configuration 
we  have  tested. 

For  linewidth  control  we  have  operated  cavities  with  a 
Galilean  telescope  expander  (Fig.  14c),  a mirror/grating  combination 
roar  reflector  (Fig.  14d)  and  a Brewster  angle  birefringent  tuning 
element  (Fig.  14o).  The  telescope  expander  is  difficult  to  align 
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and  to  collimate.  It  also  suffers  from  chromatic  aberration 
and  produces  back  reflected  focal  spots  which  may  damage  other 
components.  The  resonant  mirror/prating  setup  must  be  inter- 
ferometricly  stable  to  operate  properly.  Birefringent  tuning 
has  the  disadvantage  of  being  nonlinear.  We  find  that  the 
diffraction  grating  when  preceeded  by  a prism  beam  expander 
has  better  performance  than  any  other  configuration  tested. 
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FIGURE  CAPTIONS 


1.  Simplified  OPO  schematic -Mirror  is  highly  reflecting 
between  1.4  - 2.1  pm  • Output  coupler  has  signal 
reflectance  R.  For  DSRO  operation,  pump  high  reflector 
M„  may  be  used . 

O 

2.  Angle  tuned  LiNbO^  OPO  tuning  curve  and  crystal  gain 
bandwidth . 

3.  SRO  threshold  fluence  vs  cavity  physical  length  solid 

curve  shows  results  of  numerical  calculation.  Dashed  curve 
shows  results  of  square  pulse  model.  Z = 6 cm,  2 = 

4 mm , R = 60$  , t = 15  ns. 

# 

4.  SRO  threshold  fluence  vs  crystal  length.  L.'  = 8 cm, 

2 w = 4 mm,  R = 60$,  t = 15  ns. 

P 

5.  SRO  threshold  fluence  vs  signal  wave  reflectance.  Error 
bars  reflect  variations  due  to  differing  idler  reflectances. 
L'  = 8 crn,  Z - 6 cm,  2 w = 4 mm,  t = 15  ns. 

p 

6.  SRO  threshold  fluence  vs  pump  spot  diameter  L'  = 8 cm,  • 

Z = 6 cm,  R = 60%,  t = 15  ns. 

7.  SRO  threshold  fluence  vs  pump  pulscwidth.  Data  scaled 
by  0.37.  L 1=  8 era,  l = 6 cm,  2 w = 4 mm , R = 50%. 

o 

8.  SRO  energy  conversion  efficiency  vs  number  of  times  above 
threshold . 
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9.  Damage  limits  of  LiNbOg  vs  pump  pulsewidth.  a)  solid 
curve  shows  pump  fluence  at  3 x 0P0  threshold  for  an 
average  crystal.  Surface  damage  fluences  for  the  bare 
LiNbOg  surface,  for  the  uncoated  surface  in  an  0^  atmosphere, 
for  an  SiOg  coated  crystal  are  shown,  b)  solid  curve 

shows  pump  intensity  at  3 x 0P0  threshold  inside  an  average 
crystal.  An  approximate  experimental  self-focusing  damage 
limit  is  indicated  by  the  dashed  line. 

10.  Schematic  of  the  SRO  cavity  for  linewidth  studies.  Primary 
linewidth  control  is  provided  by  prisms  and  grating  combination. 
Line  narrowing  is  accomplished  with  a tilted  etalon.  A 
resonant  reflector  used  as  an  output  coupler  collapses  the 
linewidth  to  a single  axial  mode. 

11.  SRO  grating  linewidth  using  a ten  power  prism  beam  expander 
as  measured  using  a scanning  1 meter  spectrometer. 

12.  SRO  linewidth  measurements  using  Fabry-Perot  etalon.  0P0 
cavity  includes  a tilted  etalon  and  resonant  reflector, 
a)  resonant  reflector  is  misaligned,  producing  a FWIIM 
linewidth  of  .075  cm-^ , or  3 axial  modes,  b)  resonant 
reflector  is  aligned,  producing  stable  single  axial  mode 
operat ion . 

13.  Two  prism  expander  setup.  Incident  and  exit  angles  are 
shown  as  9^  and  ©e  . Prisms  apex  angle  is  a . Magnification 
is  given  by  W/w. 
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14.  OPO  cavity  configurations,  a)  pump  input  coupling  using 

Fresnel  reflection  of  prisms  P,  b)  Mg  is  pump  high  reflector 
dichroic  beamsplitter  BS  is  used  to  separate  signal  output, 
c)  telescope  T used  for  beam  expansion,  d)  partially 
reflecting  signal  mirror  Mg  to  act  with  grating  G as 
resonant  reflector,  e)  birefringent  filter  BF  used  as  line 
narrowing  element,  is  signal  high  reflector.  Other 

elements:  crystal  = XT,  output  coupler  = . 
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We  have  generated  50  mJ  of  lfi.9-am  radiation  by  stimulated  rotational  Raman  scattering  in  3 atm  of  H-.>  gas 
pumped  by  a CO..  TKA-laser  source.  Threshold  was  reached  by  injection  of  a few  microjoules  of  10.9-am  radiation 
generated  by  four-wave  mixing.  We  achieved  25%  peak  power,  or  40%  peak  photon  conversion  efficiency. 


In  1976  Byer  proposed  using  stimulated  rotational 
Raman  scattering  in  hydrogen  gas  as  an  efficient 
method  of  frequency  conversion  in  the  infrared.'  In 
particular,  it  was  noted  that  stimulated  Raman  scat- 
tering in  parn-Hj  gas  with  a C,0  > laser  pump  generates 
16-pm  radiation,  which  is  of  potential  use  for  U F<;  iso- 
tope enrichment.  It  was  subsequently  learned  that 
others  have  also  suggested  stimulated  rotational  Raman 
scattering  in  the  infrared.2  Recently  Frey  et  ai 3 have 
generated  16-pm  output  at  the  1 -mJ  level  by  vibrational 
Raman  scattering  in  H2  and  N2  pumped  by  a ruby- 
pumped  dye  laser. 

We  have  successfully  generated  50  mJ  of  16.9-pm 
iation  by  four-wave  mixing-assisted  stimulated  ro- 
tational Raman  scattering.  We  used  a 2-J,  70-nsec  CO< 
laser  source  to  pump  a 25-pass,  4-m-long  cell  filled  with 
3 atm  of  p-H2  gas  obtained  from  a liquid  H2  source. 
Based  on  theoretical  calculations'  and  rotational 
Raman  threshold,  scaled  from  careful  measurements 
at  1.004  pm,  we  predicted  a C02-pumped  rotational 
Raman  threshold  in  our  multipass  cell  of  1. 7-2.6  J in  a 
70-nsec  gain-switched  pulse.  Our  available  C02  pump 
energy  was  not  adequate  to  reach  stimulated  Raman 
threshold  by  amplification  of  spontaneous  emission. 
However,  when  16.9-pm  radiation  generated  by  four- 
wave  mixing  was  injected  into  the  multipass  cell,  the 
Raman  gain  was  sufficient  to  lead  to  significant  pump 
depletion  and  a peak  photon  conversion  efficiency  of 
40%  to  the  Stokes  wave  at  16.9  pm. 

The  four-wave  mixing  process  was  first  described  for 
interactions  in  crystals'  and  later  extended  to  gasesr>'fi 
as  a means  for  generating  widely  tunable  infrared  ra- 
diation.7 Sorokin  cl  al.*  have  demonstrated  four-wave 
mixing  in  H i gas  using  a CO*>  laser  source  as  a method 
for  16-pm  generation. 

The  equations  governing  both  four-wave  mixing  and 
stimulated  Raman  scattering  arc  given  by" 


^ = - -1'-  XR*< |ie„ I2  Ep  + Ua! 

02  2cn,, 

. u>„ 

2 cn, 


= + + E,*K„Kpe-iM').  (11>) 


on, 


' X it -I'd  + KpESR*-1**).  (lc) 

0 z Jc/i, 


0EU  _ _j_  M|) 

dz  2 cn„ 


XH-mr-Eg  + Ep*EsE,c^kn,  (Id! 


where  up  — ws  = <■>,*  — with  wp  — ws  = a>p  and  Ak  = 
— (kp  — ks)  + (ki  — kn)  with  a>n  the  characteristic 
Raman  frequency  of  the  medium  with  a peak  Raman 
susceptibility -xk".9  In  our  case  Ep  is  the  NdrYAG 
pump  field  at  1.064  pm,  Es  is  the  generated  Stokes  field 
at  1.10  pm,  E,  is  the  input  field  at  10.6  pm,  and  E„  is  the 
generated  output  at  16.95  pm. 

If  we  neglect  four-wave  mixing  for  the  moment,  Eqs. 
(la)— (Id)  separate  into  two  pairs  of  equations  that  de- 
scribe stimulated  Raman  scattering  pumped  at  Ep  and 
E,.  In  the  absence  of  pump  depletion,  the  Stokes  power 
increases  exponentially  as  /-*.,(/)  = Ps(0)  expO 7*0,  where 
I\(0)  and  l\(l)  are  the  input  and  output  Stokes  power, 
respectively,  and  is  the  Raman  power-gain  coefficient 
given  by 

_ <*»XR’\En\2  _ r\H"Ip  (2) 

nsc  |A snsnpioC 

The  input  power  to  the  Raman  cell  P»(0)  is  usually 
provided  by  hlackbody  or  spontaneous  Raman  noise. 
In  our  case  the  spontaneous  noise  is  dominant  so  that 
the  noise  in  a single  polarization  and  single  spatial  mode 
within  the  Raman  linewidth  Aep  is  Ps(0)  = lu>sAi>x  ~ 
10"12  W."  Thus,  for  a CO>  pump  power  of  10  MW,  the 
net  gain  required  to  amplify  the  spontaneous  power  up 
to  the  order  of  the  pump  power  is g„l  = ln|P*(i)//J*(0)) 
= ln[107 W/IO-'-  W)  = 44. 

The  16.95-pm  radiation  generated  by  four-wave 
mixing  can  also  he  amplified  by  the  stimulated  Raman 
gain  process.  It  has  been  previously  shown7''’*"  and  is 
evident  from  Eqs.  (la)  (Id)  that  the  four-wave  mixing 
conversion  efficiency  in  the  absence  of  significant 
Raman  amplification  and  pump  depletion  is  given  by 


4s«/2a)Supl2  5i2e^,  (3) 

/,  Vw»/  n„n, 


where  /„  /,„  /„  and  /,,  are  the  intensities  and  «„  n,„  n„, 
and  np  are  the  refractive  indices  at  the  four  fields.  In 
the  present  case  the  Nd:Y  AG  laser  source  at  a:,,  converts 
40%  of  its  energy  to  the  Stokes  wave  by  stimulated 
Raman  scattering  on  the  first  transit  of  the  25-pass  cell. 
The  incident  ('()••  laser  intensity  at  w,  then  generates 
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output  at  u >„  by  four-wave  mixing.  The  gener- 
ated output  power  at  16.95  nm  is  given  by  /•’„  «= 
0.4(1. 10pm/16.95 /im)z  /',  ~ l.fiSX  U)-  ' l\.  A 10-MW 
C02  laser  produces  Hi. 8 kW  or  184  p.J  of  energy  in  an 
8-nsec  pulse.  In  practice,  less  energy  is  produced  be- 
cause of  the  imperfect  spatial  overlap  of  the  beams  and 
the  nonzero  phase-mismatch  factor  Sh The  gener- 
ated 16-pm  energy  is  then  amplified  on  the  remaining 
passes  by  the  COr  laser-pumped  Raman  gain. 

'I'he  generated  four-wave  mixing  power  is  16  orders 
of  magnitude  greater  than  the  spontaneous  Raman 
power  and  thus  significantly  reduces  the  Raman  gain 
required  to  reach  threshold.10  The  required  net  Raman 
gain  with  four-wave  mixing  injection  is  gj  = ln[Ps(/)/ 
PJO)\  = ln(107W/10'W)  ~ 7.  The  required  Raman 
gain  and  thus  input  pump  power  is,  therefore,  reduced 
by  a factor  of  44/7  or  6.2  times  compared  to  the  ampli- 
fied spontaneous-emission  case. 

Figure  1 shows  a schematic  of  the  experimental  ap- 
paratus. The  CO2  laser  source  consists  of  a Lumonics 
102  oscillator  followed  by  a Lumonics  102  amplifier.  A 
Q-switched  unstable  resonator  Nd:YAG  laser  operating 
at  1.0G4  pm  provided  the  pumping  for  generating  the 
four-wave  mixing  signal.  The  Nd:YAG  laser  pulse 
energy  was  limited  to  70  mJ,  which  was  35  times  the 
measured  2-m.J  threshold  because  of  the  onset  of  dam- 
age of  the  metal  turning  mirror.  The  Nd:YAG  and  COr 
beams  were  combined  on  a ZnSe  plate  and  were  focused 
into  the  multipass  cell  by  a 9-m-radius  inode-matching 
mirror.  This  optical  arrangement  assured  that  the 
four-wave  mixing -generated  16.95-pm  radiation  had  the 
proper  spectral  and  spatial  profile  for  amplification  in 
the  remaining  passes  of  the  cell. 

The  multipass  cell  consists  of  a 3.77-m-long,  15.24- 
cm-diameter  stainless  steel  tube,  with  steerable 
12.70-cm-diameter,  2-m-radius  copper  mirrors  mounted 
at  each  end.  Light  is  coupled  into  and  out  of  the  cell 
through  salt  windows  mounted  on  the  sides  of  the  tube. 
Light  entering  the  cell  is  directed  down  the  axis  of  the 
tube  by  a 1 ,27-cm-diameter  steos-ng  mirror.  The  beam 
then  reflects  back  and  forth  between  the  12.70-cm- 
diameter  mirrors  refocusing  on  each  transit  and  walks 
in  a circle  around  the  rims  of  the  mirrors.  After  25 
transits,  the  beam  is  coupled  out  of  the  mult  ipass  cell 
by  a second  1.27-cm-diameter  steering  mirror  and  en- 
ters the  screen  room  where  spectral  and  energy  mea- 
surements arc  made.  A pyroelectric  energy  meter 
monitors  the  COr  energy.  The  16-pm  signal  is  sepa- 
rated from  the  10-pm  signal  by  a LiF  rcststrahlen  re- 
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Fig.  2.  Dependence  of  rotational  Raman  threshold  at  1.06 
pm  and  pressure-induced  absorption  at  10.6  pm  on  P-H2 
pressure. 


fleeter  followed  by  a dielectric  bandpass  filter  and  is  also 
monitored  by  a pyroelectric  energy  meter.  The  wave- 
length of  the  Raman  signal  measured  with  a grating 
monochrometer  is  16.95  pm,  as  expected. 

Figure  2 shows  the  rotational  Raman  threshold  de- 
pendence on  p- H2  pressure  at  1.06  pm  in  the  multipass 
cell.  As  the  pressure  increases,  the  threshold  rapidly 
drops  to  about  2 atm,  where  it  levels  off  at  about  2 mJ. 
In  a long  cell,  pressure-induced  absorption  in  Ho  gas  at 
10.6  pm  is  significant  and  limits  the  maximum  operat  ing 
pressure.  'I'he  second  curve  in  Fig.  2 is  a theoretical 
curve12  showing  the  fraction  of  CO?  energy  absorbed 
versus  Ho  pressure.  Measurements  of  CO2  transmit- 
tance versus  pressure  were  in  agreement  with  this  curve. 
Taking  both  pressure-induced  absorption  and  Nd:YAG 
threshold  data  into  account,  it  is  apparent  that  there  is 
an  optimum  operating  pressure  of  between  1 and  4 
atm. 

The  Nd:YAG  laser  linewidlh  was  also  varied  from 
single  axial  mode  to  about  0.5  cm- 1 by  inserting  etalons 
into  the  Nd:YAG  resonator.  The  Raman  threshold  was 
found  to  be  independent  of  laser  linewidth,  in  agree- 
ment with  theory.111-14  Our  COr  laser  source  was  thus 
operated  without  linewidth  control  to  obtain  the  max- 
imum output  pulse  energy.  Based  on  the  Nd:YAG 
threshold  measurements,  we  predicted  a Raman 
threshold  for  the  CO-2  laser  of  between  1 .7  and  2.6 .1  in 
a 70-nsec  pulse  for  3 atm  of  P-H2  in  the  multipass  cell. 
The  gas-breakdown  limit  in  our  present  E-ell  is  4.5-J 
input  CO2  energy. 

With  hot  h the  YAG  and  CO2  lasers  properly  aligned 
and  synchronized  in  the  IT  cell,  significant  16-pm 
energies  were  generated.  Figure  3 shows  the  16-pm 
output  energy  versus  CO2  laser  peak  input  power.  'I'he 
±50-nsec  t iming  jitter  between  the  YAG  and  (’Or  lasers 
resulted  in  large  pulse-to-pulse  variations  in  16-pm 
energy.  What  is  displayed  are  clusters  of  points  that 
represent  the  best  results  at  each  CO-2  laser-power 
setting.  'I’he  dots  represent  data  taken  at  1.5  atm,  and 
the  X’s  represent  data  taken  at  3-atin  Hr  pressure  in  tin- 
cell.  An  important  point  to  notice  is  that  the  16-pm 
energy  is  displayed  on  a log  scale.  An  increase  of  3 in 
COr  laser  power  results  in  an  increase  of  3 orders  of 
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PEAK  C02  LASER  POWER 

Kip.  3.  Output  energy  at  10.95  rim  versus  COg-laser  peak 
input  power. 

magnitude  in  16-pm  energy.  In  a four-wave  mixing 
process,  the  16-pm  energy  should  be  limited  to  the 
100  p-J  level  and  increase  linearly  with  CO^  power.  The 
exponential  increase  in  16-pm  energy  is  due  to  COi- 
pumped  Haman  gain.  The  slopes  of  the  lines  give  the 
Raman  exponential  gain  coefficient  for  the  multipass 
cell,  which  is  0.28/MW.  This  gain  coefficient  is  a factor 
of  3 less  than  the  plane-wave  gain  value  predicted  from 
Nd:YACi  laser  threshold  scaling  and  is  accounted  for  by 
the  spatial  mode  overlap  factor  of  the  10.6-  and  16.9-pm 
beams  proportional  to  (A,  + \p)/\l>. 

The  maximum  observed  16-pm  pulse  energies  were 
50  mJ.  The  8-nsec-long  16-pm  pulse  generated  by 
four-wave  mixing  was  significantly  shorter  than  the 
70-nsec  COg-laser  pulse,  and,  therefore,  depleted  a small 
portion  of  the  energy  in  the  CO-/  pulse.  The  50-mJ 
pulse  energy  corresponded  to  a 27%  peak  power  deple- 
tion or  to  a 40%  peak  photon  conversion  from  10.6  to 
16.9  pm. 

We  have  successfully  used  four-wave  mixing  injection 
to  reduce  the  COg  power  required  to  reach  stimulated 
rotational  Haman  threshold  on  the  S(U)  transition  in 
p- H i gas.  This  technique  should  permit  wider  appli- 
cation of  stimulated  Haman  scattering  processes  in  the 
infrared  as  an  efficient  means  of  frequency  conversion. 
'The  use  of  a tunable  high-pressure  COj-laser  source 
should  permit  continuous  tuning  over  wide  regions  of 
the  infrared  extending  from  8 pm  to  beyond  20  pm  by 
the  generation  of  Stokes  and  anti-Stokes  output  from 
various  rotational  levels  of  hydrogen  and  its  isotopes. 

The  high  conversion  efficiency  and  direct  scalability 
to  high  peak  and  average  powers  of  CO_>-pumped 
stimulated  rotational  Haman  scattering  should  permit 
this  approach  to  meet  the  large-scale  isotope  enrich- 
ment source  requirements. 
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